def test_io_surface(): """Test reading and writing of Freesurfer surface mesh files """ tempdir = _TempDir() fname_quad = op.join(data_path, 'subjects', 'bert', 'surf', 'lh.inflated.nofix') fname_tri = op.join(data_path, 'subjects', 'fsaverage', 'surf', 'lh.inflated') for fname in (fname_quad, fname_tri): pts, tri, vol_info = read_surface(fname, read_metadata=True) write_surface(op.join(tempdir, 'tmp'), pts, tri, volume_info=vol_info) c_pts, c_tri, c_vol_info = read_surface(op.join(tempdir, 'tmp'), read_metadata=True) assert_array_equal(pts, c_pts) assert_array_equal(tri, c_tri) assert_true(_is_equal_dict([vol_info, c_vol_info]))
def test_io_surface(): """Test reading and writing of Freesurfer surface mesh files.""" tempdir = _TempDir() fname_quad = op.join(data_path, 'subjects', 'bert', 'surf', 'lh.inflated.nofix') fname_tri = op.join(data_path, 'subjects', 'fsaverage', 'surf', 'lh.inflated') for fname in (fname_quad, fname_tri): with warnings.catch_warnings(record=True) as w: pts, tri, vol_info = read_surface(fname, read_metadata=True) assert_true(all('No volume info' in str(ww.message) for ww in w)) write_surface(op.join(tempdir, 'tmp'), pts, tri, volume_info=vol_info) with warnings.catch_warnings(record=True) as w: # No vol info c_pts, c_tri, c_vol_info = read_surface(op.join(tempdir, 'tmp'), read_metadata=True) assert_array_equal(pts, c_pts) assert_array_equal(tri, c_tri) assert_true(_is_equal_dict([vol_info, c_vol_info]))