A.I.M.S algorithms


harmonicCageResampler.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_HARMONICCAGERESAMPLER_H
36 #define AIMS_RESAMPLING_HARMONICCAGERESAMPLER_H
37 
38 #include <aims/data/data.h>
39 #include <aims/mesh/surface.h>
40 #include <aims/distancemap/stlsort.h> //for Point3dfCompare
41 
42 namespace aims
43 {
44 
46 {
47 public:
49  float spacing);
51  std::vector<Point3df> &controls,
52  short background=0, short border=1, short inside=2);
54 
55 public:
56  const Point3df &getControl(int ind) const;
57  const AimsData<float> &get_image_coords(int ind) const;
58  void set_keep_image_coords(bool state)
59  { _keep_image_coords = state; }
60 
61 public:
62  void moveControl(int ind, const Point3df &p);
64  Point3df coordinate(const Point3df &p) const;
67  void computeCoordinates(float threshold=10e-4);
68 
69 protected:
71  std::vector<Point3df> _controls;
72  short _background;
73  short _border;
74  short _inside;
76  std::vector<AimsData<float> > _image_coords;
77  std::map<Point3df, std::map<int, float>, Point3dfCompare > _weights;
78 };
79 
80 
81 inline
83 {
84  int ind;
85  Point3df p2(rint(p[0]), rint(p[1]), rint(p[2]));
86  const std::map<int, float> &w = (*(_weights.find(p2))).second;
87  std::map<int, float>::const_iterator i, e;
88 
89  Point3df res(0, 0, 0);
90  for (i = w.begin(), e = w.end(); i != e; ++i)
91  {
92  ind = (*i).first;
93  res += (*i).second * _controls[ind];
94  }
95  return res;
96 }
97 
98 
99 inline const AimsData<float> &
101 {
102  return _image_coords[ind];
103 }
104 
105 };
106 
107 #endif
108 
HarmonicCageMeshResampler(const AimsSurfaceTriangle &cage, float spacing)
void computeCoordinates(float threshold=10e-4)
AIMSDATA_API AimsTimeSurface< 3, Void > AimsSurfaceTriangle
const Point3df & getControl(int ind) const
Point3df coordinate(const Point3df &p) const
from native coordinate to resample coordinates (point version)
AimsSurfaceTriangle * coordinates(const AimsSurfaceTriangle &mesh) const
from native coordinate to resample coordinates (mesh version)
std::vector< AimsData< float > > _image_coords
const AimsData< float > & get_image_coords(int ind) const
std::map< Point3df, std::map< int, float >, Point3dfCompare > _weights
void moveControl(int ind, const Point3df &p)