[1339] | 1 | #include <duchamp/ATrous/TriangleFilter.hh> |
---|
[1343] | 2 | #include <duchamp/ATrous/filter.hh> |
---|
[1339] | 3 | |
---|
| 4 | namespace duchamp { |
---|
| 5 | |
---|
| 6 | TriangleFilter::TriangleFilter(): |
---|
| 7 | Filter() |
---|
| 8 | { |
---|
| 9 | const unsigned int filtersize=3; |
---|
| 10 | double filter[filtersize] = {1./4., 1./2., 1./4.}; |
---|
| 11 | this->filter1D = std::vector<double>(filter,filter+filtersize); |
---|
| 12 | this->name = "Triangle function"; |
---|
| 13 | |
---|
| 14 | const unsigned int size1D=18; |
---|
| 15 | this->maxNumScales[0] = size1D; |
---|
| 16 | double sigmaFactors1D[size1D+1] = |
---|
| 17 | {1.00000000000,6.12372436e-1,3.30718914e-1,2.11947812e-1, |
---|
| 18 | 1.45740298e-1,1.02310944e-1,7.22128185e-2,5.10388224e-2, |
---|
| 19 | 3.60857673e-2,2.55157615e-2,1.80422389e-2,1.27577667e-2, |
---|
| 20 | 9.02109930e-3,6.37887978e-3,4.51054902e-3,3.18942978e-3, |
---|
| 21 | 2.25527449e-3,1.59471988e-3,1.12763724e-4}; |
---|
| 22 | |
---|
| 23 | const unsigned int size2D=12; |
---|
| 24 | this->maxNumScales[1] = size2D; |
---|
| 25 | double sigmaFactors2D[size2D+1] = |
---|
| 26 | {1.00000000000,8.00390530e-1,2.72878894e-1,1.19779282e-1, |
---|
| 27 | 5.77664785e-2,2.86163283e-2,1.42747506e-2,7.13319703e-3, |
---|
| 28 | 3.56607618e-3,1.78297280e-3,8.91478237e-4,4.45738098e-4, |
---|
| 29 | 2.22868922e-4}; |
---|
| 30 | |
---|
| 31 | const unsigned int size3D=12; |
---|
| 32 | this->maxNumScales[2] = size3D; |
---|
| 33 | double sigmaFactors3D[size3D+1] = |
---|
| 34 | {1.00000000000,8.959544490e-1,1.92033014e-1,5.76484078e-2, |
---|
| 35 | 1.94912393e-2,6.812783870e-3,2.40175885e-3,8.48538128e-4, |
---|
| 36 | 2.99949455e-4}; |
---|
| 37 | |
---|
| 38 | this->sigmaFactors[0] = new std::vector<double>(sigmaFactors1D,sigmaFactors1D+this->maxNumScales[0]); |
---|
| 39 | this->sigmaFactors[1] = new std::vector<double>(sigmaFactors2D,sigmaFactors2D+this->maxNumScales[1]); |
---|
| 40 | this->sigmaFactors[2] = new std::vector<double>(sigmaFactors3D,sigmaFactors3D+this->maxNumScales[2]); |
---|
| 41 | |
---|
| 42 | } |
---|
| 43 | |
---|
| 44 | TriangleFilter::TriangleFilter(const TriangleFilter& other): |
---|
| 45 | Filter(other) |
---|
| 46 | { |
---|
| 47 | } |
---|
| 48 | |
---|
| 49 | TriangleFilter& TriangleFilter::operator= (const TriangleFilter& other) |
---|
| 50 | { |
---|
| 51 | if(this == &other) return *this; |
---|
| 52 | ((Filter &) *this) = other; |
---|
| 53 | return *this; |
---|
| 54 | } |
---|
| 55 | |
---|
| 56 | } |
---|