ITK
6.0.0
Insight Toolkit
|
#include <itkWindowedSincInterpolateImageFunction.h>
Use the windowed sinc function to interpolate.
This function is intended to provide an interpolation function that has minimum aliasing artifacts, in contrast to linear interpolation. According to sampling theory, the infinite-support sinc filter, whose Fourier transform is the box filter, is optimal for resampling a function. In practice, the infinite support sinc filter is approximated using a limited support 'windowed' sinc filter.
\[ I(x,y) = \sum_{i = \lfloor x \rfloor + 1 - m}^{\lfloor x \rfloor + m} \sum_{j = \lfloor y \rfloor + 1 - m}^{\lfloor y \rfloor + m} I_{i,j} K(x-i) K(y-j), \]
\[ K(t) = w(t) \textrm{sinc}(t) = w(t) \frac{\sin(\pi t)}{\pi t} \]
Use this filter the way you would use any ImageInterpolationFunction, so for instance, you can plug it into the ResampleImageFilter class. In order to initialize the filter you must choose several template parameters.
There are a few improvements that an enthusiastic ITK developer could make to this filter. One issue is with the way that the kernel is applied. The computational expense comes from two sources: computing the kernel weights K(t) and multiplying the pixels in the window by the kernel weights. The first is done more or less efficiently in \( 2 m d \) operations (where d is the dimensionality of the image). The second can be done better. Presently, each pixel \( I(i,j,k) \) is multiplied by the weights \( K(x-i), K(y-j), K(z-k) \) and added to the running total. This results in \( d (2m)^d \) multiplication operations. However, by keeping intermediate sums, it would be possible to do the operation in \( O ( (2m)^d ) \) operations. This would require some creative coding. In addition, in the case when one of the coordinates is integer, the computation could be reduced by an order of magnitude.
Definition at line 271 of file itkWindowedSincInterpolateImageFunction.h.
Public Types | |
using | ConstPointer = SmartPointer< const Self > |
using | ImageType = TInputImage |
using | Pointer = SmartPointer< Self > |
using | Self = WindowedSincInterpolateImageFunction |
using | Superclass = InterpolateImageFunction< TInputImage, TCoordRep > |
Public Types inherited from itk::InterpolateImageFunction< TInputImage, double > | |
using | ConstPointer = SmartPointer< const Self > |
using | Pointer = SmartPointer< Self > |
using | RealType = typename NumericTraits< typename TInputImage::PixelType >::RealType |
using | Self = InterpolateImageFunction |
using | SizeType = typename InputImageType::SizeType |
using | Superclass = ImageFunction< TInputImage, typename NumericTraits< typename TInputImage::PixelType >::RealType, double > |
Public Types inherited from itk::ImageFunction< TInputImage, NumericTraits< TInputImage::PixelType >::RealType, double > | |
using | ConstPointer = SmartPointer< const Self > |
using | ContinuousIndexType = ContinuousIndex< double, Self::ImageDimension > |
using | CoordRepType = double |
using | IndexType = typename InputImageType::IndexType |
using | IndexValueType = typename InputImageType::IndexValueType |
using | InputImageConstPointer = typename InputImageType::ConstPointer |
using | InputImageType = TInputImage |
using | InputPixelType = typename InputImageType::PixelType |
using | OutputType = NumericTraits< TInputImage::PixelType >::RealType |
using | Pointer = SmartPointer< Self > |
using | PointType = Point< double, Self::ImageDimension > |
using | Self = ImageFunction |
using | Superclass = FunctionBase< Point< double, Self::ImageDimension >, NumericTraits< TInputImage::PixelType >::RealType > |
Public Types inherited from itk::FunctionBase< Point< double, TInputImage::ImageDimension >, NumericTraits< TInputImage::PixelType >::RealType > | |
using | ConstPointer = SmartPointer< const Self > |
using | InputType = Point< double, TInputImage::ImageDimension > |
using | OutputType = NumericTraits< TInputImage::PixelType >::RealType |
using | Pointer = SmartPointer< Self > |
using | Self = FunctionBase |
using | Superclass = Object |
Public Types inherited from itk::Object | |
using | ConstPointer = SmartPointer< const Self > |
using | Pointer = SmartPointer< Self > |
using | Self = Object |
using | Superclass = LightObject |
Public Types inherited from itk::LightObject | |
using | ConstPointer = SmartPointer< const Self > |
using | Pointer = SmartPointer< Self > |
using | Self = LightObject |
Public Member Functions | |
OutputType | EvaluateAtContinuousIndex (const ContinuousIndexType &index) const override |
const char * | GetNameOfClass () const override |
SizeType | GetRadius () const override |
void | SetInputImage (const ImageType *image) override |
Public Member Functions inherited from itk::InterpolateImageFunction< TInputImage, double > | |
OutputType | EvaluateAtContinuousIndex (const ContinuousIndexType &index) const override=0 |
OutputType | EvaluateAtIndex (const IndexType &index) const override |
const char * | GetNameOfClass () const override |
virtual SizeType | GetRadius () const=0 |
OutputType | Evaluate (const PointType &point) const override |
Public Member Functions inherited from itk::ImageFunction< TInputImage, NumericTraits< TInputImage::PixelType >::RealType, double > | |
void | ConvertContinuousIndexToNearestIndex (const ContinuousIndexType &cindex, IndexType &index) const |
void | ConvertPointToContinuousIndex (const PointType &point, ContinuousIndexType &cindex) const |
OutputType | Evaluate (const PointType &point) const override=0 |
virtual NumericTraits< TInputImage::PixelType >::RealType | EvaluateAtContinuousIndex (const ContinuousIndexType &index) const=0 |
virtual NumericTraits< TInputImage::PixelType >::RealType | EvaluateAtIndex (const IndexType &index) const=0 |
virtual const ContinuousIndexType & | GetEndContinuousIndex () const |
virtual const IndexType & | GetEndIndex () const |
const InputImageType * | GetInputImage () const |
const char * | GetNameOfClass () const override |
virtual const ContinuousIndexType & | GetStartContinuousIndex () const |
virtual const IndexType & | GetStartIndex () const |
virtual void | SetInputImage (const InputImageType *ptr) |
virtual bool | IsInsideBuffer (const IndexType &index) const |
virtual bool | IsInsideBuffer (const ContinuousIndexType &index) const |
virtual bool | IsInsideBuffer (const PointType &point) const |
void | ConvertPointToNearestIndex (const PointType &point, IndexType &index) const |
Public Member Functions inherited from itk::FunctionBase< Point< double, TInputImage::ImageDimension >, NumericTraits< TInputImage::PixelType >::RealType > | |
virtual OutputType | Evaluate (const InputType &input) const=0 |
const char * | GetNameOfClass () const override |
Public Member Functions inherited from itk::Object | |
unsigned long | AddObserver (const EventObject &event, Command *cmd) const |
unsigned long | AddObserver (const EventObject &event, std::function< void(const EventObject &)> function) const |
LightObject::Pointer | CreateAnother () const override |
virtual void | DebugOff () const |
virtual void | DebugOn () const |
Command * | GetCommand (unsigned long tag) |
bool | GetDebug () const |
MetaDataDictionary & | GetMetaDataDictionary () |
const MetaDataDictionary & | GetMetaDataDictionary () const |
virtual ModifiedTimeType | GetMTime () const |
const char * | GetNameOfClass () const override |
virtual const TimeStamp & | GetTimeStamp () const |
bool | HasObserver (const EventObject &event) const |
void | InvokeEvent (const EventObject &) |
void | InvokeEvent (const EventObject &) const |
virtual void | Modified () const |
void | Register () const override |
void | RemoveAllObservers () |
void | RemoveObserver (unsigned long tag) const |
void | SetDebug (bool debugFlag) const |
void | SetReferenceCount (int) override |
void | UnRegister () const noexcept override |
void | SetMetaDataDictionary (const MetaDataDictionary &rhs) |
void | SetMetaDataDictionary (MetaDataDictionary &&rrhs) |
virtual void | SetObjectName (std::string _arg) |
virtual const std::string & | GetObjectName () const |
Public Member Functions inherited from itk::LightObject | |
Pointer | Clone () const |
virtual Pointer | CreateAnother () const |
virtual void | Delete () |
virtual const char * | GetNameOfClass () const |
virtual int | GetReferenceCount () const |
void | Print (std::ostream &os, Indent indent=0) const |
virtual void | Register () const |
virtual void | SetReferenceCount (int) |
virtual void | UnRegister () const noexcept |
Static Public Member Functions | |
static Pointer | New () |
Static Public Member Functions inherited from itk::Object | |
static bool | GetGlobalWarningDisplay () |
static void | GlobalWarningDisplayOff () |
static void | GlobalWarningDisplayOn () |
static Pointer | New () |
static void | SetGlobalWarningDisplay (bool val) |
Static Public Member Functions inherited from itk::LightObject | |
static void | BreakOnError () |
static Pointer | New () |
Static Public Attributes | |
static constexpr unsigned int | ImageDimension = Superclass::ImageDimension |
Static Public Attributes inherited from itk::InterpolateImageFunction< TInputImage, double > | |
static constexpr unsigned int | ImageDimension |
Static Public Attributes inherited from itk::ImageFunction< TInputImage, NumericTraits< TInputImage::PixelType >::RealType, double > | |
static constexpr unsigned int | ImageDimension |
Protected Member Functions | |
void | PrintSelf (std::ostream &os, Indent indent) const override |
WindowedSincInterpolateImageFunction ()=default | |
~WindowedSincInterpolateImageFunction () override=default | |
Protected Member Functions inherited from itk::InterpolateImageFunction< TInputImage, double > | |
InterpolateImageFunction ()=default | |
void | PrintSelf (std::ostream &os, Indent indent) const override |
~InterpolateImageFunction () override=default | |
Protected Member Functions inherited from itk::ImageFunction< TInputImage, NumericTraits< TInputImage::PixelType >::RealType, double > | |
ImageFunction () | |
void | PrintSelf (std::ostream &os, Indent indent) const override |
~ImageFunction () override=default | |
Protected Member Functions inherited from itk::FunctionBase< Point< double, TInputImage::ImageDimension >, NumericTraits< TInputImage::PixelType >::RealType > | |
FunctionBase ()=default | |
~FunctionBase () override=default | |
Protected Member Functions inherited from itk::Object | |
Object () | |
bool | PrintObservers (std::ostream &os, Indent indent) const |
void | PrintSelf (std::ostream &os, Indent indent) const override |
virtual void | SetTimeStamp (const TimeStamp &timeStamp) |
~Object () override | |
Protected Member Functions inherited from itk::LightObject | |
virtual LightObject::Pointer | InternalClone () const |
LightObject () | |
virtual void | PrintHeader (std::ostream &os, Indent indent) const |
virtual void | PrintSelf (std::ostream &os, Indent indent) const |
virtual void | PrintTrailer (std::ostream &os, Indent indent) const |
virtual | ~LightObject () |
Private Types | |
using | IteratorType = ConstNeighborhoodIterator< ImageType, TBoundaryCondition > |
Private Member Functions | |
double | Sinc (double x) const |
Private Attributes | |
unsigned int | m_OffsetTable [m_OffsetTableSize] {} |
unsigned int | m_WeightOffsetTable [m_OffsetTableSize][ImageDimension] {} |
TWindowFunction | m_WindowFunction {} |
Static Private Attributes | |
static constexpr unsigned int | m_OffsetTableSize = Math::UnsignedPower(m_WindowSize, ImageDimension) |
static constexpr unsigned int | m_WindowSize { 2 * VRadius } |
Additional Inherited Members | |
Protected Attributes inherited from itk::ImageFunction< TInputImage, NumericTraits< TInputImage::PixelType >::RealType, double > | |
ContinuousIndexType | m_EndContinuousIndex |
IndexType | m_EndIndex |
InputImageConstPointer | m_Image |
ContinuousIndexType | m_StartContinuousIndex |
IndexType | m_StartIndex |
Protected Attributes inherited from itk::LightObject | |
std::atomic< int > | m_ReferenceCount {} |
using itk::WindowedSincInterpolateImageFunction< TInputImage, VRadius, TWindowFunction, TBoundaryCondition, TCoordRep >::ConstPointer = SmartPointer<const Self> |
Definition at line 282 of file itkWindowedSincInterpolateImageFunction.h.
using itk::WindowedSincInterpolateImageFunction< TInputImage, VRadius, TWindowFunction, TBoundaryCondition, TCoordRep >::ImageType = TInputImage |
Image type definition
Definition at line 310 of file itkWindowedSincInterpolateImageFunction.h.
|
private |
Definition at line 342 of file itkWindowedSincInterpolateImageFunction.h.
using itk::WindowedSincInterpolateImageFunction< TInputImage, VRadius, TWindowFunction, TBoundaryCondition, TCoordRep >::Pointer = SmartPointer<Self> |
Definition at line 281 of file itkWindowedSincInterpolateImageFunction.h.
using itk::WindowedSincInterpolateImageFunction< TInputImage, VRadius, TWindowFunction, TBoundaryCondition, TCoordRep >::Self = WindowedSincInterpolateImageFunction |
Standard class type aliases.
Definition at line 278 of file itkWindowedSincInterpolateImageFunction.h.
using itk::WindowedSincInterpolateImageFunction< TInputImage, VRadius, TWindowFunction, TBoundaryCondition, TCoordRep >::Superclass = InterpolateImageFunction<TInputImage, TCoordRep> |
Definition at line 279 of file itkWindowedSincInterpolateImageFunction.h.
|
protecteddefault |
|
overrideprotecteddefault |
|
overridevirtual |
Evaluate the function at a ContinuousIndex position
Returns the interpolated image intensity at a specified point position. Bounds checking is based on the type of the TBoundaryCondition specified.
Implements itk::InterpolateImageFunction< TInputImage, double >.
|
overridevirtual |
Reimplemented from itk::Object.
|
inlineoverridevirtual |
Get the radius required for interpolation.
This defines the number of surrounding pixels required to interpolate at a given point.
Implements itk::InterpolateImageFunction< TInputImage, double >.
Definition at line 328 of file itkWindowedSincInterpolateImageFunction.h.
|
static |
Method for creation through the object factory.
|
overrideprotectedvirtual |
Methods invoked by Print() to print information about the object including superclasses. Typically not called by the user (use Print() instead) but used in the hierarchical print process to combine the output of several classes.
Reimplemented from itk::Object.
|
overridevirtual |
Set the input image.
Reimplemented from itk::ImageFunction< TInputImage, NumericTraits< TInputImage::PixelType >::RealType, double >.
|
inlineprivate |
The sinc function
Definition at line 362 of file itkWindowedSincInterpolateImageFunction.h.
References itk::Math::pi.
|
staticconstexpr |
Dimension underlying input image.
Definition at line 300 of file itkWindowedSincInterpolateImageFunction.h.
|
private |
The offset array, used to keep a list of relevant offsets in the neighborhoodIterator
Definition at line 355 of file itkWindowedSincInterpolateImageFunction.h.
|
staticconstexprprivate |
Size of the offset table
Definition at line 351 of file itkWindowedSincInterpolateImageFunction.h.
|
private |
Index into the weights array for each offset
Definition at line 358 of file itkWindowedSincInterpolateImageFunction.h.
|
private |
The function object, used to compute window
Definition at line 348 of file itkWindowedSincInterpolateImageFunction.h.
|
staticconstexprprivate |
Definition at line 345 of file itkWindowedSincInterpolateImageFunction.h.