def test_get_root_mean_square_symmetric_surface_distance(self):
# same image
s = Surface(self._imaged1, self._imaged1)
self.assertEqual(s.get_root_mean_square_symmetric_surface_distance(),
0.)
# similar image
s = Surface(self._imaged1, self._imaged2)
self.assertAlmostEqual(
s.get_root_mean_square_symmetric_surface_distance(), 1.732050808)
# shifted image
s = Surface(self._imaged1, self._imaged1, (1, 1, 1), (0, 0, 0),
(1, 1, 1))
self.assertAlmostEqual(
s.get_root_mean_square_symmetric_surface_distance(), 1.732050808)
# shifte image \w non-one physical pixel spacing
s = Surface(self._imaged1, self._imaged1, (2, 2, 2), (0, 0, 0),
(1, 1, 1))
self.assertAlmostEqual(
s.get_root_mean_square_symmetric_surface_distance(), 3.464101615)
# different image
s = Surface(self._imaged1, self._imaged4)
self.assertAlmostEqual(
s.get_root_mean_square_symmetric_surface_distance(), 2.898275349)
# cube images A->B
s = Surface(self._imagedA, self._imagedB)
self.assertAlmostEqual(
s.get_root_mean_square_symmetric_surface_distance(), 1.4272480643)
# cube images B->A
s = Surface(self._imagedB, self._imagedA)
self.assertAlmostEqual(
s.get_root_mean_square_symmetric_surface_distance(), 1.4272480643)
def test_get_maximum_symmetric_surface_distance(self):
# same image
s = Surface(self._imaged1, self._imaged1)
self.assertEqual(s.get_maximum_symmetric_surface_distance(), 0.)
# similar image
s = Surface(self._imaged1, self._imaged2)
self.assertAlmostEqual(s.get_maximum_symmetric_surface_distance(),
1.732050808)
# shifted image
s = Surface(self._imaged1, self._imaged1, (1, 1, 1), (0, 0, 0),
(1, 1, 1))
self.assertAlmostEqual(s.get_maximum_symmetric_surface_distance(),
1.732050808)
# shifte image \w non-one physical pixel spacing
s = Surface(self._imaged1, self._imaged1, (2, 2, 2), (0, 0, 0),
(1, 1, 1))
self.assertAlmostEqual(s.get_maximum_symmetric_surface_distance(),
3.464101615)
# different image
s = Surface(self._imaged1, self._imaged4)
self.assertAlmostEqual(s.get_maximum_symmetric_surface_distance(),
5.196152423)
# cube images A->B
s = Surface(self._imagedA, self._imagedB)
self.assertAlmostEqual(s.get_maximum_symmetric_surface_distance(),
1.73205080757)
# cube images B->A
s = Surface(self._imagedB, self._imagedA)
self.assertAlmostEqual(s.get_maximum_symmetric_surface_distance(),
1.73205080757)
