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
