/usr/include/mia-2.4/mia/core/sparse_histogram.hh is in libmia-2.4-dev 2.4.3-5.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 | /* -*- mia-c++ -*-
*
* This file is part of MIA - a toolbox for medical image analysis
* Copyright (c) Leipzig, Madrid 1999-2016 Gert Wollny
*
* MIA is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation; either version 3 of the License, or
* (at your option) any later version.
*
* This program is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with MIA; if not, see <http://www.gnu.org/licenses/>.
*
*/
#include <mia/core/filter.hh>
#include <vector>
#include <type_traits>
#include <limits>
NS_MIA_BEGIN
/**
\ingroup utils
\brief A sparse histogram
This class implements a sparse histogram. As input only
8 and 16 bit (un)signed data is allowed.
*/
class EXPORT_CORE CSparseHistogram : public TFilter<size_t> {
public:
typedef std::vector<std::pair<int, unsigned long>> Compressed;
CSparseHistogram();
/**
Fees pixels to the histogram
\tparam InIterator input iterator, must be a forward iterator
\param begin start of range
\param end end of range
\returns number of added pixels
*/
template <typename InIterator>
size_t operator ()(InIterator begin, InIterator end);
template <typename Image>
size_t operator ()(const Image& image) {
return (*this)(image.begin(), image.end());
}
/**
\returns the histogram as a vector of <value, count> pairs
*/
Compressed get_compressed_histogram()const;
private:
std::vector<unsigned long> m_histogram;
int m_shift;
EPixelType m_pixeltype;
};
EXPORT_CORE std::ostream& operator << (std::ostream& os, const std::pair<short, unsigned long>& pair);
// Implementation
/// @cond INTERNAL
template <typename InIterator, bool sig>
struct dispatch_by_pixeltype {
static size_t apply(InIterator MIA_PARAM_UNUSED(begin), InIterator MIA_PARAM_UNUSED(end),
std::vector<unsigned long>& MIA_PARAM_UNUSED(histogram)){
throw std::invalid_argument("Input pixel type not supported");
}
};
template <typename InIterator>
struct dispatch_by_pixeltype<InIterator, false> {
static size_t apply(InIterator begin, InIterator end, std::vector<unsigned long>& histogram){
size_t n = 0;
while ( begin != end) {
++histogram[*begin];
++begin;
++n;
}
return n;
}
};
template <typename InIterator>
struct dispatch_by_pixeltype<InIterator, true> {
static size_t apply(InIterator begin, InIterator end, std::vector<unsigned long>& histogram){
typedef typename InIterator::value_type in_pixels;
int shift = -std::numeric_limits<in_pixels>::min();
size_t n = 0;
while ( begin != end) {
++histogram[*begin + shift];
++begin;
++n;
}
return n;
}
};
template <typename InIterator>
size_t CSparseHistogram::operator ()(InIterator begin, InIterator end)
{
typedef typename InIterator::value_type in_pixeltype;
if (m_pixeltype ==it_none) {
m_pixeltype = pixel_type<in_pixeltype>::value;
m_shift = -std::numeric_limits<in_pixeltype>::min();
switch (m_pixeltype) {
case it_sbyte:
case it_ubyte:
m_histogram.resize(256);
break;
case it_sshort:
case it_ushort:
m_histogram.resize(65536);
break;
default:
throw create_exception<std::invalid_argument>("Input pixel type '",
CPixelTypeDict.get_name(m_pixeltype),
"' not supported.");
}
} else if (m_pixeltype != pixel_type<in_pixeltype>::value){
throw create_exception<std::invalid_argument>("Input pixels not of consisted type, started with ",
CPixelTypeDict.get_name(m_pixeltype), ", but got now ",
CPixelTypeDict.get_name(pixel_type<in_pixeltype>::value));
}
const bool is_signed = std::is_signed<in_pixeltype>::value;
size_t n = 0;
n += dispatch_by_pixeltype<InIterator, is_signed>::apply(begin, end, m_histogram);
return n;
}
/// @endcond
NS_MIA_END
|