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
|
