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)