/usr/include/OTB-6.4/otbMRFEnergyFisherClassification.h is in libotb-dev 6.4.0+dfsg-1.
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 | /*
* Copyright (C) 2008 Jan Wegner
*
* This file is part of Orfeo Toolbox
*
* https://www.orfeo-toolbox.org/
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*/
#ifndef otbMRFEnergyFisherClassification_h
#define otbMRFEnergyFisherClassification_h
#include "otbMRFEnergy.h"
#include "otbMath.h"
#include "otb_boost_math_gamma_header.h"
namespace otb
{
/**
* \class MRFEnergyFisherClassification
* \brief Fisher model for Markov classification
*
* This is the implementation of the Fisher model for Markov classification
* to be used for the fidelity (likelihood) term for the classification.
*
* This class is meant to be used in the MRF framework with the
* otb::MarkovRandomFieldFilter
*
*
* \ingroup OTBMarkov
*/
template< class TInput1, class TInput2>
class ITK_EXPORT MRFEnergyFisherClassification: public MRFEnergy< TInput1, TInput2>
{
public:
typedef MRFEnergyFisherClassification Self;
typedef MRFEnergy< TInput1, TInput2> Superclass;
typedef itk::SmartPointer<Self> Pointer;
typedef itk::SmartPointer<const Self> ConstPointer;
typedef TInput1 InputImageType;
typedef TInput2 LabelledImageType;
typedef typename InputImageType::PixelType InputImagePixelType;
typedef typename LabelledImageType::PixelType LabelledImagePixelType;
typedef itk::Array < double > ParametersType;
itkNewMacro(Self);
itkTypeMacro(MRFEnergyFisherClassification, MRFEnergy);
void SetNumberOfParameters(const unsigned int nParameters) ITK_OVERRIDE
{
Superclass::SetNumberOfParameters(nParameters);
this->m_Parameters.SetSize(nParameters);
this->Modified();
}
double GetSingleValue(const InputImagePixelType & value1, const LabelledImagePixelType & value2) ITK_OVERRIDE
{
if ((unsigned int)value2 >= this->GetNumberOfParameters()/3)
{
itkExceptionMacro(<<"Number of parameters does not correspond to number of classes" );
}
double val1 = static_cast<double>(value1);
double mu = this->m_Parameters[3*value2];
double l = this->m_Parameters[3*value2+1];
double m = this->m_Parameters[3*value2+2];
double result = -vcl_log((boost::math::tgamma(l+m)/(boost::math::tgamma(l)*boost::math::tgamma(m)))
* (2/(mu)) * (vcl_sqrt(l/m)) *
((vcl_pow((vcl_sqrt(l/m)*(val1/mu)), ((2*l)-1))) /
(vcl_pow(1+(vcl_sqrt(l/m)*(val1/mu)*vcl_sqrt(l/m)*(val1/mu)), (l+m)))));
return result;
}
protected:
// The constructor and destructor.
MRFEnergyFisherClassification() {};
~MRFEnergyFisherClassification() ITK_OVERRIDE {};
};
}
#endif
|