aimsalgo 6.0.0
Neuroimaging image processing
sixthorderresampler_d.h
Go to the documentation of this file.
1/* This software and supporting documentation are distributed by
2 * Institut Federatif de Recherche 49
3 * CEA/NeuroSpin, Batiment 145,
4 * 91191 Gif-sur-Yvette cedex
5 * France
6 *
7 * This software is governed by the CeCILL-B license under
8 * French law and abiding by the rules of distribution of free software.
9 * You can use, modify and/or redistribute the software under the
10 * terms of the CeCILL-B license as circulated by CEA, CNRS
11 * and INRIA at the following URL "http://www.cecill.info".
12 *
13 * As a counterpart to the access to the source code and rights to copy,
14 * modify and redistribute granted by the license, users are provided only
15 * with a limited warranty and the software's author, the holder of the
16 * economic rights, and the successive licensors have only limited
17 * liability.
18 *
19 * In this respect, the user's attention is drawn to the risks associated
20 * with loading, using, modifying and/or developing or reproducing the
21 * software by the user in light of its specific status of free software,
22 * that may mean that it is complicated to manipulate, and that also
23 * therefore means that it is reserved for developers and experienced
24 * professionals having in-depth computer knowledge. Users are therefore
25 * encouraged to load and test the software's suitability as regards their
26 * requirements in conditions enabling the security of their systems and/or
27 * data to be ensured and, more generally, to use and operate it in the
28 * same conditions as regards security.
29 *
30 * The fact that you are presently reading this means that you have had
31 * knowledge of the CeCILL-B license and that you accept its terms.
32 */
33
34
35#ifndef AIMS_RESAMPLING_SIXTHORDERRESAMPLER_D_H
36#define AIMS_RESAMPLING_SIXTHORDERRESAMPLER_D_H
37
38#include <aims/resampling/sixthorderresampler.h>
39
40#include <cmath>
41
42namespace aims
43{
44
45template < class T >
46SixthOrderResampler< T >::SixthOrderResampler()
47 : SplineResampler< T >()
48{
49
50 this->_poles.push_back(
51 -0.488294589303044755130118038883789062112279161239377608394 );
52 this->_poles.push_back(
53 -0.081679271076237512597937765737059080653379610398148178525368 );
54 this->_poles.push_back(
55 -0.00141415180832581775108724397655859252786416905534669851652709 );
56 this->_gain = 2.598975999348577818390170516255374207847876853191217652822;
57
58}
59
60
61template < class T >
65
66
67template < class T >
68int SixthOrderResampler< T >::getOrder() const
69{
70
71 return 6;
72
73}
74
75
76template < class T >
77double SixthOrderResampler< T >::getBSplineWeight( int i, double x ) const
78{
79
80 x = std::fabs( x - ( double )i );
81 if ( x < 0.5 )
82 {
83
84 x *= x;
85 return x * ( x * ( 7.0 / 48.0 - x * ( 1.0 / 36.0 ) ) - 77.0 / 192.0 ) +
86 5887.0 / 11520.0;
87
88 }
89 if ( x < 1.5 )
90 {
91
92 return x * ( x * ( x * ( x * ( x * ( x * ( 1.0 / 48.0 ) - 7.0 / 48.0 ) +
93 0.328125 ) - 35.0 / 288.0 ) - 91.0 / 256.0 ) - 7.0 / 768.0 ) +
94 7861.0 / 15360.0;
95
96 }
97 if ( x < 2.5 )
98 {
99
100 return x * ( x * ( x * ( x * ( x * ( 7.0 / 60.0 - x * ( 1.0 / 120.0 ) ) -
101 0.65625 ) + 133.0 / 72.0 ) - 2.5703125 ) + 1267.0 / 960.0 ) +
102 1379.0 / 7680.0;
103
104 }
105 if ( x < 3.5 )
106 {
107
108 x -= 3.5;
109 x *= x * x;
110 return x * x * ( 1.0 / 720.0 );
111
112 }
113 return 0.0;
114
115}
116
117} // namespace aims
118
119#endif
aims::SixthOrderResampler::getOrder
int getOrder() const CARTO_OVERRIDE
Spline order (1 to 7)
Definition sixthorderresampler_d.h:68
aims::SixthOrderResampler::getBSplineWeight
double getBSplineWeight(int i, double x) const CARTO_OVERRIDE
Returns .
Definition sixthorderresampler_d.h:77
aims::SplineResampler::_poles
std::vector< double > _poles
Definition splineresampler.h:189