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:
  itkStaticConstMacro(ImageDimension, unsigned int, TInputImage::ImageDimension);

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

  typedef StrainImageFilter                                     Self;
  typedef ImageToImageFilter< InputImageType, OutputImageType > Superclass;

  typedef SmartPointer< Self >                                  Pointer;
  typedef SmartPointer< const Self >                            ConstPointer;


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

  typedef ImageToImageFilter< InputImageType, GradientOutputImageType > VectorGradientFilterType;

  itkNewMacro( Self );

  itkTypeMacro( StrainImageFilter, ImageToImageFilter );

  itkSetObjectMacro( GradientFilter, GradientFilterType );
  itkGetObjectMacro( GradientFilter, GradientFilterType );

  itkSetObjectMacro( VectorGradientFilter, VectorGradientFilterType );
  itkGetObjectMacro( VectorGradientFilter, VectorGradientFilterType );

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

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

protected:
  typedef typename OutputImageType::RegionType OutputRegionType;

  StrainImageFilter();

  virtual void BeforeThreadedGenerateData() ITK_OVERRIDE;

  virtual void ThreadedGenerateData( const OutputRegionType& outputRegion, ThreadIdType threadId ) ITK_OVERRIDE;

  typedef itk::SplitComponentsImageFilter< InputImageType, OperatorImageType >
    InputComponentsImageFilterType;

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

private:
  ITK_DISALLOW_COPY_AND_ASSIGN(StrainImageFilter);

  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