def test_no_op(self):
vox_order = [dcmstack.closest_ortho_pat_axis(self.affine[:3, idx])[0] for idx in range(3)]
vox_order = "".join(vox_order)
vox_array, affine, perm = dcmstack.reorder_voxels(self.vox_array, self.affine, vox_order)
ok_((vox_array == self.vox_array).all())
ok_((affine == self.affine).all())
eq_(perm, (0, 1, 2))
eq_(np.may_share_memory(affine, self.affine), False)
def test_no_op(self):
vox_order = [
dcmstack.closest_ortho_pat_axis(self.affine[:3, idx])[0]
for idx in range(3)
]
vox_order = ''.join(vox_order)
vox_array, affine, perm = dcmstack.reorder_voxels(
self.vox_array, self.affine, vox_order)
ok_((vox_array == self.vox_array).all())
ok_((affine == self.affine).all())
eq_(perm, (0, 1, 2))
eq_(np.may_share_memory(affine, self.affine), False)
def test_closest_ortho_pat_axis():
eq_(dcmstack.closest_ortho_pat_axis((0.9, 0.1, 0.1)), 'lr')
eq_(dcmstack.closest_ortho_pat_axis((-0.9, 0.1, 0.1)), 'rl')
eq_(dcmstack.closest_ortho_pat_axis((0.1, 0.9, 0.1)), 'pa')
eq_(dcmstack.closest_ortho_pat_axis((0.1, -0.9, 0.1)), 'ap')
eq_(dcmstack.closest_ortho_pat_axis((0.1, 0.1, 0.9)), 'is')
eq_(dcmstack.closest_ortho_pat_axis((0.1, 0.1, -0.9)), 'si')
def test_closest_ortho_pat_axis():
eq_(dcmstack.closest_ortho_pat_axis((0.9, 0.1, 0.1)), "lr")
eq_(dcmstack.closest_ortho_pat_axis((-0.9, 0.1, 0.1)), "rl")
eq_(dcmstack.closest_ortho_pat_axis((0.1, 0.9, 0.1)), "pa")
eq_(dcmstack.closest_ortho_pat_axis((0.1, -0.9, 0.1)), "ap")
eq_(dcmstack.closest_ortho_pat_axis((0.1, 0.1, 0.9)), "is")
eq_(dcmstack.closest_ortho_pat_axis((0.1, 0.1, -0.9)), "si")
def test_closest_ortho_pat_axis():
eq_(dcmstack.closest_ortho_pat_axis((0.9, 0.1, 0.1)), 'lr')
eq_(dcmstack.closest_ortho_pat_axis((-0.9, 0.1, 0.1)), 'rl')
eq_(dcmstack.closest_ortho_pat_axis((0.1, 0.9, 0.1)), 'pa')
eq_(dcmstack.closest_ortho_pat_axis((0.1, -0.9, 0.1)), 'ap')
eq_(dcmstack.closest_ortho_pat_axis((0.1, 0.1, 0.9)), 'is')
eq_(dcmstack.closest_ortho_pat_axis((0.1, 0.1, -0.9)), 'si')
