Function SetDirection

Wrapping/Python/mitk/GeometryHelpers.cpp:97–115

}

void SetDirection(mitk::TimeGeometry* tg, py::array_t<double, py::array::c_style | py::array::forcecast> direction, mitk::TimeStepType t)
{
  if (direction.ndim() != 2 || direction.shape(0) != 3 || direction.shape(1) != 3)
    throw py::value_error("direction must be a 3x3 matrix");

  auto geom = GetGeometryForTimeStep(tg, t);
  auto buf = direction.unchecked<2>();
  const auto spacing = geom->GetSpacing();

  auto* transform = geom->GetIndexToWorldTransform();
  auto matrix = transform->GetMatrix();

  for (int row = 0; row < 3; ++row)
    for (int col = 0; col < 3; ++col)
      matrix[row][col] = buf(row, col) * spacing[col];

  transform->SetMatrix(matrix);
  geom->Modified();
}

std::array<double, 9> GetDirectionCosines(const mitk::TimeGeometry* tg, mitk::TimeStepType t)
{