/usr/include/OTB-5.8/otbBinaryImageDensityFunction.txx 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 | /*=========================================================================
Program: ORFEO Toolbox
Language: C++
Date: $Date$
Version: $Revision$
Copyright (c) Centre National d'Etudes Spatiales. All rights reserved.
See OTBCopyright.txt for details.
Copyright (c) CS Systemes d'information. All rights reserved.
See CSCopyright.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 otbBinaryImageDensityFunction_txx
#define otbBinaryImageDensityFunction_txx
#include "otbBinaryImageDensityFunction.h"
#include "itkConstNeighborhoodIterator.h"
namespace otb
{
/**
* Constructor
*/
template <class TInputImage, class TCoordRep>
BinaryImageDensityFunction<TInputImage, TCoordRep>
::BinaryImageDensityFunction()
{
m_NeighborhoodRadius.Fill(1);
}
/**
*
*/
template <class TInputImage, class TCoordRep>
void
BinaryImageDensityFunction<TInputImage, TCoordRep>
::PrintSelf(std::ostream& os, itk::Indent indent) const
{
this->Superclass::PrintSelf(os, indent);
os << indent << "NeighborhoodRadius: " << m_NeighborhoodRadius << std::endl;
}
/**
*
*/
template <class TInputImage, class TCoordRep>
typename BinaryImageDensityFunction<TInputImage, TCoordRep>
::RealType
BinaryImageDensityFunction<TInputImage, TCoordRep>
::EvaluateAtIndex(const IndexType& index) const
{
RealType sum;
RealType var;
sum = itk::NumericTraits<RealType>::Zero;
if (!this->GetInputImage())
{
return (itk::NumericTraits<RealType>::max());
}
if (!this->IsInsideBuffer(index))
{
return (itk::NumericTraits<RealType>::max());
}
// Create an N-d neighborhood kernel, using a zeroflux boundary condition
typename InputImageType::SizeType kernelSize = m_NeighborhoodRadius;
itk::ConstNeighborhoodIterator<InputImageType>
it(kernelSize, this->GetInputImage(), this->GetInputImage()->GetBufferedRegion());
// Set the iterator at the desired location
it.SetLocation(index);
// Walk the neighborhood
const unsigned int size = it.Size();
for (unsigned int i = 0; i < size; ++i)
{
const RealType value = static_cast<RealType>(it.GetPixel(i));
sum += value;
}
const RealType num = static_cast<RealType>(size);
var = sum / num;
return var;
}
} // end namespace otb
#endif
|