WikiExamples/PointSet/BoundingBox.cxx

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

int main(int, char*[])
{
  typedef float PixelType;
  const unsigned int Dimension = 3;
  typedef itk::PointSet< PixelType, Dimension >   PointSetType;
  typedef PointSetType::PointType PointType;
  typedef PointSetType::PointsContainerPointer 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);
  
  typedef itk::BoundingBox<unsigned long, Dimension, PixelType> BoundingBoxType;
  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;
}