elif args.im : tranfo = -mid_point mesh.set_vertices(mesh.vertices_translation(tranfo)) print "imz : " + str(tranfo) elif args.mv : tranfo = mid_point*z mesh.set_vertices(mesh.vertices_translation(tranfo)) print "mz : " + str(tranfo) elif args.imv : tranfo = -mid_point*z mesh.set_vertices(mesh.vertices_translation(tranfo)) print "imz : " + str(tranfo) if args.RAS2LPS : mesh.set_vertices(mesh.vertices_flip([-1, -1, 1])) mesh.update_polydata() mesh.save(surface_file_out) """ ### FreeSurfer info out1 = os.popen("mri_info --ras2vox-tkr " + file_name ).read() out2 = os.popen("mri_info --vox2ras " + file_name ).read() r2v = np.array([float(x) for x in out1.split()]).reshape(4,4) v2r = np.array([float(x) for x in out2.split()]).reshape(4,4) transfo = np.dot(r2v, v2r) print transfo