spacial_kernel.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-2017 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_spacial_kernel_hh
#define mia_core_spacial_kernel_hh


#include <vector>

// MIA specific
#include <mia/core/defines.hh>
#include <mia/core/factory.hh>

NS_MIA_BEGIN

struct EXPORT_CORE spacial_kernel_data {
        static const char *data_descr;
};


struct EXPORT_CORE kernel_plugin_type {
        static const char *type_descr;
};

class EXPORT_CORE C1DFilterKernel : public CProductBase {
protected:
        
        C1DFilterKernel(unsigned fsize);

        virtual ~C1DFilterKernel();
public:
        typedef spacial_kernel_data plugin_data; 
        typedef kernel_plugin_type  plugin_type; 

        int get_fsize() const;

        size_t size() const;

        void apply_inplace(std::vector<double>& data) const;

        std::vector<double> apply(const std::vector<double>& data) const;

private:
        virtual size_t do_size() const = 0;

        virtual std::vector<double> do_apply(const std::vector<double>& data) const = 0;

        unsigned m_fsize;

};

class EXPORT_CORE C1DFoldingKernel : public C1DFilterKernel {
public:

        
        typedef std::vector<double> vec_mask;

        typedef vec_mask::const_iterator const_iterator;

        C1DFoldingKernel(int fsize);


        const_iterator begin()const;

        const_iterator end()const;

        const_iterator dbegin()const;

        const_iterator dend()const;

        double& operator[](int i) {
                return m_mask[i];
        }
        
        double operator[](int i)const {
                return m_mask[i];
        }
protected:
        typedef vec_mask::iterator iterator;

        iterator begin();

        iterator end();

        iterator dbegin();

        iterator dend();

private:
        virtual size_t do_size() const;

        vec_mask m_mask;
        vec_mask m_derivative;

};

typedef TFactory<C1DFoldingKernel> C1DSpacialKernelPlugin;

typedef std::shared_ptr<C1DFoldingKernel> P1DSpacialKernel; 

template<>  const char * const 
TPluginHandler<TFactory<C1DFoldingKernel>>::m_help; 

typedef THandlerSingleton<TFactoryPluginHandler<C1DSpacialKernelPlugin> > C1DSpacialKernelPluginHandler;


inline P1DSpacialKernel produce_spacial_kernel(const std::string& descr) 
{
        return C1DSpacialKernelPluginHandler::instance().produce(descr); 
}

FACTORY_TRAIT(C1DSpacialKernelPluginHandler); 

NS_MIA_END

#endif