groundtruthproblem.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_2d_ground_truth_problem_hh
#define mia_2d_ground_truth_problem_hh

#include <mia/core/minimizer.hh>
#include <mia/2d/correlation_weight.hh>


NS_MIA_BEGIN

class  EXPORT_2D GroundTruthProblem : public CMinimizer::Problem {
public:

        GroundTruthProblem(double a, double b, 
                           const C2DBounds& slice_size, 
                           size_t nframes, 
                           const CDoubleVector& left_side, 
                           const  CCorrelationEvaluator::result_type& corr); 

        void set_alpha_beta(double a, double b); 

protected: 
        
        const std::vector<double>& get_spacial_gradient() const; 

        const std::vector<double>& get_time_derivative() const; 
private: 
        virtual double  do_f(const CDoubleVector& x); 
        virtual void    do_df(const CDoubleVector&  x, CDoubleVector&  g); 
        virtual double  do_fdf(const CDoubleVector&  x, CDoubleVector&  g); 
        size_t do_size() const; 
        
        double evaluate_spacial_gradients(const CDoubleVector& x); 
        double  evaluate_time_gradients(const CDoubleVector& x); 
        double evaluate_slice_gradient(CDoubleVector::const_iterator ii,  std::vector<double>::iterator iout); 

        std::vector<double> m_spacial_gradient; 
        std::vector<double> m_time_derivative;
        
        double m_a;
        double m_b;
        const C2DBounds& m_slice_size; 
        size_t m_nframes; 
        size_t m_frame_size; 
        const CDoubleVector& m_left_side; 
        const  CCorrelationEvaluator::result_type& m_corr; 
}; 

NS_MIA_END

#endif