/*=========================================================================
 *
 *  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.
 *
 *=========================================================================*/
#ifndef itkBSplineTransformInitializer_h
#define itkBSplineTransformInitializer_h

#include "itkObject.h"
#include "itkObjectFactory.h"

namespace itk
{
/** \class BSplineTransformInitializer
 * \brief BSplineTransformInitializer is a helper class intended to
 * initialize the control point grid such that it has a physically consistent
 * definition.
 * It sets the transform domain origin, physical dimensions and direction from
 * information obtained from the image. It also sets the mesh size if asked
 * to do so by calling SetTransformDomainMeshSize() before calling
 * InitializeTransform().
 *
 * \author Luis Ibanez
 * \author Nick Tustison
 *
 * \ingroup ITKTransform
 */
template <typename TTransform, typename TImage>
class ITK_TEMPLATE_EXPORT BSplineTransformInitializer : public Object
{
public:
  ITK_DISALLOW_COPY_AND_MOVE(BSplineTransformInitializer);

  /** Standard class type aliases. */
  using Self = BSplineTransformInitializer;
  using Superclass = Object;
  using Pointer = SmartPointer<Self>;
  using ConstPointer = SmartPointer<const Self>;

  /** New macro for creation of through a Smart Pointer. */
  itkNewMacro(Self);

  /** \see LightObject::GetNameOfClass() */
  itkOverrideGetNameOfClassMacro(BSplineTransformInitializer);

  /** Type of the transform to initialize. */
  using TransformType = TTransform;

  /** Types defined from the input image traits. */
  using ImageType = TImage;
  using ImagePointer = typename ImageType::ConstPointer;
  using IndexType = typename ImageType::IndexType;
  using ImagePointType = typename ImageType::PointType;
  using ImagePointCoordRepType = typename ImagePointType::CoordRepType;

  /** Types defined from transform traits. */
  using TransformPointer = typename TransformType::Pointer;
  using RegionType = typename TransformType::RegionType;
  using SizeType = typename RegionType::SizeType;
  using SpacingType = typename TransformType::SpacingType;
  using OriginType = typename TransformType::OriginType;
  using DirectionType = typename TransformType::DirectionType;
  using PhysicalDimensionsType = typename TransformType::PhysicalDimensionsType;
  using MeshSizeType = typename TransformType::MeshSizeType;
  using SpacingComponentType = typename SpacingType::ComponentType;

  /** Dimension of parameters. */
  static constexpr unsigned int SpaceDimension = TransformType::SpaceDimension;

  /** Set/Get the transform to be initialized. */
  itkGetConstObjectMacro(Transform, TransformType);
  itkSetObjectMacro(Transform, TransformType);

  /** Set/Get the image to initialize the domain. */
  itkGetConstObjectMacro(Image, ImageType);
  itkSetConstObjectMacro(Image, ImageType);

  /** Allow the user to set the mesh size of the transform via the initializer,
   * even though the initializer does not do anything with that information.
   * Default size is 1^ImageDimension. */
  itkGetConstMacro(TransformDomainMeshSize, MeshSizeType);
  void
  SetTransformDomainMeshSize(const MeshSizeType);

  /** Initialize the transform using the specified transformation domain. */
  virtual void
  InitializeTransform() const;

protected:
  BSplineTransformInitializer() = default;
  ~BSplineTransformInitializer() override = default;

  void
  PrintSelf(std::ostream & os, Indent indent) const override;

private:
  ImagePointer     m_Image{};
  TransformPointer m_Transform{};

  MeshSizeType m_TransformDomainMeshSize{ MeshSizeType::Filled(1) };
  bool         m_SetTransformDomainMeshSizeViaInitializer{ false };

}; // class BSplineTransformInitializer
} // namespace itk

#ifndef ITK_MANUAL_INSTANTIATION
#  include "itkBSplineTransformInitializer.hxx"
#endif

#endif /* itkBSplineTransformInitializer_h */