SphinxExamples/src/Core/Common/ImageRegionIntersection/Code.cxx

/*=========================================================================
 *
 *  Copyright NumFOCUS
 *
 *  Licensed under the Apache License, Version 2.0 (the "License");
 *  you may not use this file except in compliance with the License.
 *  You may obtain a copy of the License at
 *
 *         https://www.apache.org/licenses/LICENSE-2.0.txt
 *
 *  Unless required by applicable law or agreed to in writing, software
 *  distributed under the License is distributed on an "AS IS" BASIS,
 *  WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 *  See the License for the specific language governing permissions and
 *  limitations under the License.
 *
 *=========================================================================*/
 
#include "itkImage.h"
 
int
main()
{
  constexpr unsigned int Dimension = 2;
  using PixelType = unsigned char;
 
  using ImageType = itk::Image<PixelType, Dimension>;
 
  // Big region
  ImageType::RegionType bigRegion;
 
  ImageType::SizeType bigSize;
  bigSize[0] = 100;
  bigSize[1] = 100;
 
  ImageType::IndexType bigStart;
  bigStart[0] = 0;
  bigStart[1] = 0;
 
  bigRegion.SetSize(bigSize);
  bigRegion.SetIndex(bigStart);
 
  // Small inside region
  ImageType::RegionType smallInsideRegion;
 
  ImageType::SizeType smallInsideSize;
  smallInsideSize[0] = 10;
  smallInsideSize[1] = 10;
 
  ImageType::IndexType smallInsideStart;
  smallInsideStart[0] = 50;
  smallInsideStart[1] = 50;
 
  smallInsideRegion.SetSize(smallInsideSize);
  smallInsideRegion.SetIndex(smallInsideStart);
 
  std::cout << "Small inside region is " << bigRegion.IsInside(smallInsideRegion) << std::endl;
 
  // Small outside region
  ImageType::RegionType smallOutsideRegion;
 
  ImageType::SizeType smallOutsideSize;
  smallOutsideSize[0] = 10;
  smallOutsideSize[1] = 10;
 
  ImageType::IndexType smallOutsideStart;
  smallOutsideStart[0] = 110;
  smallOutsideStart[1] = 110;
 
  smallOutsideRegion.SetSize(smallOutsideSize);
  smallOutsideRegion.SetIndex(smallOutsideStart);
 
  std::cout << "Small outside region is " << bigRegion.IsInside(smallOutsideRegion) << std::endl;
 
  // Small overlap region
  ImageType::RegionType smallOverlapRegion;
 
  ImageType::SizeType smallOverlapSize;
  smallOverlapSize[0] = 10;
  smallOverlapSize[1] = 10;
 
  ImageType::IndexType smallOverlapStart;
  smallOverlapStart[0] = 97;
  smallOverlapStart[1] = 97;
 
  smallOverlapRegion.SetSize(smallOverlapSize);
  smallOverlapRegion.SetIndex(smallOverlapStart);
 
  std::cout << "Small overlap region is " << bigRegion.IsInside(smallOverlapRegion) << std::endl;
 
  return EXIT_SUCCESS;
}