itkStrainImageFilter.h

Go to the documentation of this file.

/*=========================================================================
 *
 *  Copyright Insight Software Consortium
 *
 *  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
 *
 *         http://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 itkStrainImageFilter_h
#define itkStrainImageFilter_h

#include "itkCovariantVector.h"
#include "itkImageToImageFilter.h"
#include "itkSymmetricSecondRankTensor.h"
#include "itkSplitComponentsImageFilter.h"

namespace itk
{

template< typename TInputImage, typename TOperatorValueType=float, typename TOutputValueType=float >
class StrainImageFilter : public
  ImageToImageFilter< TInputImage,
    Image< SymmetricSecondRankTensor< TOutputValueType, TInputImage::ImageDimension >,
                                  TInputImage::ImageDimension > >
{
public:
  ITK_DISALLOW_COPY_AND_ASSIGN(StrainImageFilter);

  static constexpr unsigned int ImageDimension = TInputImage::ImageDimension;

  using InputImageType = TInputImage;
  using OutputPixelType = SymmetricSecondRankTensor< TOutputValueType, ImageDimension >;
  using OutputImageType = Image< OutputPixelType, ImageDimension >;
  using OperatorImageType = Image< TOperatorValueType, ImageDimension >;

  using Self = StrainImageFilter;
  using Superclass = ImageToImageFilter< InputImageType, OutputImageType >;

  using Pointer = SmartPointer< Self >;
  using ConstPointer = SmartPointer< const Self >;


  using GradientOutputPixelType = CovariantVector< TOperatorValueType, ImageDimension >;
  using GradientOutputImageType = Image< GradientOutputPixelType, ImageDimension >;
  using GradientFilterType = ImageToImageFilter< OperatorImageType, GradientOutputImageType >;

  using VectorGradientFilterType = ImageToImageFilter< InputImageType, GradientOutputImageType >;

  itkNewMacro( Self );

  itkTypeMacro( StrainImageFilter, ImageToImageFilter );

  itkSetObjectMacro( GradientFilter, GradientFilterType );
  itkGetConstObjectMacro( GradientFilter, GradientFilterType );

  itkSetObjectMacro( VectorGradientFilter, VectorGradientFilterType );
  itkGetConstObjectMacro( VectorGradientFilter, VectorGradientFilterType );

  enum StrainFormType {INFINITESIMAL = 0, GREENLAGRANGIAN = 1, EULERIANALMANSI = 2};

  itkSetMacro( StrainForm, StrainFormType );
  itkGetConstMacro( StrainForm, StrainFormType );

protected:
  using OutputRegionType = typename OutputImageType::RegionType;

  StrainImageFilter();

  void BeforeThreadedGenerateData() override;

  void DynamicThreadedGenerateData( const OutputRegionType& outputRegion ) override;

  using InputComponentsImageFilterType = itk::SplitComponentsImageFilter< InputImageType, OperatorImageType >;

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

private:
  typename InputComponentsImageFilterType::Pointer m_InputComponentsFilter;

  typename GradientFilterType::Pointer m_GradientFilter;

  typename VectorGradientFilterType::Pointer m_VectorGradientFilter;

  StrainFormType m_StrainForm;
};

} // end namespace itk

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

#endif