def test_diff_ct():
d = np.zeros([5, 5, 5, 8])
for z in range(8):
d[..., z] = z
img = AslImage(name="asldata", image=d, tis=[1], iaf="ct", order='lrt')
img = img.diff()
assert img.ntis == 1
assert img.tis == [1]
assert not img.have_plds
assert img.rpts == [4]
assert img.ntc == 1
assert img.order == "rt"
data = img.nibImage.get_data()
assert list(data.shape) == [5, 5, 5, 4]
assert np.all(data == -1)
def test_multite_diff():
d = np.zeros([5, 5, 5, 16])
for z in range(16):
d[..., z] = z
img = AslImage(name="d",
image=d,
plds=[1, 2],
iaf="tc",
ibf="tis",
tes=[8, 9],
casl=True,
bolus=[3, 4])
imgdiff = img.diff()
assert imgdiff.ntes == 2
assert imgdiff.tes == [8, 9]
assert imgdiff.ntis == 2
assert imgdiff.have_plds
assert imgdiff.plds == [1, 2]
assert imgdiff.tis == [4, 6]
assert imgdiff.rpts == [2, 2]
assert imgdiff.ntc == 1
assert imgdiff.order == "ert"
assert imgdiff.nvols == 8
np.all(imgdiff.data == 2)