/usr/include/ITK-4.5/itkLabelImageGaussianInterpolateImageFunction.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 | /*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: $RCSfile: itkLabelImageGaussianInterpolateImageFunction.h,v $
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 __itkLabelImageGaussianInterpolateImageFunction_h
#define __itkLabelImageGaussianInterpolateImageFunction_h
#include "itkGaussianInterpolateImageFunction.h"
namespace itk
{
/** \class LabelImageGaussianInterpolateImageFunction
* \brief Interpolation function for multi-label images that implicitly smooths each
* unique value in the image corresponding to each label set element and returns the
* corresponding label set element with the largest wieght.
*
* This filter is an alternative to nearest neighbor interpolation for multi-label
* images. Given a multi-label image \c I with label set \c L, this function returns a
* label at the non-voxel position \c I(x), based on the following rule
*
* \f[
* I(x) = \arg\max_{l \in L} (G_\sigma * I_l)(x)
* \f]
*
* Where \f$ I_l \f$ is the \c l-th binary component of the multilabel image. In other words,
* each label in the multi-label image is convolved with a Gaussian, and the label
* for which the response is largest is returned. For sigma=0, this is just nearest
* neighbor interpolation.
*
* This class defines an N-dimensional Gaussian interpolation function for label
* using the vnl error function. The two parameters associated with this function
* are:
* \li \c Sigma - a scalar array of size ImageDimension determining the width
* of the interpolation function.
* \li \c Alpha - a scalar specifying the cutoff distance over which the function
* is calculated.
*
* \note The input image can be of any type, but the number of unique intensity values
* in the image will determine the amount of memory needed to complete each interpolation.
*
*
* \author Paul Yushkevich
* \author Nick Tustison
*
* \ingroup ITKImageFunction
*/
template <typename TInputImage, typename TCoordRep = double,
typename TPixelCompare = std::less<typename itk::NumericTraits<typename TInputImage::PixelType>::RealType> >
class LabelImageGaussianInterpolateImageFunction :
public GaussianInterpolateImageFunction<TInputImage, TCoordRep>
{
public:
/** Standard class typedefs. */
typedef LabelImageGaussianInterpolateImageFunction Self;
typedef GaussianInterpolateImageFunction<TInputImage, TCoordRep> Superclass;
typedef SmartPointer<Self> Pointer;
typedef SmartPointer<const Self> ConstPointer;
typedef typename TInputImage::PixelType InputPixelType;
/** Run-time type information (and related methods). */
itkTypeMacro( LabelImageGaussianInterpolateImageFunction, GaussianInterpolateImageFunction );
/** Method for creation through the object factory. */
itkNewMacro( Self );
/** ImageDimension constant */
itkStaticConstMacro( ImageDimension, unsigned int, TInputImage::ImageDimension );
/** OutputType typedef support. */
typedef typename Superclass::OutputType OutputType;
/** InputImageType typedef support. */
typedef typename Superclass::InputImageType InputImageType;
/** RealType typedef support. */
typedef typename Superclass::RealType RealType;
/** Index typedef support. */
typedef typename Superclass::IndexType IndexType;
/** ContinuousIndex typedef support. */
typedef typename Superclass::ContinuousIndexType ContinuousIndexType;
/** Array typedef support */
typedef typename Superclass::ArrayType ArrayType;
/**
* Evaluate at the given index
*/
virtual OutputType EvaluateAtContinuousIndex(
const ContinuousIndexType & cindex ) const
{
return this->EvaluateAtContinuousIndex( cindex, NULL );
}
protected:
LabelImageGaussianInterpolateImageFunction();
~LabelImageGaussianInterpolateImageFunction(){};
private:
LabelImageGaussianInterpolateImageFunction( const Self& ); //purposely not implemented
void operator=( const Self& ); //purposely not implemented
/**
* Evaluate function value at the given index
*/
virtual OutputType EvaluateAtContinuousIndex(
const ContinuousIndexType &, OutputType * ) const;
};
} // end namespace itk
#ifndef ITK_MANUAL_INSTANTIATION
#include "itkLabelImageGaussianInterpolateImageFunction.hxx"
#endif
#endif
|