itkScalarImageToTextureFeaturesFilter.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 itkScalarImageToTextureFeaturesFilter_h
#define itkScalarImageToTextureFeaturesFilter_h
 
#include "itkDataObjectDecorator.h"
 
#include "itkHistogramToTextureFeaturesFilter.h"
#include "itkScalarImageToCooccurrenceMatrixFilter.h"
 
namespace itk::Statistics
{
 
template <typename TImageType,
          typename THistogramFrequencyContainer = DenseFrequencyContainer2,
          typename TMaskImageType = TImageType>
class ITK_TEMPLATE_EXPORT ScalarImageToTextureFeaturesFilter : public ProcessObject
{
public:
  using Self = ScalarImageToTextureFeaturesFilter;
  using Superclass = ProcessObject;
  using Pointer = SmartPointer<Self>;
  using ConstPointer = SmartPointer<const Self>;

  itkOverrideGetNameOfClassMacro(ScalarImageToTextureFeaturesFilter);

  itkNewMacro(Self);
 
  using FrequencyContainerType = THistogramFrequencyContainer;
  using ImageType = TImageType;
  using ImagePointer = typename ImageType::Pointer;
  using MaskImageType = TMaskImageType;
  using MaskPointer = typename MaskImageType::Pointer;
  using MaskPixelType = typename MaskImageType::PixelType;
 
  using PixelType = typename ImageType::PixelType;
  using OffsetType = typename ImageType::OffsetType;
  using OffsetVector = VectorContainer<unsigned char, OffsetType>;
  using OffsetVectorPointer = typename OffsetVector::Pointer;
  using OffsetVectorConstPointer = typename OffsetVector::ConstPointer;
 
  using CooccurrenceMatrixFilterType =
    ScalarImageToCooccurrenceMatrixFilter<ImageType, FrequencyContainerType, MaskImageType>;
 
  using HistogramType = typename CooccurrenceMatrixFilterType::HistogramType;
  using TextureFeaturesFilterType = HistogramToTextureFeaturesFilter<HistogramType>;
 
  // More work needs to be done to fix wrapping
  // using TextureFeatureName = itk::Statistics::TextureFeatureEnum;
  using TextureFeatureName = uint8_t;
  using FeatureNameVector = VectorContainer<unsigned char, TextureFeatureName>;
 
  using FeatureNameVectorPointer = typename FeatureNameVector::Pointer;
  using FeatureNameVectorConstPointer = typename FeatureNameVector::ConstPointer;
  using FeatureValueVector = VectorContainer<unsigned char, double>;
  using FeatureValueVectorPointer = typename FeatureValueVector::Pointer;

  using DataObjectPointer = DataObject::Pointer;

  using FeatureValueVectorDataObjectType = DataObjectDecorator<FeatureValueVector>;
 
  const FeatureValueVectorDataObjectType *
  GetFeatureMeansOutput() const;
 
  const FeatureValueVectorDataObjectType *
  GetFeatureStandardDeviationsOutput() const;

  using Superclass::SetInput;
  void
  SetInput(const ImageType *);
 
  const ImageType *
  GetInput() const;

  itkGetConstReferenceObjectMacro(FeatureMeans, FeatureValueVector);
  itkGetConstReferenceObjectMacro(FeatureStandardDeviations, FeatureValueVector);

  itkSetConstObjectMacro(RequestedFeatures, FeatureNameVector);
  itkGetConstObjectMacro(RequestedFeatures, FeatureNameVector);

  itkSetConstObjectMacro(Offsets, OffsetVector);
  itkGetConstObjectMacro(Offsets, OffsetVector);

  void
  SetNumberOfBinsPerAxis(unsigned int);

  void
  SetPixelValueMinMax(PixelType min, PixelType max);

  void
  SetMaskImage(const MaskImageType *);
 
  const MaskImageType *
  GetMaskImage() const;

  void
  SetInsidePixelValue(MaskPixelType insidePixelValue);
 
  itkGetConstMacro(FastCalculations, bool);
  itkSetMacro(FastCalculations, bool);
  itkBooleanMacro(FastCalculations);
 
protected:
  ScalarImageToTextureFeaturesFilter();
  ~ScalarImageToTextureFeaturesFilter() override = default;
  void
  PrintSelf(std::ostream & os, Indent indent) const override;
 
  void
  FastCompute();
 
  void
  FullCompute();

  void
  GenerateData() override;

  using DataObjectPointerArraySizeType = ProcessObject::DataObjectPointerArraySizeType;
  using Superclass::MakeOutput;
  DataObjectPointer MakeOutput(DataObjectPointerArraySizeType) override;
 
private:
  typename CooccurrenceMatrixFilterType::Pointer m_GLCMGenerator{};
 
  typename TextureFeaturesFilterType::Pointer m_GLCMCalculator{};
 
  FeatureValueVectorPointer     m_FeatureMeans{};
  FeatureValueVectorPointer     m_FeatureStandardDeviations{};
  FeatureNameVectorConstPointer m_RequestedFeatures{};
  OffsetVectorConstPointer      m_Offsets{};
  bool                          m_FastCalculations{};
};
} // namespace itk::Statistics
 
#ifndef ITK_MANUAL_INSTANTIATION
#  include "itkScalarImageToTextureFeaturesFilter.hxx"
#endif
 
#endif