def test_convert():
td = get_testdir()
s = Surface(op.join(td, 'in.surf.gii'))
s.save('test.vtk')
s2 = Surface('test.vtk')
assert np.allclose(s.points, s2.points)
os.remove('test.vtk')
def convert_surface():
"""
CLI for converting surface formats
"""
parser = CommonParser(
'surf',
'coords',
'struct',
'out',
description="""Convert a surface file (.white/.pial/.vtk/.surf.gii).
NB FreeSurfer files will have the c_ras offset automatically
applied during conversion.""")
parsed = parser.parse_args()
if parsed.coords == 'fsl' and parsed.struct:
insurf = Surface(parsed.surf, 'fsl', parsed.struct)
else:
insurf = Surface(parsed.surf)
insurf.save(parsed.out)