Method GetRotatedId

Filters/General/vtkAxisAlignedTransformFilter.cxx:660–694

----------------------------------------------------------------------------

//----------------------------------------------------------------------------
int vtkAxisAlignedTransformFilter::GetRotatedId(
  int id, int rotationMatrix[3][3], int newDims[3], int dims[3], int Tvec[3], bool transposed)
{
  int x, y, z;
  if (transposed)
  {
    z = id % std::max(newDims[2] - 1, 1);
    y = (id / (std::max(newDims[2] - 1, 1))) % std::max(newDims[1] - 1, 1);
    x = id / std::max(newDims[2] - 1, 1) / std::max(newDims[1] - 1, 1);
  }
  else
  {
    x = id % std::max(newDims[0] - 1, 1);
    y = (id / (std::max(newDims[0] - 1, 1))) % std::max(newDims[1] - 1, 1);
    z = id / std::max(newDims[0] - 1, 1) / std::max(newDims[1] - 1, 1);
  }

  int newX =
    rotationMatrix[0][0] * x + rotationMatrix[1][0] * y + rotationMatrix[2][0] * z + Tvec[0];
  int newY =
    rotationMatrix[0][1] * x + rotationMatrix[1][1] * y + rotationMatrix[2][1] * z + Tvec[1];
  int newZ =
    rotationMatrix[0][2] * x + rotationMatrix[1][2] * y + rotationMatrix[2][2] * z + Tvec[2];

  if (transposed)
  {
    return newX * std::max(dims[2] - 1, 1) * std::max(dims[1] - 1, 1) +
      newY * std::max(dims[2] - 1, 1) + newZ;
  }
  else
  {
    return newZ * std::max(dims[0] - 1, 1) * std::max(dims[1] - 1, 1) +
      newY * std::max(dims[0] - 1, 1) + newX;
  }
}

//----------------------------------------------------------------------------
vtkAxisAlignedTransformFilter::Axis vtkAxisAlignedTransformFilter::FindNormalAxis(int dims[3])