/usr/include/OTB-5.8/otbContinuousMinimumMaximumImageCalculator.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 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 | /*=========================================================================
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 otbContinuousMinimumMaximumImageCalculator_h
#define otbContinuousMinimumMaximumImageCalculator_h
#include "itkObject.h"
#include "itkObjectFactory.h"
#include "itkNumericTraits.h"
namespace otb
{
/** \class ContinuousMinimumMaximumImageCalculator
* \brief Compute the continuous maximum and the minimum of an image
*
* To compute the continuous maximum and the minimum of an image, this
* class uses a simple polynomial interpolation.
*
* First traditional discrete maximum and minimum are found. Then a second
* order polynomial is fitted between these extrema and their neighboring
* pixels.
*
* The continuous extrema is assumed to be at the zero of the first order
* derivative of this polynom.
*
* If we denote \f$ (x_0, y_0) \f$ the extrema and \f$ (x_{-1}, y_{-1}) \f$ and
* \f$ (x_1, y_1) \f$ its neighbor, the second degree polynom verify the
* following equations:
*
*
* \f[
* y_{-1} = a*x_{-1}^2 + b*x_{-1} +c
* y_0 = a*x_0^2 + b*x_0 +c
* y_1 = a*x_1^2 + b*x_1 +c
* \f]
*
* The maximum is at \f$ -b/2a \f$ with is
* \f$ -\frac{(y_1-y_{-1})}{(2*(y_{-1}+y_1-2*y_0))} \f$
*
* Remark: image is assumed to be 2 dimensionnal
*
*
*
* \ingroup OTBStatistics
*/
template <class TInputImage>
class ITK_EXPORT ContinuousMinimumMaximumImageCalculator :
public itk::Object
{
public:
/** Standard class typedefs. */
typedef ContinuousMinimumMaximumImageCalculator Self;
typedef itk::Object Superclass;
typedef itk::SmartPointer<Self> Pointer;
typedef itk::SmartPointer<const Self> ConstPointer;
/** Method for creation through the object factory. */
itkNewMacro(Self);
/** Run-time type information (and related methods). */
itkTypeMacro(ContinuousMinimumMaximumImageCalculator, Object);
/** Type definition for the input image. */
typedef TInputImage ImageType;
/** Pointer type for the image. */
typedef typename TInputImage::Pointer ImagePointer;
/** Const Pointer type for the image. */
typedef typename TInputImage::ConstPointer ImageConstPointer;
/** Type definition for the input image pixel type. */
typedef typename TInputImage::PixelType PixelType;
/** Type definition for the input image real pixel type. */
typedef typename itk::NumericTraits<PixelType>::RealType RealPixelType;
/** Type definition for the input image index type. */
typedef typename TInputImage::IndexType IndexType;
/** Type definition for the input image index type. */
typedef typename TInputImage::PointType ContinuousIndexType;
/** Type definition for the input image region type. */
typedef typename TInputImage::RegionType RegionType;
/** Set the input image. */
itkSetConstObjectMacro(Image, ImageType);
/** Compute the minimum value of intensity of the input image. */
void ComputeMinimum(void);
/** Compute the maximum value of intensity of the input image. */
void ComputeMaximum(void);
/** Compute the minimum and maximum values of intensity of the input image. */
void Compute(void);
/** Return the minimum intensity value. */
itkGetMacro(Minimum, PixelType);
/** Return the maximum intensity value. */
itkGetMacro(Maximum, PixelType);
/** Return the index of the minimum intensity value. */
itkGetConstReferenceMacro(IndexOfMinimum, IndexType);
/** Return the index of the maximum intensity value. */
itkGetConstReferenceMacro(IndexOfMaximum, IndexType);
/** Return the index of the minimum intensity value. */
itkGetConstReferenceMacro(ContinuousIndexOfMinimum, ContinuousIndexType);
/** Return the index of the maximum intensity value. */
itkGetConstReferenceMacro(ContinuousIndexOfMaximum, ContinuousIndexType);
/** Set the region over which the values will be computed */
void SetRegion(const RegionType& region);
protected:
ContinuousMinimumMaximumImageCalculator();
~ContinuousMinimumMaximumImageCalculator() ITK_OVERRIDE {}
void PrintSelf(std::ostream& os, itk::Indent indent) const ITK_OVERRIDE;
PixelType m_Minimum;
PixelType m_Maximum;
ImageConstPointer m_ContinuousImage;
IndexType m_IndexOfMinimum;
IndexType m_IndexOfMaximum;
ContinuousIndexType m_ContinuousIndexOfMinimum;
ContinuousIndexType m_ContinuousIndexOfMaximum;
ImageConstPointer m_Image;
RegionType m_Region;
bool m_RegionSetByUser;
private:
ContinuousMinimumMaximumImageCalculator(const Self &); //purposely not implemented
void operator =(const Self&); //purposely not implemented
};
} // end namespace otb
#ifndef OTB_MANUAL_INSTANTIATION
#include "otbContinuousMinimumMaximumImageCalculator.txx"
#endif
#endif
|