A.I.M.S algorithms


medianfilter.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_SIGNALFILTER_MEDIANFILTER_H
36 #define AIMS_SIGNALFILTER_MEDIANFILTER_H
37 
41 
42 namespace aims {
43 
44  template <typename T>
46  {
47  public:
49  MedianFilter( const StructuringElement & se = strel::Cube(1.0),
50  const carto::Object & options = carto::none() ):
51  ElementFilteringImageAlgorithm<T>( MedianFilterFunc<ChannelType>( options ),
52  se )
53  {}
54  MedianFilter( const carto::Object & options ):
55  ElementFilteringImageAlgorithm<T>( MedianFilterFunc<ChannelType>( options ),
56  strel::Cube(1.0) )
57  {}
58  MedianFilter( const MedianFilter<T> & other ):
60  {}
63  {
65  return *this;
66  }
68  return new MedianFilter<T>(*this); }
69  };
70 
71  template <typename T>
73  {
74  public:
76  NotNullMedianFilter( const StructuringElement & se = strel::Cube(1.0),
77  const carto::Object & options = carto::none() ):
78  ElementFilteringImageAlgorithm<T>( NotNullMedianFilterFunc<ChannelType>( options ),
79  se )
80  {}
81  NotNullMedianFilter( const carto::Object & options ):
82  ElementFilteringImageAlgorithm<T>( NotNullMedianFilterFunc<ChannelType>( options ),
83  strel::Cube(1.0) )
84  {}
87  {}
90  {
92  return *this;
93  }
94  NotNullMedianFilter<T> *clone() const { return new NotNullMedianFilter<T>(*this); }
95  };
96 
97 } // namespace aims
98 
99 //============================================================================
100 // Backward compatibility bindings
101 //============================================================================
102 
103 #include <aims/data/data_g.h>
104 #include <vector>
105 
106 template <typename T>
108 {
109  public:
110  MedianSmoothing( int sx = 3, int sy = 3, int sz = 3 );
111  virtual ~MedianSmoothing();
112  virtual AimsData<T> doit( const AimsData<T>& in ) const;
113  private:
114  MedianSmoothing<T> & operator = ( const MedianSmoothing<T> & );
115  int _sx;
116  int _sy;
117  int _sz;
118 };
119 
120 //----------------------------------------------------------------------------
121 // DEFINITIONS
122 //----------------------------------------------------------------------------
123 
124 template <typename T>
125 MedianSmoothing<T>::MedianSmoothing( int sx, int sy, int sz ):
126  _sx(sx), _sy(sy), _sz(sz)
127 {}
128 
129 template <typename T>
131 {}
132 
133 template <typename T>
135  const MedianSmoothing<T> & other
136 )
137 {
138  _sx = other._sx;
139  _sy = other._sy;
140  _sz = other._sz;
141  return (*this);
142 }
143 
144 template <typename T>
146 {
147  std::vector<double> amplitude(3,0.);
148  amplitude[0] = .5 * (double)_sx;
149  amplitude[1] = .5 * (double)_sy;
150  amplitude[2] = .5 * (double)_sz;
151  aims::strel::Cube se( amplitude, true );
152  aims::MedianFilter<T> f( se );
153  return f.execute( in );
154 }
155 
156 
157 #endif
NotNullMedianFilter< T > & operator=(const NotNullMedianFilter< T > &other)
Definition: medianfilter.h:89
virtual AimsData< T > doit(const AimsData< T > &in) const
Definition: medianfilter.h:145
virtual ~MedianSmoothing()
Definition: medianfilter.h:130
NotNullMedianFilter(const StructuringElement &se=strel::Cube(1.0), const carto::Object &options=carto::none())
Definition: medianfilter.h:76
MedianSmoothing(int sx=3, int sy=3, int sz=3)
Definition: medianfilter.h:125
MedianFilter(const carto::Object &options)
Definition: medianfilter.h:54
carto::DataTypeTraits< T >::ChannelType ChannelType
Definition: medianfilter.h:48
MedianFilter(const MedianFilter< T > &other)
Definition: medianfilter.h:58
carto::DataTypeTraits< T >::ChannelType ChannelType
Definition: medianfilter.h:75
MedianFilter< T > * clone() const
ImageAlgorithmInterface Pure virtual method.
Definition: medianfilter.h:67
NotNullMedianFilter(const carto::Object &options)
Definition: medianfilter.h:81
ElementFilteringImageAlgorithm & operator=(ElementFilteringImageAlgorithm &other)
virtual void execute(const carto::VolumeRef< T > &in, carto::VolumeRef< T > &out) const
ImageAlgorithmInterface Pure virtual method.
MedianFilter< T > & operator=(const MedianFilter< T > &other)
Definition: medianfilter.h:62
NotNullMedianFilter< T > * clone() const
ImageAlgorithmInterface Pure virtual method.
Definition: medianfilter.h:94
Object none()
MedianFilter(const StructuringElement &se=strel::Cube(1.0), const carto::Object &options=carto::none())
Definition: medianfilter.h:49