A.I.M.S


mask_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 #ifndef AIMS_RESAMPLING_MASK_D_H
35 #define AIMS_RESAMPLING_MASK_D_H
36 
37 #include <aims/resampling/mask.h>
38 
39 namespace aims {
40 
41  template <typename T>
44  const T & val, bool insideval )
45  {
46  BucketMap<Void>::const_iterator ib, eb = src.end();
49  std::vector<float> vs;
50  src.header().getProperty( "voxel_size", vs );
51  while( vs.size() < 4 )
52  vs.push_back( 1. );
53  out->setSizeXYZT( vs[0], vs[1], vs[2], vs[3] );
54  out->setHeader( src.header() ); // maybe a bad idea ?
55  int t;
56  for( ib=src.begin(); ib!=eb; ++ib )
57  {
58  // find timestep
59  t = ib->first;
60  if( t < m->getSizeT() )
61  {
62  const BucketMap<Void>::Bucket & bk = ib->second;
63  BucketMap<Void>::Bucket & bkout = (*out)[ ib->first ];
64 
65  if( insideval )
66  {
67  for( i=bk.begin(), e=bk.end(); i!=e; ++i )
68  {
69  const Point3d & pos = i->first;
70  if( m->at( pos[0], pos[1], pos[2], t ) == val )
71  bkout[ i->first ] = Void();
72  }
73  }
74  else
75  {
76  for( i=bk.begin(), e=bk.end(); i!=e; ++i )
77  {
78  const Point3d & pos = i->first;
79  if( m->at( pos[0], pos[1], pos[2], t ) != val )
80  bkout[ i->first ] = Void();
81  }
82  }
83  }
84  }
85  return out;
86  }
87 
88 }
89 
90 #endif
BucketMap< Void > * maskWithVolume(const BucketMap< Void > &src, const carto::rc_ptr< carto::Volume< T > > m, const T &val=0, bool insideval=false)
mask src object with m, avoiding value val in m if insideval is true, keeping only val otherwise ...
Definition: mask_d.h:42
void setSizeXYZT(float sizex, float sizey, float sizez, float sizet)
sets X,Y,Z and T resolutions of the data
Definition: bucketMap.h:309
virtual bool getProperty(const std::string &, Object &) const
void setHeader(const aims::PythonHeader &hdr)
Definition: bucketMap.h:150
const aims::PythonHeader & header() const
Definition: bucketMap.h:148