def test_blot_with_default(self):
"""
Test do_blot with default grid image
"""
input = os.path.basename(
self.get_input_file('input', 'j8bt06nyq_flt.fits'))
output = 'output_blot_default.fits'
output_difference = os.path.basename(
self.get_data('truth', 'difference_blot_default.txt'))
output_template = os.path.basename(
self.get_data('truth', 'reference_blot_default.fits'))
insci = self.read_image(input)
insci = self.make_grid_image(insci, 64, 100.0)
input_wcs = self.read_wcs(input)
output_wcs = self.read_wcs(output_template)
expin = 1.0
outsci = ablot.do_blot(insci,
input_wcs,
output_wcs,
expin,
coeffs=False)
self.write_image(output, output_wcs, outsci)
template_data = self.read_image(output_template)
(min_diff, med_diff,
max_diff) = self.centroid_statistics("blot with defaults",
output_difference, outsci,
template_data, 20.0, 16)
assert (med_diff < 1.0e-6)
assert (max_diff < 1.0e-5)
def test_blot_with_lan5(self):
"""
Test do_blot with lan5 grid image
"""
input = os.path.basename(self.get_input_file('input', 'j8bt06nyq_unit.fits'))
output = 'output_blot_lan5.fits'
output_difference = os.path.basename(self.get_data('truth',
'difference_blot_lan5.txt'))
output_template = os.path.basename(self.get_data('truth',
'reference_blot_lan5.fits'))
insci = self.read_image(input)
insci = self.make_grid_image(insci, 64, 100.0)
input_wcs = self.read_wcs(input)
output_wcs = self.read_wcs(output_template)
expin = 1.0
outsci = ablot.do_blot(insci, input_wcs, output_wcs, expin, interp="lan5", coeffs = False)
self.write_image(output, output_wcs, outsci)
template_data = self.read_image(output_template)
(min_diff, med_diff, max_diff) = self.centroid_statistics("blot with lan5", output_difference,
outsci, template_data, 20.0, 16)
assert(med_diff < 1.0e-6)
assert(max_diff < 1.0e-5)
def test_blot_with_point(self):
"""
Test do_blot with point image
"""
input = os.path.basename(self.get_input_file('input', 'j8bt06nyq_unit.fits'))
output = 'output_blot_point.fits'
output_difference = os.path.basename(self.get_data('truth',
'difference_blot_point.txt'))
output_template = os.path.basename(self.get_data('truth',
'reference_blot_point.fits'))
insci = self.read_image(input)
input_wcs = self.read_wcs(input)
insci = self.make_point_image(insci, (500, 200), 40.0)
output_wcs = self.read_wcs(output_template)
expin = 1.0
outsci = ablot.do_blot(insci, input_wcs, output_wcs, expin, coeffs = False)
output_bounds = self.bound_image(outsci)
self.write_image(output, output_wcs, outsci)
template_data = self.read_image(output_template)
template_bounds = self.bound_image(template_data)
(min_diff, med_diff, max_diff) = self.centroid_statistics("blot with point", output_difference,
outsci, template_data, 20.0, 16)
assert(med_diff < 1.0e-6)
assert(max_diff < 1.0e-5)
def test_blot_with_image(self):
"""
Test do_blot with full image
"""
input = os.path.basename(self.get_input_file('input', 'j8bt06nyq_unit.fits'))
output = 'output_blot_image.fits'
output_template = os.path.basename(self.get_data('truth',
'reference_blot_image.fits'))
insci = self.read_image(input)
input_wcs = self.read_wcs(input)
output_wcs = self.read_wcs(output_template)
expin = 1.0
outsci = ablot.do_blot(insci, input_wcs, output_wcs, expin, coeffs = False)
self.write_image(output, output_wcs, outsci)
self.compare_outputs([(output, output_template)])
def test_blot_with_image(self):
"""
Test do_blot with full image
"""
input = os.path.basename(
self.get_input_file('input', 'j8bt06nyq_flt.fits'))
output = 'output_blot_image.fits'
output_template = os.path.basename(
self.get_data('truth', 'reference_blot_image.fits'))
insci = self.read_image(input)
input_wcs = self.read_wcs(input)
output_wcs = self.read_wcs(output_template)
expin = 1.0
outsci = ablot.do_blot(insci,
input_wcs,
output_wcs,
expin,
coeffs=False)
self.write_image(output, output_wcs, outsci)
self.compare_outputs([(output, output_template)])
def test_blot_with_point(self):
"""
Test do_blot with point image
"""
input = os.path.basename(
self.get_input_file('input', 'j8bt06nyq_unit.fits'))
output = 'output_blot_point.fits'
output_difference = os.path.basename(
self.get_data('truth', 'difference_blot_point.txt'))
output_template = os.path.basename(
self.get_data('truth', 'reference_blot_point.fits'))
insci = self.read_image(input)
input_wcs = self.read_wcs(input)
insci = self.make_point_image(insci, (500, 200), 40.0)
output_wcs = self.read_wcs(output_template)
expin = 1.0
outsci = ablot.do_blot(insci,
input_wcs,
output_wcs,
expin,
coeffs=False)
output_bounds = self.bound_image(outsci)
self.write_image(output, output_wcs, outsci)
template_data = self.read_image(output_template)
template_bounds = self.bound_image(template_data)
(min_diff, med_diff,
max_diff) = self.centroid_statistics("blot with point",
output_difference, outsci,
template_data, 20.0, 16)
assert (med_diff < 1.0e-6)
assert (max_diff < 1.0e-5)
