def visualize_axial(niidir, savedir, min=0.75, max=1.0): gl.viewaxial(1) for method in ['grad', 'cam']: for network in ['Vis', 'SomMot', 'DorsAttn', 'SalVentAttn', 'Limbic', 'Cont', 'Default']: for gender in ['female', 'male']: gl.overlayload(os.path.join(niidir, 'network', 'saliency_{}_{}_top20_{}_lh'.format(method, gender, network))) gl.overlayload(os.path.join(niidir, 'network', 'saliency_{}_{}_top20_{}_rh'.format(method, gender, network))) gl.minmax(1, min, max) gl.minmax(2, min, max) gl.colorname(1, colorname) gl.colorname(2, colorname) gl.savebmp(os.path.join(savedir, '{}_{}_{}_axial.png'.format(method, network, gender))) gl.overlaycloseall()
def visualize_sagittal(niidir, savedir, min=0.75, max=1.0): for hemisphere in ['lh', 'rh']: if hemisphere =='lh': gl.clipazimuthelevation(0.49, 90, 0) elif hemisphere =='rh': gl.clipazimuthelevation(0.49, 270, 0) for network in ['Vis', 'SomMot', 'DorsAttn', 'SalVentAttn', 'Limbic', 'Cont', 'Default']: for gender in ['female', 'male']: gl.overlayload(os.path.join(niidir, 'network', 'saliency_{}_top20_{}_{}'.format(method, gender, network, hemisphere))) gl.minmax(1, min, max) gl.colorname(1, colorname) gl.viewsagittal(1) gl.savebmp(os.path.join(savedir, '{}_{}_sagittal_{}_lt.png'.format(network, gender, hemisphere))) gl.viewsagittal(0) gl.savebmp(os.path.join(savedir, '{}_{}_sagittal_{}_rt.png'.format(network, gender, hemisphere))) gl.overlaycloseall()
# %% overlay ---------------------- FILESPATH='D:\\GD_UNICAMP\\IC_NeuroFisica\\Projetos\\fMRI_TLE_2020\\Seed_based_analysis' os.chdir(FILESPATH) networks_names=['DMN','SMOTOR','VISUAL','AUDIT','ECN_L','ECN_R'] netfile_sufx= ['_Group1X','_Group2L','_Group3R','_Group4N'] for net in networks_names: savefolder=FILESPATH + '\\' + net + '\\' for i in netfile_sufx: FILENAME=net+i SEEDFILENAME=net+"_seed" gl.overlaycloseall() gl.overlayload(FILENAME) gl.overlayload(SEEDFILENAME) gl.colorname (2,"3blue") #gl.opacity(1,50) # changing view ---------------------- gl.azimuthelevation(140,20) ktime= 1 ksteps= 72 for x in range(0, ksteps): gl.azimuthelevation(140+(x*5),20) gl.wait(ktime) if x < 10: