def test_find_sh_hole_center(self):
"""
Test FindCircleCenter for holes
"""
# Real image from the DELPHI
data = tiff.read_data(os.path.join(TEST_IMAGE_PATH, "sh_hole_up.tiff"))
hole_coordinates = delphi.FindCircleCenter(data[0], delphi.HOLE_RADIUS, 6, darkest=True)
# FIXME: it fails (but not that important for calibration)
expected_coordinates = (-0.00014212, 9.405e-05) # about: 888, 934 = -0.00014212, 9.405e-05
numpy.testing.assert_almost_equal(hole_coordinates, expected_coordinates)
def test_find_hole_center(self):
"""
Test FindCircleCenter for holes
"""
data = hdf5.read_data("sem_hole.h5")
C, T, Z, Y, X = data[0].shape
data[0].shape = Y, X
hole_coordinates = delphi.FindCircleCenter(data[0], 0.02032, 3)
expected_coordinates = (390.5, 258.5)
numpy.testing.assert_almost_equal(hole_coordinates,
expected_coordinates)
def test_find_lens_center(self):
"""
Test FindCircleCenter for lenses
"""
data = hdf5.read_data("navcam-calib2.h5")
Z, Y, X = data[0].shape
lens_coordinates = delphi.FindCircleCenter(data[0][0],
delphi.LENS_RADIUS, 6)
expected_coordinates = (450.5, 445.5)
numpy.testing.assert_almost_equal(lens_coordinates,
expected_coordinates)
def test_find_hole_center(self):
"""
Test FindCircleCenter for holes
"""
# Note: this hole image has a better contrast and less noise than typical
# hole images on the DELPHI
data = hdf5.read_data(os.path.join(TEST_IMAGE_PATH, "sem_hole.h5"))
C, T, Z, Y, X = data[0].shape
data[0].shape = Y, X
hole_coordinates = delphi.FindCircleCenter(data[0], 0.02032, 6, darkest=True)
expected_coordinates = (0.0052705, -0.0018415) # (391.5, 257.5) px
numpy.testing.assert_almost_equal(hole_coordinates, expected_coordinates)
