atlas_stattk.cal_between_subj_relation('dice', atlastk.between_dice)
#########################################################################
dist_func = lambda x, y:atlastk.dist(x, y, atlas_deposite.img_resolution)
atlas_stattk.cal_multi_roi_relation('peak_distance', dist_func, atlastk.RoiStatTk.get_subject_peak)
hausdorff_func = lambda x, y: atlastk.hausdorff(x, y, atlas_deposite.img_resolution)
atlas_stattk.cal_multi_roi_relation('hausdorff_distance', hausdorff_func, atlastk.RoiStatTk.get_subject_label)
t4 = time.time() - t3
print 'atlas object construction done: %f sec' % t4
#########################################################################
atlas_stattk.cal_prob_maps()
atlas_stattk.save_prob_maps(probmapdir)
atlas_stattk.cal_max_prob_map(0, mapno=0, save=True, filename=probmapdir+'maxprob_thr0.nii.gz')
#########################################################################
mat_dir = '/nfs/j3/userhome/huangtaicheng/workingdir/parcellation_MT/BAA/results/yang_test/sub/data'
atlastk.output_mats(atlas_stattk, mat_dir, 'MT', MT_con, MT_con_key,'motion-fix' , 5.0, 0, 160)
atlas_stattk.cal_anat_relations('overlap', atlastk.multi_left_overlap, 'HOcort_thr0', anat_img)
#atlas_stattk.cal_anat_relations('jaccard', atlastk.multi_jaccard, 'HOcombo', anat_img)
#atlas_stattk.cal_anat_relations('dice', atlastk.multi_dice, 'HOcombo', anat_img)
#multi_hausdorff_func = lambda label, anat, aids, nid: atlastk.multi_hausdorff(label, anat, aids, nid, atlas_deposite.img_resolution)
#atlas_stattk.cal_anat_relations('dice', multi_hausdorff_func, 'HOcombo', anat_img)
##--------------------------------- Relationship between different subject (for each ROI) ------------------------------
atlas_stattk.cal_between_subj_relation('dice', atlastk.between_dice)
##--------------------------------- Relationship between different ROIs (for each subject)------------------------------
dist_func = lambda x, y:atlastk.dist(x, y, atlas_deposite.img_resolution)
atlas_stattk.cal_multi_roi_relation('peak_distance', dist_func, atlastk.RoiStatTk.get_subject_peak)
hausdorff_func = lambda x, y: atlastk.hausdorff(x, y, atlas_deposite.img_resolution)
atlas_stattk.cal_multi_roi_relation('hausdorff_distance', hausdorff_func, atlastk.RoiStatTk.get_subject_label)
t4 = time.time() - t3
print 'ROI relationship done: %f sec' % t4
#--------------------------------- Probmaps and MaxProbMap -------------------------------------------------------
atlas_stattk.cal_prob_maps()
atlas_stattk.save_prob_maps(probmapdir)
atlas_stattk.cal_max_prob_map(0, mapno=0, save=True, filename=probmapdir+'maxprob_thr0.nii.gz')
#---------------------------------- Output -------------------------------------------------------
mat_dir = '/nfs/t3/workingshop/yangzetian/Atlas/2006ObjectAtlas/data/mats'
atlastk.output_mats(atlas_stattk, mat_dir)
rint 'Total time: %f sec' % (time.tim:e()-t0)
