/usr/include/OTB-5.8/otbVirtualDimensionality.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 | /*=========================================================================
Program: ORFEO Toolbox
Language: C++
Date: $Date$
Version: $Revision$
Copyright (c) Centre National d'Etudes Spatiales. All rights reserved.
See OTBCopyright.txt for details.
Some parts of this code are derived from ITK. See ITKCopyright.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 otbVirtualDimensionality_txx
#define otbVirtualDimensionality_txx
#include "otbVirtualDimensionality.h"
#include <vcl_algorithm.h>
#include "otb_boost_math_normal_header.h"
namespace otb
{
template<class TPrecision>
VirtualDimensionality<TPrecision>
::VirtualDimensionality()
: m_NumberOfPixels(0),
m_NumberOfEndmembers(0),
m_FAR(1.0E-3)
{
}
template<class TInputImage>
void
VirtualDimensionality<TInputImage>
::Compute()
{
// TODO check size
const unsigned int nbBands = m_Covariance.rows();
// Compute diagonalisation of sample covariance and correlation matrices
vnl_symmetric_eigensystem<PrecisionType> eigenK(m_Covariance);
VectorType eigenCovariance = eigenK.D.diagonal();
vcl_sort(eigenCovariance.begin(), eigenCovariance.end());
eigenCovariance.flip();
vnl_symmetric_eigensystem<PrecisionType> eigenR(m_Correlation);
VectorType eigenCorrelation = eigenR.D.diagonal();
vcl_sort(eigenCorrelation.begin(), eigenCorrelation.end());
eigenCorrelation.flip();
m_NumberOfEndmembers = 0;
for(unsigned int i = 0; i < nbBands; ++i)
{
double sigma = vcl_sqrt( 2.0 / m_NumberOfPixels
* (eigenCovariance[i] + eigenCorrelation[i]
+ eigenCovariance[i] * eigenCorrelation[i]) );
boost::math::normal normalDist(0, sigma);
double tau = -boost::math::quantile(normalDist, m_FAR);
if (eigenCorrelation[i] - eigenCovariance[i] > tau )
m_NumberOfEndmembers++;
}
}
template <class TImage>
void
VirtualDimensionality<TImage>
::PrintSelf(std::ostream& os, itk::Indent indent) const
{
Superclass::PrintSelf(os, indent);
os << indent << "Covariance: " << m_Covariance << std::endl;
os << indent << "Correlation: " << m_Correlation << std::endl;
os << indent << "NumberOfEndmembers: " << m_NumberOfEndmembers << std::endl;
}
} // end namespace otb
#endif
|