In [1]:
import numpy as np
# Set numpy to print only 2 decimal digits for neatness
np.set_printoptions(precision=2, suppress=True)
import os
import nibabel as nb
from nibabel.testing import data_path
In [2]:
from nilearn import image
In [3]:
nib.save(img1c, '/home/sophie/Documents/100412ss2onc.nii.gz')
In [ ]:
filename1 = '/home/sophie/Documents/100411ss2on.nii.gz'
img1 = nb.load(filename1)
data = img1.get_data()
S=data.shape
In [ ]:
S
In [ ]:
MeanImage=np.sum(data(:,:,:,1),axis=2)
In [ ]:
for x in range(1,100):
for y in range(1,200):
F(x,y)=np.sum(MeanImage(x:x+180,y:y+140))
[Id,M]=max(F)