tests/testthat/test-write_fs_weight.R

In freesurferformats: Read and Write 'FreeSurfer' Neuroimaging File Formats

test_that("One can write and re-read binary weight data", {
  vertex_indices = as.integer(c(1L, 45L, 255L));
  values = c(-0.56, 0.34, 44000.45);


  tmp_file = tempfile(fileext=".w");
  format_written = write.fs.weight(tmp_file, vertex_indices, values);


  res = read.fs.weight(tmp_file);
  expect_equal(res$vertex_indices, vertex_indices);
  expect_equal(res$values, values, tolerance=1e-2);
})


test_that("One can write and re-read ASCII weight data", {
  vertex_indices = as.integer(c(1L, 45L, 255L));
  values = c(-0.56, 0.34, 44000.45);


  tmp_file = tempfile(fileext=".w.asc");
  format_written = write.fs.weight(tmp_file, vertex_indices, values, format = 'asc');

  res = read.fs.weight(tmp_file);
  expect_equal(res$vertex_indices, vertex_indices);
  expect_equal(res$values, values, tolerance=1e-2);
})


test_that("Invalid arguments to write.fs.weight lead to errors", {
  vertex_indices = as.integer(c(1L, 45L, 255L));
  values = c(-0.56, 0.34, 44000.45);
  tmp_file = tempfile(fileext=".w");

  expect_error(write.fs.weight(tmp_file, (vertex_indices - 1L), values)); # vertex index 0 not allowed
  expect_error(write.fs.weight(tmp_file, (vertex_indices - 1L), values, format = 'asc')); # vertex index 0 not allowe

  expect_error(write.fs.weight(tmp_file, c(2, 24), values)); # 3 values but only 2 vertex indices
  expect_error(write.fs.weight(tmp_file, c(2, 24), values, format = 'asc')); # 3 values but only 2 vertex indices

  expect_error(write.fs.weight(tmp_file, vertex_indices, values, format = "no such format")); # invalid format
})