WikiExamples/PointSet/BoundingBox.cxx

#include "itkPoint.h"
#include "itkPointSet.h"
#include "itkBoundingBox.h"

int main(int, char*[])
{
  using PixelType = float;
  constexpr unsigned int Dimension = 3;
  using PointSetType = itk::PointSet< PixelType, Dimension >;
  using PointType = PointSetType::PointType;
  using PointsContainerPointer = PointSetType::PointsContainerPointer;
  
  PointSetType::Pointer  pointSet = PointSetType::New();
  PointsContainerPointer  points = pointSet->GetPoints();
  
  // Create points
  PointType p0, p1, p2;

  p0[0]=  0.0; p0[1]= 0.0; p0[2]= 0.0;
  p1[0]=  0.1; p1[1]= 0.0; p1[2]= 0.0;
  p2[0]=  0.0; p2[1]= 0.1; p2[2]= 0.0;
  
  points->InsertElement(0, p0);
  points->InsertElement(1, p1);
  points->InsertElement(2, p2);
  
  using BoundingBoxType = itk::BoundingBox<itk::IdentifierType, Dimension, PixelType>;
  BoundingBoxType::Pointer boundingBox = BoundingBoxType::New();
  boundingBox->SetPoints(points);
  boundingBox->ComputeBoundingBox();
  std::cout << "bounds: " << boundingBox->GetBounds() << std::endl;
  std::cout << "center: " << boundingBox->GetCenter() << std::endl;
  std::cout << "diagonal length squared: " << boundingBox->GetDiagonalLength2() << std::endl;
  
  return EXIT_SUCCESS;
}