approxgm.hh

Go to the documentation of this file.

/* -*- 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/>.
 *
 */


#ifndef mia_core_approxgm_hh
#define mia_core_approxgm_hh

#include <mia/core/minimizer.hh>

#include <vector>
#include <ostream>

NS_MIA_BEGIN

class EXPORT_CORE CApproximateGausianMixture {
public:
        struct EXPORT_CORE Result {
                Result(const std::vector<double>& centers);
                Result(const std::vector<double>& centers,
                       const std::vector<double>& weights,
                       double sigma);

                Result(const CDoubleVector& x); 
                
                Result(const Result& orig) = default;
                Result& operator = (const Result& orig) = default; 

                float value(unsigned x) const;

                void print(std::ostream& os) const; 
                
                std::vector<double> m_centers;
                std::vector<double> m_alpha;
                double m_sigma; 
        }; 
        
        CApproximateGausianMixture(double epsilon, unsigned max_iter);
        ~CApproximateGausianMixture();

        
        Result run(const std::vector<float>& data, const std::vector<double>& init_centers) const;

        template <typename InIterator>
        Result run(InIterator begin, InIterator end, const std::vector<double>& init_centers) const;

private:
        struct CApproximateGausianMixtureImpl *impl; 
}; 

std::ostream& operator << (std::ostream& os, const CApproximateGausianMixture::Result& r)
{
        r.print(os);
        return os; 
}

template <typename InIterator>
CApproximateGausianMixture::Result
CApproximateGausianMixture::run(InIterator begin, InIterator end,
                                const std::vector<double>& init_centers) const
{
        std::vector<float> data(begin, end);
        return run(data, init_centers); 
}

NS_MIA_END
#endif