/usr/include/OTB-5.8/otbWindowedSincInterpolateImageBlackmanFunction.h is in libotb-dev 5.8.0+dfsg-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 | /*=========================================================================
Program: ORFEO Toolbox
Language: C++
Date: $Date$
Version: $Revision$
Copyright (c) Centre National d'Etudes Spatiales. All rights reserved.
See OTBCopyright.txt 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 otbWindowedSincInterpolateImageBlackmanFunction_h
#define otbWindowedSincInterpolateImageBlackmanFunction_h
#include "otbWindowedSincInterpolateImageFunctionBase.h"
#include "vnl/vnl_math.h"
namespace otb
{
namespace Function
{
/**
* \class BlackmanWindowFunction
* \brief Window function for sinc interpolation.
* \f[ w(x) = 0.42 + 0.5 cos(\frac{\pi x}{m}) + 0.08 cos(\frac{2 \pi x}{m}) \f]
* \sa WindowedSincInterpolateImageBlackmanFunction
*
* \ingroup OTBInterpolation
*/
template<class TInput = double, class TOutput = double>
class BlackmanWindowFunction
{
public:
void SetRadius(unsigned int radius)
{
m_Radius = radius;
m_Factor1 = CONST_PI / static_cast<double>(radius);
m_Factor2 = 2.0 * CONST_PI / static_cast<double>(radius);
}
unsigned int GetRadius() const
{
return m_Radius;
}
double GetFactor1()
{
return m_Factor1;
}
double GetFactor2()
{
return m_Factor2;
}
inline TOutput operator ()(const TInput& A) const
{
double x = static_cast<double>(A);
double px = CONST_PI * x;
double temp = 0.42 + 0.5 * vcl_cos(x * m_Factor1) + 0.08 * vcl_cos(x * m_Factor2);
return (x == 0.0) ? static_cast<TOutput>(temp) : static_cast<TOutput>(temp * vcl_sin(px) / px);
}
private:
unsigned int m_Radius;
// Equal to \f$ \frac{\pi}{m} \f$
double m_Factor1;
// Equal to \f$ \frac{2 \pi}{m} \f$
double m_Factor2;
};
} //namespace Function
/**
* \class WindowedSincInterpolateImageBlackmanFunction
* \brief Use the WindowedSincInterpolateImageFunctionBase with a Blackman Function.
*
* \sa GenericInterpolatorImageFunction
* \sa Function::GaussianWindowFunction
* \sa Function::HammingWindowFunction
* \sa Function::CosineWindowFunction
* \sa Function::WelchWindowFunction
* \sa Function::LanczosWindowFunction
* \sa Function::BlackmanWindowFunction
* \ingroup ImageFunctionBases ImageInterpolators
*
* \ingroup OTBInterpolation
*/
template<class TInputImage, class TBoundaryCondition = itk::ConstantBoundaryCondition<TInputImage>, class TCoordRep =
double, class TInputInterpolator = double, class TOutputInterpolator = double>
class ITK_EXPORT WindowedSincInterpolateImageBlackmanFunction :
public WindowedSincInterpolateImageFunctionBase<TInputImage,
typename Function::BlackmanWindowFunction<TInputInterpolator,
TOutputInterpolator>,
TBoundaryCondition,
TCoordRep>
{
public:
/** Standard class typedefs. */
typedef WindowedSincInterpolateImageBlackmanFunction Self;
typedef WindowedSincInterpolateImageFunctionBase<TInputImage,
typename Function::BlackmanWindowFunction<TInputInterpolator,
TOutputInterpolator>,
TBoundaryCondition,
TCoordRep>
Superclass;
typedef itk::SmartPointer<Self> Pointer;
typedef itk::SmartPointer<const Self> ConstPointer;
/** Run-time type information (and related methods). */
itkTypeMacro(WindowedSincInterpolateImageBlackmanFunction, WindowedSincInterpolateImageFunctionBase);
/** Method for creation through the object factory. */
itkNewMacro(Self);
/** Input and output images typedef definition. */
typedef typename Superclass::InputImageType InputImageType;
typedef typename Superclass::OutputType OutputType;
/** Dimension underlying input image. */
itkStaticConstMacro(ImageDimension, unsigned int, Superclass::ImageDimension);
/** Superclass typedef inheritance. */
typedef typename Superclass::IndexType IndexType;
typedef typename Superclass::SizeType SizeType;
typedef typename Superclass::RealType RealType;
typedef typename Superclass::IteratorType IteratorType;
typedef typename Superclass::ContinuousIndexType ContinuousIndexType;
protected:
WindowedSincInterpolateImageBlackmanFunction() {};
~WindowedSincInterpolateImageBlackmanFunction() ITK_OVERRIDE {}
void PrintSelf(std::ostream& os, itk::Indent indent) const ITK_OVERRIDE
{
Superclass::PrintSelf(os, indent);
}
private:
WindowedSincInterpolateImageBlackmanFunction(const Self &); //purposely not implemented
void operator =(const Self&); //purposely not implemented
};
} // end namespace otb
#endif
|