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)