seriesstats.hh
Go to the documentation of this file.
1 /* -*- mia-c++ -*-
2  *
3  * This file is part of MIA - a toolbox for medical image analysis
4  * Copyright (c) Leipzig, Madrid 1999-2013 Gert Wollny
5  *
6  * MIA is free software; you can redistribute it and/or modify
7  * it under the terms of the GNU General Public License as published by
8  * the Free Software Foundation; either version 3 of the License, or
9  * (at your option) any later version.
10  *
11  * This program is distributed in the hope that it will be useful,
12  * but WITHOUT ANY WARRANTY; without even the implied warranty of
13  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14  * GNU General Public License for more details.
15  *
16  * You should have received a copy of the GNU General Public License
17  * along with MIA; if not, see <http://www.gnu.org/licenses/>.
18  *
19  */
20 
21 #ifndef mia_2d_seriesstats_hh
22 #define mia_2d_seriesstats_hh
23 
24 
25 #include <mia/core/filter.hh>
26 
28 
35  double sum;
37  double sumsq;
39  double mean;
41  double variation;
43  double min;
45  double max;
47  size_t n;
48 };
49 
60 public:
62 
68  template <typename Container>
69  void operator () ( const Container& data);
70 
71 
76  const SIntensityStats& get_result() const;
77 private:
78  mutable SIntensityStats m_stats;
79  mutable bool m_stats_valid;
80 };
81 
82 template <typename Container>
83 void FIntensityStatsAccumulator::operator () ( const Container& data)
84 {
85  m_stats_valid = false;
86  m_stats.n += data.size();
87  for (auto i = data.begin(); i != data.end(); ++i) {
88  m_stats.sum += *i;
89  m_stats.sumsq += *i * *i;
90  if (m_stats.min > *i)
91  m_stats.min = *i;
92  if (m_stats.max < *i)
93  m_stats.max = *i;
94  }
95 }
96 
98 
99 #endif