def test_present_options_2():
data = np.arange(100000)
data = data.reshape((100,100,10))
full=present_3d_options(data,axis=2)
first=np.ceil(np.sqrt(10))
second=np.ceil(10/first)
assert full.shape == (100*first,100*second)
from Image_Visualizing import present_3d, make_mask, present_3d_options
i="sub001"
brain=nib.load(path_to_data+i+'/BOLD/task001_run001/bold.nii.gz')
data=brain.get_data()
brain_hi= nib.load(path_to_data + i+'/anatomy/inplane001_brain.nii.gz')
joy_hi=brain_hi.get_data()
joy=data[...,7]
plt.imshow(present_3d(data[...,7]),cmap="gray",interpolation="nearest")
plt.figure()
plt.imshow(present_3d_options(joy,2),cmap="gray",interpolation="nearest")
plt.title("2")
plt.figure()
plt.imshow(present_3d_options(joy,1),cmap="gray",interpolation="nearest")
plt.title("1")
plt.figure()
plt.imshow(present_3d_options(joy,0),cmap="gray",interpolation="nearest")
plt.title("0")
plt.close()
plt.imshow(present_3d_options(joy_hi,2),cmap="gray",interpolation="nearest")
plt.title("2")
plt.savefig(location_of_images+"kent_brain_2.png")
sys.path.append(location_of_functions)
from Image_Visualizing import present_3d, make_mask, present_3d_options
i = "sub001"
brain = nib.load(path_to_data + i + '/BOLD/task001_run001/bold.nii.gz')
data = brain.get_data()
brain_hi = nib.load(path_to_data + i + '/anatomy/inplane001_brain.nii.gz')
joy_hi = brain_hi.get_data()
joy = data[..., 7]
plt.imshow(present_3d(data[..., 7]), cmap="gray", interpolation="nearest")
plt.figure()
plt.imshow(present_3d_options(joy, 2), cmap="gray", interpolation="nearest")
plt.title("2")
plt.figure()
plt.imshow(present_3d_options(joy, 1), cmap="gray", interpolation="nearest")
plt.title("1")
plt.figure()
plt.imshow(present_3d_options(joy, 0), cmap="gray", interpolation="nearest")
plt.title("0")
plt.close()
plt.imshow(present_3d_options(joy_hi, 2), cmap="gray", interpolation="nearest")
plt.title("2")
plt.savefig(location_of_images + "kent_brain_2.png")
plt.close()
plt.imshow(present_3d_options(joy_hi, 1), cmap="gray", interpolation="nearest")
def test_present_options_0():
data = np.arange(100000)
data = data.reshape((100,100,10))
full=present_3d_options(data,axis=0)
assert full.shape == (10*10,100*10)
