itkTriangleMeshToBinaryImageFilter.h

Go to the documentation of this file.

/*=========================================================================
 *
 *  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 itkTriangleMeshToBinaryImageFilter_h
#define itkTriangleMeshToBinaryImageFilter_h
 
#include "itkImageSource.h"
 
#include "itkPolygonCell.h"
#include "itkMapContainer.h"
#include "itkVectorContainer.h"
#include "itkAutomaticTopologyMeshSource.h"
#include "itkPointSet.h"
 
#include <vector>
 
namespace itk
{
class Point1D
{
public:
  double m_X;
  int    m_Sign;
 
  Point1D() = default;
  Point1D(const double p, const int s)
    : m_X(p)
    , m_Sign(s)
  {}
 
  Point1D(const Point1D & point) = default;
 
  [[nodiscard]] double
  getX() const
  {
    return m_X;
  }
 
  [[nodiscard]] int
  getSign() const
  {
    return m_Sign;
  }
};

template <typename TInputMesh, typename TOutputImage>
class ITK_TEMPLATE_EXPORT TriangleMeshToBinaryImageFilter : public ImageSource<TOutputImage>
{
public:
  ITK_DISALLOW_COPY_AND_MOVE(TriangleMeshToBinaryImageFilter);

  using Self = TriangleMeshToBinaryImageFilter;
  using Superclass = ImageSource<TOutputImage>;
  using Pointer = SmartPointer<Self>;
  using ConstPointer = SmartPointer<const Self>;
 
  using IndexType = typename TOutputImage::IndexType;
  using SizeType = typename TOutputImage::SizeType;
  using OutputImageType = TOutputImage;
  using OutputImagePointer = typename OutputImageType::Pointer;
  using ValueType = typename OutputImageType::ValueType;
  using SpacingType = typename OutputImageType::SpacingType;
  using DirectionType = typename OutputImageType::DirectionType;

  itkNewMacro(Self);

  itkOverrideGetNameOfClassMacro(TriangleMeshToBinaryImageFilter);

  using typename Superclass::OutputImageRegionType;

  using InputMeshType = TInputMesh;
  using InputMeshPointer = typename InputMeshType::Pointer;
  using InputPointType = typename InputMeshType::PointType;
  using InputPixelType = typename InputMeshType::PixelType;
  using InputCellTraitsType = typename InputMeshType::MeshTraits::CellTraits;
  using CellType = typename InputMeshType::CellType;
  using CellsContainerPointer = typename InputMeshType::CellsContainerPointer;
  using CellsContainerIterator = typename InputMeshType::CellsContainerIterator;
 
  using InputPointsContainer = typename InputMeshType::PointsContainer;
  using InputPointsContainerPointer = typename InputPointsContainer::Pointer;
  using InputPointsContainerIterator = typename InputPointsContainer::Iterator;
 
  using PointSetType = itk::PointSet<double, 3>;
  using PointsContainer = typename PointSetType::PointsContainer;
 
  using PointType = itk::Point<double, 3>;
  using Point2DType = itk::Point<double, 2>;
 
  using DoubleArrayType = itk::Array<double>;
 
  using Point1DVector = std::vector<Point1D>;
  using Point1DArray = std::vector<std::vector<Point1D>>;
 
  using Point2DVector = std::vector<Point2DType>;
  using Point2DArray = std::vector<std::vector<Point2DType>>;
 
  using PointVector = std::vector<PointType>;
  using PointArray = std::vector<std::vector<PointType>>;

  itkSetMacro(Spacing, SpacingType);
  virtual void
  SetSpacing(const double spacing[3]);
  virtual void
  SetSpacing(const float spacing[3]);
 
  itkGetConstReferenceMacro(Spacing, SpacingType);

  itkSetMacro(Direction, DirectionType);
  itkGetConstMacro(Direction, DirectionType);
  itkSetMacro(InsideValue, ValueType);
  itkGetConstMacro(InsideValue, ValueType);
  itkSetMacro(OutsideValue, ValueType);
  itkGetConstMacro(OutsideValue, ValueType);
  itkSetMacro(Origin, PointType);
  virtual void
  SetOrigin(const double origin[3]);
  virtual void
  SetOrigin(const float origin[3]);
 
  itkGetConstReferenceMacro(Origin, PointType);

  itkSetMacro(Index, IndexType);
  itkGetConstMacro(Index, IndexType);
  itkSetMacro(Size, SizeType);
  itkGetConstMacro(Size, SizeType);
  using Superclass::SetInput;
  void
  SetInput(InputMeshType * input);
 
  void
  SetInfoImage(OutputImageType * InfoImage)
  {
    if (InfoImage != m_InfoImage)
    {
      this->Modified();
      m_InfoImage = InfoImage;
    }
  }

  InputMeshType *
  GetInput();
 
  InputMeshType *
  GetInput(unsigned int idx);
 
  /* Set the tolerance for doing spatial searches of the polydata. */
  itkSetMacro(Tolerance, double);
  itkGetConstMacro(Tolerance, double);
 
protected:
  TriangleMeshToBinaryImageFilter();
  ~TriangleMeshToBinaryImageFilter() override = default;
 
  void
  GenerateOutputInformation() override
  {} // do nothing
  void
  GenerateData() override;
 
  virtual void
  RasterizeTriangles();

  static int
  PolygonToImageRaster(PointVector coords, Point1DArray & zymatrix, int extent[6]);
 
  OutputImageType * m_InfoImage{};
 
  IndexType m_Index{};
 
  SizeType m_Size{};
 
  SpacingType m_Spacing{};
 
  PointType m_Origin{}; // start value
 
  double m_Tolerance{};
 
  ValueType m_InsideValue{};
  ValueType m_OutsideValue{};
 
  DirectionType m_Direction{};
 
  void
  PrintSelf(std::ostream & os, Indent indent) const override;
 
private:
  static bool
  ComparePoints2D(Point2DType a, Point2DType b);
 
  static bool
  ComparePoints1D(Point1D a, Point1D b);
};
} // end namespace itk
 
#ifndef ITK_MANUAL_INSTANTIATION
#  include "itkTriangleMeshToBinaryImageFilter.hxx"
#endif
 
#endif