/*=========================================================================
 *
 *  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.
 *
 *=========================================================================*/

/*
This file test the various inclusion strategies available
to itkFloodFilledSpatialFunctionConditionalIterator.
*/


#include "itkImageRegionIterator.h"
#include "itkSphereSpatialFunction.h"
#include "itkFloodFilledSpatialFunctionConditionalIterator.h"

int
itkFloodFilledSpatialFunctionTest(int, char *[])
{
  constexpr unsigned int dim = 2;

  // Image type alias
  using ImageType = itk::Image<bool, dim>;

  using SizeValueType = ImageType::SizeValueType;
  using SpacingValueType = ImageType::SpacingValueType;
  using PointValueType = ImageType::PointValueType;

  // Image size and spacing parameters
  SizeValueType    sourceImageSize[] = { 5, 5 };
  SpacingValueType sourceImageSpacing[] = { 1.0, 1.0 };
  PointValueType   sourceImageOrigin[] = { 0, 0 };

  // Creates the sourceImage (but doesn't set the size or allocate memory)
  auto sourceImage = ImageType::New();
  sourceImage->SetOrigin(sourceImageOrigin);
  sourceImage->SetSpacing(sourceImageSpacing);

  // Create a size object native to the sourceImage type
  ImageType::SizeType sourceImageSizeObject;

  // Set the size object to the array defined earlier
  sourceImageSizeObject.SetSize(sourceImageSize);

  // Create a region object native to the sourceImage type
  ImageType::RegionType largestPossibleRegion;

  // Resize the region
  largestPossibleRegion.SetSize(sourceImageSizeObject);

  // Set the largest legal region size (i.e. the size of the whole sourceImage), the buffered, and
  // the requested region to what we just defined.
  sourceImage->SetRegions(largestPossibleRegion);

  // Now allocate memory for the sourceImage
  sourceImage->Allocate();

  // Loop over all available iterator strategies
  for (int strat = 0; strat < 4; ++strat)
  {
    // Initialize the image to hold all 0's
    itk::ImageRegionIterator<ImageType> it(sourceImage, largestPossibleRegion);

    for (it.GoToBegin(); !it.IsAtEnd(); ++it)
    {
      it.Set(false);
    }

    // Create and initialize a spatial function
    using FunctionType = itk::SphereSpatialFunction<dim>;
    using FunctionPositionType = FunctionType::InputType;

    auto spatialFunc = FunctionType::New();
    spatialFunc->SetRadius(1.0);

    FunctionPositionType center;
    center[0] = 2.5;
    center[1] = 2.5;
    spatialFunc->SetCenter(center);

    // Create and initialize a spatial function iterator
    ImageType::IndexType            seedPos;
    const ImageType::IndexValueType pos[] = { 2, 2 };
    seedPos.SetIndex(pos);

    using ItType = itk::FloodFilledSpatialFunctionConditionalIterator<ImageType, FunctionType>;

    ItType sfi(sourceImage, spatialFunc, seedPos);

    switch (strat)
    {
      case 0:
      {
        sfi.SetOriginInclusionStrategy();
      }
      break;
      case 1:
      {
        sfi.SetCenterInclusionStrategy();
      }
      break;
      case 2:
      {
        sfi.SetCompleteInclusionStrategy();
      }
      break;
      case 3:
      {
        sfi.SetIntersectInclusionStrategy();
      }
    } // end switch inclusion strategy

    // Iterate through the entire image and set interior pixels to 1
    for (sfi.GoToBegin(); !(sfi.IsAtEnd()); ++sfi)
    {
      sfi.Set(true);
    }

  } // end loop over iterator strategies

  return EXIT_SUCCESS;
}
