Exemple #1
0
def test_labels1():
    img = load_image(funcfile)
    data = img.get_data()
    parcelmap = Image(img[0].get_data(), AfT('kji', 'zyx', np.eye(4)))
    parcelmap = (parcelmap.get_data() * 100).astype(np.int32)
    v = 0
    for i, d in axis0_generator(data, parcels(parcelmap)):
        v += d.shape[1]
    assert_equal(v, parcelmap.size)
Exemple #2
0
def test_labels1():
    img = load_image(funcfile)
    data = img.get_data()
    parcelmap = fromarray(np.asarray(img[0]), 'kji', 'zyx')
    parcelmap = (np.asarray(parcelmap) * 100).astype(np.int32)
    v = 0
    for i, d in axis0_generator(data, parcels(parcelmap)):
        v += d.shape[1]
    assert_equal(v, parcelmap.size)
Exemple #3
0
def test_labels1():
    img = load_image(funcfile)
    data = img.get_data()
    parcelmap = Image(img[0].get_data(), AfT("kji", "zyx", np.eye(4)))
    parcelmap = (parcelmap.get_data() * 100).astype(np.int32)
    v = 0
    for i, d in axis0_generator(data, parcels(parcelmap)):
        v += d.shape[1]
    assert_equal(v, parcelmap.size)
Exemple #4
0
def test_iter():
    img = load_image(funcfile)
    img_shape = img.shape
    exp_shape = (img_shape[0],) + img_shape[2:]
    j = 0
    for i, d in axis0_generator(img.get_data()):
        j += 1
        assert_equal(d.shape, exp_shape)
        del(i); gc.collect()
    assert_equal(j, img_shape[1])
Exemple #5
0
def test_iter():
    img = load_image(funcfile)
    img_shape = img.shape
    exp_shape = (img_shape[0],) + img_shape[2:]
    j = 0
    for i, d in axis0_generator(img.get_data()):
        j += 1
        assert_equal(d.shape, exp_shape)
        del(i); gc.collect()
    assert_equal(j, img_shape[1])
Exemple #6
0
def test_iterate_over_image():
    # Fit a model, iterating over the slices of an array
    # associated to an FmriImage.
    c = np.zeros(FDATA.shape[1:]) + 0.5
    res_gen = result_generator(flatten_generator(axis0_generator(FDATA)))
    write_data(c, unflatten_generator(contrast_generator(res_gen)))
    # Fit a model, iterating over the array associated to an
    # FmriImage, iterating over a list of ROIs defined by binary
    # regions of the same shape as a frame of FmriImage

    # this might really be an anatomical image or AR(1) coefficients
    a = np.asarray(FDATA[0])
    p = np.greater(a, a.mean())
    d = np.ones(FDATA.shape[1:]) * 2.0
    flat_gen = flatten_generator(axis0_generator(FDATA, parcels(p)))
    write_data(d, contrast_generator(result_generator(flat_gen)))
    assert_array_almost_equal(d, c)

    e = np.zeros(FDATA.shape[1:]) + 3.0
    flat_gen2 = flatten_generator(axis0_generator(FDATA, parcels(p)))
    write_data(e, f_generator(contrast, result_generator(flat_gen2)))
    assert_array_almost_equal(d, e)
Exemple #7
0
def test_iterate_over_image():
    # Fit a model, iterating over the slices of an array
    # associated to an FmriImage.
    c = np.zeros(FDATA.shape[1:]) + 0.5
    res_gen = result_generator(flatten_generator(axis0_generator(FDATA)))
    write_data(c, unflatten_generator(contrast_generator(res_gen)))
    # Fit a model, iterating over the array associated to an
    # FmriImage, iterating over a list of ROIs defined by binary
    # regions of the same shape as a frame of FmriImage

    # this might really be an anatomical image or AR(1) coefficients
    a = np.asarray(FDATA[0])
    p = np.greater(a, a.mean())
    d = np.ones(FDATA.shape[1:]) * 2.0
    flat_gen = flatten_generator(axis0_generator(FDATA, parcels(p)))
    write_data(d, contrast_generator(result_generator(flat_gen)))
    assert_array_almost_equal(d, c)

    e = np.zeros(FDATA.shape[1:]) + 3.0
    flat_gen2 = flatten_generator(axis0_generator(FDATA, parcels(p)))
    write_data(e, f_generator(contrast, result_generator(flat_gen2)))
    assert_array_almost_equal(d, e)