# get the masks
mask_path = op.join(func_path, "%s_mask.gii" % side)
# read the texture and curvature data
tex = read(tex_path).darrays[0].data.ravel()
curv = (read(curv_path).darrays[0].data < 0).astype(np.int)
mask = read(mask_path).darrays[0].data
### Plot meshes
f = mlab.figure(bgcolor=(0.05, 0, 0.1), size=(400, 400))
mlab.clf()
vmin, vmax = -np.pi, np.pi
plot_retino_image(
inflated_mesh_path,
name="%s hemisphere" % side,
mask=mask,
tf=None,
tex=tex,
curv=curv,
vmin=vmin,
vmax=vmax,
)
if side == "left":
mlab.view(280, 120)
else:
mlab.view(250, 120)
mlab.savefig(op.join(func_path, "%s_%s.png") % (subject, side))
mlab.show()
tex_path = op.join(func_path, '%s_phase_wedge.gii' % side)
inflated_mesh_path = op.join(anat_path, '%sh.inflated.gii' % side[0])
curv_path = op.join(anat_path, '%sh.avg_curv.gii' % side[0])
# get the masks
mask_path = op.join(func_path, '%s_mask.gii' % side)
# read the texture and curvature data
tex = read(tex_path).darrays[0].data.ravel()
curv = (read(curv_path).darrays[0].data < 0).astype(np.int)
mask = read(mask_path).darrays[0].data
### Plot meshes
f = mlab.figure(bgcolor=(.05, 0, .1), size=(400, 400))
mlab.clf()
vmin, vmax = -np.pi, np.pi
plot_retino_image(inflated_mesh_path,
name="%s hemisphere" % side,
mask=mask,
tf=None,
tex=tex,
curv=curv,
vmin=vmin, vmax=vmax)
if side == 'left':
mlab.view(280, 120)
else:
mlab.view(250, 120)
mlab.savefig(op.join(func_path, '%s_%s.png') % (subject, side))
mlab.show()
ltex = read(ltex_path).darrays[0].data.ravel()
lcurv = (read(lcurv_path).darrays[0].data < 0).astype(np.int)
lmask = read(lmask_path).darrays[0].data
rtex = read(rtex_path).darrays[0].data.ravel()
rcurv = (read(rcurv_path).darrays[0].data < 0).astype(np.int)
rmask = read(rmask_path).darrays[0].data
### Plot meshes
# left hemisphere
f = mlab.figure(bgcolor=(.05, 0, .1), size=(400, 400))
mlab.clf()
plot_retino_image(lmesh_path_inflated,
name="LeftHemisphere",
mask=lmask,
tf=None,
tex=ltex,
curv=lcurv,
vmin=-np.pi,
vmax=np.pi)
#plot_retino_image(lmesh_path_inflated, name="LeftHemisphere", mask=lmask,
# tf=None, tex=ltex, curv=lcurv, vmin=-np.pi, vmax=0)
mlab.view(280, 120)
mlab.savefig(op.join(func_path, '%s_%s.png') % (subject, 'left'))
# right hemisphere
f = mlab.figure(bgcolor=(.05, 0, .1), size=(400, 400))
mlab.clf()
plot_retino_image(rmesh_path_inflated,
name="RightHemisphere",
mask=rmask,
tf=None,
rmask_path = op.join(func_path, 'right_mask.gii')
# read the texture and curvature data
ltex = read(ltex_path).darrays[0].data.ravel()
lcurv = (read(lcurv_path).darrays[0].data < 0).astype(np.int)
lmask = read(lmask_path).darrays[0].data
rtex = read(rtex_path).darrays[0].data.ravel()
rcurv = (read(rcurv_path).darrays[0].data < 0).astype(np.int)
rmask = read(rmask_path).darrays[0].data
### Plot meshes
# left hemisphere
f = mlab.figure(bgcolor=(.05, 0, .1), size=(400, 400))
mlab.clf()
plot_retino_image(lmesh_path_inflated, name="LeftHemisphere", mask=lmask,
tf=None, tex=ltex, curv=lcurv, vmin=-np.pi, vmax=np.pi)
#plot_retino_image(lmesh_path_inflated, name="LeftHemisphere", mask=lmask,
# tf=None, tex=ltex, curv=lcurv, vmin=-np.pi, vmax=0)
mlab.view(280, 120)
mlab.savefig(op.join(func_path,'%s_%s.png') % (subject, 'left'))
# right hemisphere
f = mlab.figure(bgcolor=(.05, 0, .1), size=(400, 400))
mlab.clf()
plot_retino_image(rmesh_path_inflated, name="RightHemisphere", mask=rmask,
tf=None, tex=rtex, curv=rcurv, vmin=-np.pi, vmax=np.pi)
#plot_retino_image(rmesh_path_inflated, name="RightHemisphere", mask=rmask,
# tf=None, tex=rtex, curv=rcurv, vmin=0, vmax=np.pi)
mlab.view(250, 120)
mlab.savefig(op.join(func_path, '%s_%s.png') % (subject, 'right'))
mlab.show()