/usr/include/mia-2.2/mia/template/binarize.cxx is in libmia-2.2-dev 2.2.2-1+b1.
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 | /* -*- mia-c++ -*-
*
* This file is part of MIA - a toolbox for medical image analysis
* Copyright (c) Leipzig, Madrid 1999-2014 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 <limits>
#include <stdexcept>
#include <boost/type_traits.hpp>
NS_MIA_BEGIN
template <typename T, bool is_integral>
struct FBinarize {
FBinarize(float fmin, float fmax):
m_min(std::numeric_limits<T>::min()),
m_max(std::numeric_limits<T>::max()){
if (fmin > m_min)
m_min = static_cast<T>(fmin < std::numeric_limits<T>::max() ? fmin : std::numeric_limits<T>::max());
if (fmax < m_max)
m_max = static_cast<T>(fmax > std::numeric_limits<T>::min() ? fmax : std::numeric_limits<T>::min());
}
bool operator ()(T x)const {
return !(x > m_max || x < m_min);
}
private:
T m_min;
T m_max;
};
template <typename T>
struct FBinarize<T, false> {
FBinarize(float fmin, float fmax):
m_min(fmin),
m_max(fmax){
}
bool operator ()(T x)const {
return !(x > m_max || x < m_min);
}
private:
T m_min;
T m_max;
};
template <class Image>
template <template <typename> class Data, typename T>
typename TBinarize<Image>::result_type TBinarize<Image>::operator () (const Data<T>& data) const
{
const bool is_integral = ::boost::is_integral<T>::value;
Data<bool> *result = new Data<bool>(data.get_size(), data);
if (!result)
throw create_exception<std::runtime_error>("binarize: unable to allocate image of size ", data.get_size());
std::transform(data.begin(), data.end(), result->begin(),
FBinarize<T, is_integral>(m_min, m_max));
return result_type(result);
}
/* The actual filter dispacth function calls the filter by selecting the right pixel type through wrap_filter */
template <class Image>
typename TBinarize<Image>::result_type TBinarize<Image>::do_filter(const Image& image) const
{
return mia::filter(*this, image);
}
template <class Image>
TBinarizeImageFilterFactory<Image>::TBinarizeImageFilterFactory():
TDataFilterPlugin<Image>("binarize"),
m_min(0),
m_max(std::numeric_limits<float>::max())
{
this->add_parameter("min", new CFloatParameter(m_min, 0, std::numeric_limits<float>::max(),
false, "minimum of accepted range"));
this->add_parameter("max", new CFloatParameter(m_max, 0, std::numeric_limits<float>::max(),
false, "maximum of accepted range"));
}
template <class Image>
TDataFilter<Image> *TBinarizeImageFilterFactory<Image>::do_create()const
{
return new TBinarize<Image>(m_min, m_max);
}
template <class Image>
const std::string TBinarizeImageFilterFactory<Image>::do_get_descr()const
{
return "image binarize filter";
}
NS_MIA_END
|