def test_aff_to_hdr(): hdr = {} affine = np.diag([1,2,3,1]) affine[:3,3] = [10,11,12] tv.aff_to_hdr(affine, hdr) yield assert_array_almost_equal(tv.aff_from_hdr(hdr), affine) # put flip into affine aff2 = affine.copy() aff2[:,2] *=-1 tv.aff_to_hdr(aff2, hdr) yield assert_array_almost_equal(tv.aff_from_hdr(hdr), aff2)
def test_aff_to_hdr(): for version in (1, 2): hdr = {'version': version} affine = np.diag([1, 2, 3, 1]) affine[:3, 3] = [10, 11, 12] tv.aff_to_hdr(affine, hdr) yield assert_array_almost_equal(tv.aff_from_hdr(hdr), affine) # put flip into affine aff2 = affine.copy() aff2[:, 2] *= -1 tv.aff_to_hdr(aff2, hdr) yield assert_array_almost_equal(tv.aff_from_hdr(hdr), aff2) if version == 2: yield assert_array_almost_equal(hdr['vox_to_ras'], aff2)
def test_aff_to_hdr(): for version in (1, 2): hdr = {'version': version} affine = np.diag([1,2,3,1]) affine[:3,3] = [10,11,12] tv.aff_to_hdr(affine, hdr) yield assert_array_almost_equal(tv.aff_from_hdr(hdr), affine) # put flip into affine aff2 = affine.copy() aff2[:,2] *=-1 tv.aff_to_hdr(aff2, hdr) yield assert_array_almost_equal(tv.aff_from_hdr(hdr), aff2) if version == 2: yield assert_array_almost_equal(hdr['vox_to_ras'], aff2)