Draft01

01. import nibabel as nib
02. import numpy as np
03. 
04. fs_seg = nib.load(find_attachment_file('seg.nii')).get_data() == 53
05. print fs_seg.sum()
06. man_seg1 = nib.load(find_attachment_file('ROBERTO_LOTUFO_VBM_6min_SENSE_3_1y.nii.gz')).get_data() > 0
07. print man_seg1.shape
08. print man_seg1.sum()
09. 
10. man_seg2 = nib.load(find_attachment_file('ROBERTO_LOTUFO_VBM_6min_SENSE_11_1y.nii.gz')).get_data() > 0
11. print man_seg1.shape
12. print man_seg2.sum()
13. 
14. man_seg3 = nib.load(find_attachment_file('ROBERTO_LOTUFO_VBM_7min-32_SENSE_12_1y.nii.gz')).get_data() > 0
15. print man_seg1.shape
16. print man_seg3.sum()*0.45*0.45*0.45
3769
(240, 240, 180)
3105
(240, 240, 180)
3042
(240, 240, 180)
2824.875