yeo_DMN_dist.py

#! /usr/bin/env python
from hcp_class import hcp_subj
import numpy as np
import gdist
import surfdist.analysis
import surfdist.utils
import surfdist
import nibabel as nib
import pickle
import sys
import os

subj=sys.argv[1]
out=f'/well/margulies/projects/pkReliability/Dist2SensoryBorder/{subj}/'
os.makedirs(out,exist_ok=True)

subj=hcp_subj(subj,4)
print(subj.subj)    

Lsrf=(subj.Lcoords,subj.Lfaces)
Rsrf=(subj.Rcoords,subj.Rfaces)

Lyeo=np.load('./32k_fs_LR/yeo/yeo.L.DMN.npy')
Lyeo=np.where(Lyeo==1)[0]
Lyeo=surfdist.analysis.dist_calc(Lsrf,subj.Lfill,Lyeo)
np.save(f'{out}/L.yeo.npy',Lyeo)


Lgrp=np.load('./32k_fs_LR/top90GrpGrad.L.32k.npy')
Lgrp=np.where(Lgrp==1)[0]
Lgrp=surfdist.analysis.dist_calc(Lsrf,subj.Lfill,Lgrp)
np.save(f'{out}/L.grpGrad90.npy',Lgrp)


Ryeo=np.load('./32k_fs_LR/yeo/yeo.R.DMN.npy')
Ryeo=np.where(Ryeo==1)[0]
Ryeo=surfdist.analysis.dist_calc(Rsrf,subj.Rfill,Ryeo)
np.save(f'{out}/R.yeo.npy',Ryeo)


Rgrp=np.load('./32k_fs_LR/top90GrpGrad.R.32k.npy')
Rgrp=np.where(Rgrp==1)[0]
Rgrp=surfdist.analysis.dist_calc(Rsrf,subj.Rfill,Rgrp)
np.save(f'{out}/R.grpGrad90.npy',Rgrp)