#include <itkRecursiveBSplineTransformImplementation.h>
This helper class contains the actual implementation of the recursive B-spline transform.
Define the end case for SpaceDimension = 0.
Compared to the RecursiveBSplineTransformImplementation class, this class works as a vector operator, and is therefore also templated over the OutputDimension.
Note: More optimized code can be found in itkRecursiveBSplineImplementation.h
Definition at line 46 of file itkRecursiveBSplineTransformImplementation.h.
Public Types | |
using | InternalFloatType = double |
using | RecursiveBSplineWeightFunctionType |
Public Member Functions | |
itkStaticConstMacro (BSplineNumberOfIndices, unsigned int, RecursiveBSplineWeightFunctionType::NumberOfIndices) | |
itkStaticConstMacro (HelperConstVariable, unsigned int,(SpaceDimension - 1) *(SplineOrder+1)) | |
Static Public Member Functions | |
static void | ComputeNonZeroJacobianIndices (unsigned long *&nzji, const unsigned long parametersPerDim, unsigned long currentIndex, const OffsetValueType *const gridOffsetTable) |
static void | EvaluateJacobianWithImageGradientProduct (TScalar *&imageJacobian, const InternalFloatType *const movingImageGradient, const double *const weights1D, const double value) |
static void | GetJacobian (TScalar *&jacobians, const double *const weights1D, const double value) |
static void | GetJacobianOfSpatialHessian (InternalFloatType *&jsh_out, const double *const weights1D, const double *const derivativeWeights1D, const double *const hessianWeights1D, const double *const directionCosines, const InternalFloatType *const jsh) |
static void | GetJacobianOfSpatialJacobian (InternalFloatType *&jsj_out, const double *const weights1D, const double *const derivativeWeights1D, const double *const directionCosines, const InternalFloatType *const jsj) |
static void | GetSpatialHessian (InternalFloatType *const sh, const TScalar *const *const mu, const OffsetValueType *const gridOffsetTable, const double *const weights1D, const double *const derivativeWeights1D, const double *const hessianWeights1D) |
static void | GetSpatialJacobian (InternalFloatType *const sj, const TScalar *const *const mu, const OffsetValueType *const gridOffsetTable, const double *const weights1D, const double *const derivativeWeights1D) |
static void | TransformPoint (TScalar *const opp, const TScalar *const *const mu, const OffsetValueType *const gridOffsetTable, const double *const weights1D) |
using itk::RecursiveBSplineTransformImplementation< OutputDimension, SpaceDimension, SplineOrder, TScalar >::InternalFloatType = double |
Definition at line 49 of file itkRecursiveBSplineTransformImplementation.h.
using itk::RecursiveBSplineTransformImplementation< OutputDimension, SpaceDimension, SplineOrder, TScalar >::RecursiveBSplineWeightFunctionType |
Typedef to know the number of indices at compile time.
Definition at line 55 of file itkRecursiveBSplineTransformImplementation.h.
|
inlinestatic |
ComputeNonZeroJacobianIndices recursive implementation.
Definition at line 125 of file itkRecursiveBSplineTransformImplementation.h.
|
inlinestatic |
EvaluateJacobianWithImageGradientProduct recursive implementation.
Definition at line 108 of file itkRecursiveBSplineTransformImplementation.h.
|
inlinestatic |
GetJacobian recursive implementation.
Definition at line 95 of file itkRecursiveBSplineTransformImplementation.h.
|
inlinestatic |
GetJacobianOfSpatialHessian recursive implementation. Multiplication with the direction cosines is performed in the end - case.
Definition at line 305 of file itkRecursiveBSplineTransformImplementation.h.
|
inlinestatic |
GetJacobianOfSpatialJacobian recursive implementation. Multiplication with the direction cosines is performed in the end-case.
Definition at line 269 of file itkRecursiveBSplineTransformImplementation.h.
|
inlinestatic |
GetSpatialHessian recursive implementation. As an (almost) free by-product this function delivers the displacement, i.e. the TransformPoint() function, as well as the SpatialJacobian.
Specifically, sh is the output argument. It should be allocated with a size OutputDimension * ( SpaceDimension + 1 ) * ( SpaceDimension + 2 ) / 2. sh should point to allocated memory, but this function initializes sh.
Upon return sh contains the spatial Hessian, spatial Jacobian and transformpoint. With Hk = [ transformPoint spatialJacobian' spatialJacobian spatialHessian ] . (Hk specifies all info of dimension (element) k (< OutputDimension) of the point and spatialJacobian is a vector of the derivative of this point with respect to the dimensions.) The i,j (both < SpaceDimension) element of Hk is stored in: i<=j : sh[ k + OutputDimension * (i + j*(j+1)/2 ) ] i>=j : sh[ k + OutputDimension * (j + i*(i+1)/2 ) ]
Note that we store only one of the symmetric halves of Hk.
Definition at line 209 of file itkRecursiveBSplineTransformImplementation.h.
|
inlinestatic |
GetSpatialJacobian recursive implementation. As an (almost) free by-product this function delivers the displacement, i.e. the TransformPoint() function.
Definition at line 147 of file itkRecursiveBSplineTransformImplementation.h.
itk::RecursiveBSplineTransformImplementation< OutputDimension, SpaceDimension, SplineOrder, TScalar >::itkStaticConstMacro | ( | BSplineNumberOfIndices | , |
unsigned int | , | ||
RecursiveBSplineWeightFunctionType::NumberOfIndices | ) |
itk::RecursiveBSplineTransformImplementation< OutputDimension, SpaceDimension, SplineOrder, TScalar >::itkStaticConstMacro | ( | HelperConstVariable | , |
unsigned int | , | ||
(SpaceDimension - 1) *(SplineOrder+1) | ) |
Helper constant variable.
|
inlinestatic |
TransformPoint recursive implementation.
Definition at line 61 of file itkRecursiveBSplineTransformImplementation.h.
Generated on 2024-07-17
for elastix by ![]() |
![]() |