source: trunk/src/ATrous/TriangleFilter.cc @ 1441

#include <duchamp/ATrous/TriangleFilter.hh>
#include <duchamp/ATrous/filter.hh>

namespace duchamp {

    TriangleFilter::TriangleFilter():
    Filter()
    {
        const unsigned int filtersize=3;
        double filter[filtersize] = {1./4., 1./2., 1./4.};
        this->filter1D = std::vector<double>(filter,filter+filtersize);
        this->name = "Triangle function";
        
        const unsigned int size1D=18;
        this->maxNumScales[0] = size1D;
        double sigmaFactors1D[size1D+1] = 
            {1.00000000000,6.12372436e-1,3.30718914e-1,2.11947812e-1,
             1.45740298e-1,1.02310944e-1,7.22128185e-2,5.10388224e-2,
             3.60857673e-2,2.55157615e-2,1.80422389e-2,1.27577667e-2,
             9.02109930e-3,6.37887978e-3,4.51054902e-3,3.18942978e-3,
             2.25527449e-3,1.59471988e-3,1.12763724e-4};

        const unsigned int size2D=12;
        this->maxNumScales[1] = size2D;
        double sigmaFactors2D[size2D+1] = 
            {1.00000000000,8.00390530e-1,2.72878894e-1,1.19779282e-1,
             5.77664785e-2,2.86163283e-2,1.42747506e-2,7.13319703e-3,
             3.56607618e-3,1.78297280e-3,8.91478237e-4,4.45738098e-4,
             2.22868922e-4};

        const unsigned int size3D=12;
        this->maxNumScales[2] = size3D;
        double sigmaFactors3D[size3D+1] = 
            {1.00000000000,8.959544490e-1,1.92033014e-1,5.76484078e-2,
             1.94912393e-2,6.812783870e-3,2.40175885e-3,8.48538128e-4,
             2.99949455e-4};

        this->sigmaFactors[0] = new std::vector<double>(sigmaFactors1D,sigmaFactors1D+this->maxNumScales[0]);
        this->sigmaFactors[1] = new std::vector<double>(sigmaFactors2D,sigmaFactors2D+this->maxNumScales[1]);
        this->sigmaFactors[2] = new std::vector<double>(sigmaFactors3D,sigmaFactors3D+this->maxNumScales[2]);

    }

    TriangleFilter::TriangleFilter(const TriangleFilter& other):
    Filter(other)
    {
    }

    TriangleFilter& TriangleFilter::operator= (const TriangleFilter& other)
    {
        if(this == &other) return *this;
        ((Filter &) *this) = other;
        return *this;
    }

}