/usr/include/ITK-4.5/itkDemonsImageToImageMetricv4.h is in libinsighttoolkit4-dev 4.5.0-3.
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 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 | /*=========================================================================
*
* 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 __itkDemonsImageToImageMetricv4_h
#define __itkDemonsImageToImageMetricv4_h
#include "itkImageToImageMetricv4.h"
#include "itkDemonsImageToImageMetricv4GetValueAndDerivativeThreader.h"
namespace itk
{
/** \class DemonsImageToImageMetricv4
*
* \brief Class implementing demons metric.
*
* The implementation is taken from itkDemonsRegistrationFunction.
*
* The metric derivative can be calculated using image derivatives
* either from the fixed or moving images. The default is to use fixed-image
* gradients. See ObjectToObjectMetric::SetGradientSource to change
* this behavior.
*
* An intensity threshold is used, below which image pixels are considered
* equal for the purpose of derivative calculation. The threshold can be
* changed by calling SetIntensityDifferenceThreshold.
*
* \note This metric supports only moving transforms with local support and
* with a number of local parameters that matches the moving image dimension.
* In particular, it's meant to be used with itkDisplacementFieldTransform and
* derived classes.
*
* See DemonsImageToImageMetricv4GetValueAndDerivativeThreader::ProcessPoint
* for algorithm implementation.
*
* \sa itkImageToImageMetricv4
* \ingroup ITKMetricsv4
*/
template <typename TFixedImage, typename TMovingImage, typename TVirtualImage = TFixedImage,
typename TInternalComputationValueType = double,
typename TMetricTraits = DefaultImageToImageMetricTraitsv4<TFixedImage,TMovingImage,TVirtualImage,TInternalComputationValueType>
>
class DemonsImageToImageMetricv4 :
public ImageToImageMetricv4<TFixedImage, TMovingImage, TVirtualImage, TInternalComputationValueType, TMetricTraits>
{
public:
/** Standard class typedefs. */
typedef DemonsImageToImageMetricv4 Self;
typedef ImageToImageMetricv4<TFixedImage, TMovingImage, TVirtualImage,
TInternalComputationValueType,TMetricTraits> Superclass;
typedef SmartPointer<Self> Pointer;
typedef SmartPointer<const Self> ConstPointer;
/** Method for creation through the object factory. */
itkNewMacro(Self);
/** Run-time type information (and related methods). */
itkTypeMacro(DemonsImageToImageMetricv4, ImageToImageMetricv4);
/** Superclass types */
typedef typename Superclass::MeasureType MeasureType;
typedef typename Superclass::DerivativeType DerivativeType;
typedef typename Superclass::FixedImagePointType FixedImagePointType;
typedef typename Superclass::FixedImagePixelType FixedImagePixelType;
typedef typename Superclass::FixedImageGradientType FixedImageGradientType;
typedef typename Superclass::MovingImagePointType MovingImagePointType;
typedef typename Superclass::MovingImagePixelType MovingImagePixelType;
typedef typename Superclass::MovingImageGradientType MovingImageGradientType;
typedef typename Superclass::MovingTransformType MovingTransformType;
typedef typename Superclass::JacobianType JacobianType;
typedef typename Superclass::VirtualImageType VirtualImageType;
typedef typename Superclass::VirtualIndexType VirtualIndexType;
typedef typename Superclass::VirtualPointType VirtualPointType;
typedef typename Superclass::VirtualPointSetType VirtualSPointSetType;
typedef typename Superclass::NumberOfParametersType NumberOfParametersType;
/** It should be possible to derive the internal computation type from the class object. */
typedef TInternalComputationValueType InternalComputationValueType;
typedef typename Superclass::ImageDimensionType ImageDimensionType;
/* Image dimension accessors */
itkStaticConstMacro(VirtualImageDimension, ImageDimensionType,
TVirtualImage::ImageDimension);
itkStaticConstMacro(FixedImageDimension, ImageDimensionType,
TFixedImage::ImageDimension);
itkStaticConstMacro(MovingImageDimension, ImageDimensionType,
TMovingImage::ImageDimension);
virtual void Initialize(void) throw ( itk::ExceptionObject );
/** Accessors for the image intensity difference threshold use
* in derivative calculation */
itkGetConstMacro(IntensityDifferenceThreshold, TInternalComputationValueType);
itkSetMacro(IntensityDifferenceThreshold, TInternalComputationValueType);
/** Get the denominator threshold used in derivative calculation. */
itkGetConstMacro(DenominatorThreshold, TInternalComputationValueType);
protected:
itkGetConstMacro(Normalizer, TInternalComputationValueType);
DemonsImageToImageMetricv4();
virtual ~DemonsImageToImageMetricv4();
friend class DemonsImageToImageMetricv4GetValueAndDerivativeThreader< ThreadedImageRegionPartitioner< Superclass::VirtualImageDimension >, Superclass, Self >;
friend class DemonsImageToImageMetricv4GetValueAndDerivativeThreader< ThreadedIndexedContainerPartitioner, Superclass, Self >;
typedef DemonsImageToImageMetricv4GetValueAndDerivativeThreader< ThreadedImageRegionPartitioner< Superclass::VirtualImageDimension >, Superclass, Self >
DemonsDenseGetValueAndDerivativeThreaderType;
typedef DemonsImageToImageMetricv4GetValueAndDerivativeThreader< ThreadedIndexedContainerPartitioner, Superclass, Self >
DemonsSparseGetValueAndDerivativeThreaderType;
void PrintSelf(std::ostream& os, Indent indent) const;
private:
/** Threshold below which the denominator term is considered zero.
* Fixed programatically in constructor. */
TInternalComputationValueType m_DenominatorThreshold;
/** Threshold below which two intensity value are assumed to match. */
TInternalComputationValueType m_IntensityDifferenceThreshold;
/* Used to normalize derivative calculation. Automatically calculated */
TInternalComputationValueType m_Normalizer;
DemonsImageToImageMetricv4(const Self &); //purposely not implemented
void operator = (const Self &); //purposely not implemented
};
} // end namespace itk
#ifndef ITK_MANUAL_INSTANTIATION
#include "itkDemonsImageToImageMetricv4.hxx"
#endif
#endif
|