def test_eval_sliver_distance_two_points(self):
"""
Two points 2 by 2 pixels diagonal.
"""
import volumetry_evaluation as ve
vol1 = np.zeros([20, 21, 22], dtype=np.int8)
vol1[10, 10, 10] = 1
vol2 = np.zeros([20, 21, 22], dtype=np.int8)
vol2[10, 12, 12] = 1
eval1 = ve.distance_matrics(vol1, vol2, [1, 1, 1])
self.assertAlmostEquals(eval1[2], np.sqrt(2) * 2)
def test_eval_sliver_distance_two_points_with_half_voxelsize(self):
"""
Two points 2 by 2 pixels diagonal and voxelsize is 0.5
"""
import volumetry_evaluation as ve
vol1 = np.zeros([20, 21, 22], dtype=np.int8)
vol1[10, 10, 10] = 1
vol2 = np.zeros([20, 21, 22], dtype=np.int8)
vol2[10, 12, 12] = 1
eval1 = ve.distance_matrics(vol1, vol2, [0.5, 0.5, 0.5])
#print eval1
self.assertAlmostEquals(eval1[2], np.sqrt(2))
def test_eval_sliver_distance_for_two_pixels_bigger_volume(self):
"""
maxd is measured in corner. It is space diagonal of 2 pixels cube.
"""
import volumetry_evaluation as ve
vol1 = np.zeros([20, 21, 22], dtype=np.int8)
vol1[10:15, 10:15, 10:15] = 1
vol2 = np.zeros([20, 21, 22], dtype=np.int8)
vol2[8:17, 8:17, 8:17] = 1
eval1 = ve.distance_matrics(vol1, vol2, [1, 1, 1])
self.assertAlmostEquals(eval1[2], 3 ** (0.5) * 2)
