/usr/include/OTB-5.8/otbWaveletOperatorBase.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 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 166 167 168 169 | /*=========================================================================
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) Institut Mines-Telecom. All rights reserved.
See IMTCopyright.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 otbWaveletOperatorBase_txx
#define otbWaveletOperatorBase_txx
#include "otbWaveletOperatorBase.h"
#include <cassert>
namespace otb {
template <Wavelet::Wavelet TMotherWaveletOperator,
class TPixel, unsigned int VDimension, class TAllocator>
void
WaveletOperatorBase<TMotherWaveletOperator, TPixel, VDimension, TAllocator>
::PrintSelf(std::ostream& os, itk::Indent i) const
{
os << i << "Wavelet kind : " << GetWaveletName() << "\n";
os << i << "Up-Sampling factor " << GetUpSampleFactor() << "\n";
Superclass::PrintSelf(os, i.GetNextIndent());
os << i << "Wavelet coeff: [ ";
for (typename Superclass::ConstIterator iter = Superclass::Begin();
iter != Superclass::End();
++iter)
{
os << *iter << ' ';
}
os << "]\n";
}
template <Wavelet::Wavelet TMotherWaveletOperator,
class TPixel, unsigned int VDimension, class TAllocator>
void
WaveletOperatorBase<TMotherWaveletOperator, TPixel, VDimension, TAllocator>
::UpSamplingCoefficients(CoefficientVector& coeff)
{
if (m_UpSampleFactor <= 1) return;
unsigned long radius = static_cast<unsigned long>(coeff.size()) >> 1;
unsigned long upSampleRadius = radius * m_UpSampleFactor;
CoefficientVector upSampledCoeff ( 2 * upSampleRadius + 1 );
upSampledCoeff.assign(2 * upSampleRadius + 1, 0.);
upSampledCoeff[upSampleRadius] = coeff[radius];
for (unsigned int i = 1; i <= radius; ++i)
{
upSampledCoeff[upSampleRadius + m_UpSampleFactor * i] = coeff[radius + i];
upSampledCoeff[upSampleRadius - m_UpSampleFactor * i] = coeff[radius - i];
}
coeff = upSampledCoeff;
}
template <Wavelet::Wavelet TMotherWaveletOperator,
class TPixel, unsigned int VDimension, class TAllocator>
void
WaveletOperatorBase<TMotherWaveletOperator, TPixel, VDimension, TAllocator>
::RevertFilter(CoefficientVector& coeff)
{
const unsigned int length = coeff.size();
const unsigned int medianPosition = length / 2;
CoefficientVector newCoeff(length);
newCoeff[medianPosition] = coeff[medianPosition];
for (unsigned int i = 1; i <= medianPosition; ++i)
{
newCoeff[medianPosition + i] = coeff[medianPosition - i];
newCoeff[medianPosition - i] = coeff[medianPosition + i];
}
coeff = newCoeff;
}
template <Wavelet::Wavelet TMotherWaveletOperator,
class TPixel, unsigned int VDimension, class TAllocator>
void
WaveletOperatorBase<TMotherWaveletOperator, TPixel, VDimension, TAllocator>
::GenerateInverseHighPassFilterFromLowPassFilter(CoefficientVector& coeff)
{
const unsigned int length = coeff.size();
CoefficientVector highPassCoeff(length + 2);
double sign = -1;
for (unsigned int i = 0; i < length; ++i)
{
highPassCoeff[i] = sign * coeff[i];
sign *= -1.;
}
highPassCoeff[length] = 0.;
highPassCoeff[length + 1] = 0.;
coeff = highPassCoeff;
// Note that the 0.0 value is obtained by affectation and not by
// computation. It is fine to do an == comparison
while ((coeff[0] == coeff[coeff.size() - 1]) && (coeff[0] == 0.0))
{
ReduceFilterLength(coeff);
}
}
template <Wavelet::Wavelet TMotherWaveletOperator,
class TPixel, unsigned int VDimension, class TAllocator>
void
WaveletOperatorBase<TMotherWaveletOperator, TPixel, VDimension, TAllocator>
::GenerateInverseLowPassFilterFromHighPassFilter(CoefficientVector& coeff)
{
const unsigned int length = coeff.size();
CoefficientVector highPassCoeff(length + 2);
double sign = 1;
for (unsigned int i = 0; i < length; ++i)
{
highPassCoeff[i] = sign * coeff[i];
sign *= -1.;
}
highPassCoeff[length] = 0.;
highPassCoeff[length + 1] = 0.;
coeff = highPassCoeff;
// Note that the 0.0 value is obtained by affectation and not by
// computation. It is fine to do an == comparison
while ((coeff[0] == coeff[coeff.size() - 1]) && (coeff[0] == 0.0))
{
assert(coeff.size() > 1);
ReduceFilterLength(coeff);
}
}
template <Wavelet::Wavelet TMotherWaveletOperator,
class TPixel, unsigned int VDimension, class TAllocator>
void
WaveletOperatorBase<TMotherWaveletOperator, TPixel, VDimension, TAllocator>
::ReduceFilterLength(CoefficientVector& coeff)
{
const unsigned int length = coeff.size();
assert(length >= 2);
CoefficientVector newFilter(length - 2);
for (unsigned int i = 0; i < newFilter.size(); ++i)
{
newFilter[i] = coeff[i + 1];
}
coeff = newFilter;
}
} // end of namespace otb
#endif
|