/usr/include/InsightToolkit/Common/itkElasticBodyReciprocalSplineKernelTransform.h is in libinsighttoolkit3-dev 3.20.1-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 | /*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: itkElasticBodyReciprocalSplineKernelTransform.h
Language: C++
Date: $Date$
Version: $Revision$
Copyright (c) Insight Software Consortium. All rights reserved.
See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.
This software is distributed WITHOUT ANY WARRANTY; without even
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
PURPOSE. See the above copyright notices for more information.
=========================================================================*/
#ifndef __itkElasticBodyReciprocalSplineKernelTransform_h
#define __itkElasticBodyReciprocalSplineKernelTransform_h
#include "itkKernelTransform.h"
namespace itk
{
/** \class ElasticBodyReciprocalSplineKernelTransform
* This class defines the elastic body spline (EBS) transformation.
* It is implemented in as straightforward a manner as possible from
* the IEEE TMI paper by Davis, Khotanzad, Flamig, and Harms,
* Vol. 16 No. 3 June 1997
* Taken from the paper:
* The EBS "is based on a physical model of a homogeneous, isotropic,
* three-dimensional elastic body. The model can approximate the way
* that some physical objects deform".
*
* \ingroup Transforms
*/
template <class TScalarType = double, // Data type for scalars (float or double)
unsigned int NDimensions = 3> // Number of dimensions
class ITK_EXPORT ElasticBodyReciprocalSplineKernelTransform :
public KernelTransform< TScalarType, NDimensions>
{
public:
/** Standard class typedefs. */
typedef ElasticBodyReciprocalSplineKernelTransform Self;
typedef KernelTransform< TScalarType,
NDimensions> Superclass;
typedef SmartPointer<Self> Pointer;
typedef SmartPointer<const Self> ConstPointer;
/** Run-time type information (and related methods). */
itkTypeMacro( ElasticBodyReciprocalSplineKernelTransform, KernelTransform );
/** New macro for creation of through a Smart Pointer */
itkNewMacro( Self );
/** Scalar type. */
typedef typename Superclass::ScalarType ScalarType;
/** Parameters type. */
typedef typename Superclass::ParametersType ParametersType;
/** Jacobian type. */
typedef typename Superclass::JacobianType JacobianType;
/** Dimension of the domain space. */
itkStaticConstMacro(SpaceDimension, unsigned int,Superclass::SpaceDimension);
/** Set alpha. Alpha is related to Poisson's Ratio (\f$\nu\f$) as
* \f$\alpha = 8 ( 1 - \nu ) - 1\f$
*/
itkSetMacro( Alpha, TScalarType );
/** Get alpha */
itkGetConstMacro( Alpha, TScalarType );
/** These (rather redundant) typedefs are needed because on SGI, typedefs
* are not inherited */
typedef typename Superclass::InputPointType InputPointType;
typedef typename Superclass::OutputPointType OutputPointType;
typedef typename Superclass::InputVectorType InputVectorType;
typedef typename Superclass::OutputVectorType OutputVectorType;
typedef typename Superclass::InputCovariantVectorType InputCovariantVectorType;
typedef typename Superclass::OutputCovariantVectorType OutputCovariantVectorType;
protected:
ElasticBodyReciprocalSplineKernelTransform();
virtual ~ElasticBodyReciprocalSplineKernelTransform();
void PrintSelf(std::ostream& os, Indent indent) const;
/** These (rather redundant) typedefs are needed because on SGI, typedefs
* are not inherited */
typedef typename Superclass::GMatrixType GMatrixType;
/** Compute G(x)
* For the elastic body spline, this is:
* G(x) = [alpha*r(x)*I - 3*x*x'/r(x)]
* \f$ G(x) = [\alpha*r(x)*I - 3*x*x'/r(x) ]\f$
* where
* \f$\alpha = 8 ( 1 - \nu ) - 1\f$
* \f$\nu\f$ is Poisson's Ratio
* r(x) = Euclidean norm = sqrt[x1^2 + x2^2 + x3^2]
* \f[ r(x) = \sqrt{ x_1^2 + x_2^2 + x_3^2 } \f]
* I = identity matrix */
virtual void ComputeG(const InputVectorType& landmarkVector, GMatrixType & gmatrix) const;
/**
* \deprecated in ITK 3.6, please use void ComputeG(vector,gmatrix) instead.
*/
itkLegacyMacro( virtual const GMatrixType & ComputeG(const InputVectorType& landmarkVector) const );
/** alpha, Poisson's ratio */
TScalarType m_Alpha;
private:
ElasticBodyReciprocalSplineKernelTransform(const Self&); //purposely not implemented
void operator=(const Self&); //purposely not implemented
};
} // namespace itk
// Define instantiation macro for this template.
#define ITK_TEMPLATE_ElasticBodyReciprocalSplineKernelTransform(_, EXPORT, x, y) namespace itk { \
_(2(class EXPORT ElasticBodyReciprocalSplineKernelTransform< ITK_TEMPLATE_2 x >)) \
namespace Templates { typedef ElasticBodyReciprocalSplineKernelTransform< ITK_TEMPLATE_2 x > \
ElasticBodyReciprocalSplineKernelTransform##y; } \
}
#if ITK_TEMPLATE_EXPLICIT
# include "Templates/itkElasticBodyReciprocalSplineKernelTransform+-.h"
#endif
#if ITK_TEMPLATE_TXX
# include "itkElasticBodyReciprocalSplineKernelTransform.txx"
#endif
#endif // __itkElasticBodyReciprocalSplineKernelTransform_h
|