2d/splinetransformpenalty.hh
Go to the documentation of this file.
1 /* -*- mia-c++ -*-
2  *
3  * This file is part of MIA - a toolbox for medical image analysis
4  * Copyright (c) Leipzig, Madrid 1999-2013 Gert Wollny
5  *
6  * MIA is free software; you can redistribute it and/or modify
7  * it under the terms of the GNU General Public License as published by
8  * the Free Software Foundation; either version 3 of the License, or
9  * (at your option) any later version.
10  *
11  * This program is distributed in the hope that it will be useful,
12  * but WITHOUT ANY WARRANTY; without even the implied warranty of
13  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14  * GNU General Public License for more details.
15  *
16  * You should have received a copy of the GNU General Public License
17  * along with MIA; if not, see <http://www.gnu.org/licenses/>.
18  *
19  */
20 
21 #ifndef mia_2d_splinetransformpenalty_hh
22 #define mia_2d_splinetransformpenalty_hh
23 
24 #include <mia/2d/interpolator.hh>
25 #include <mia/2d/datafield.hh>
26 #include <mia/core/vector.hh>
27 #include <mia/core/product_base.hh>
28 
29 
31 
40 
41 public:
43  static const char *data_descr;
45  static const char *type_descr;
46 
51  C2DSplineTransformPenalty(double weight);
52 
54  C2DSplineTransformPenalty& operator = (const C2DSplineTransformPenalty& org) = delete;
55 
60  C2DSplineTransformPenalty *clone() const __attribute__((warn_unused_result));
61 
62  virtual ~C2DSplineTransformPenalty();
63 
70  void initialize(const C2DBounds& size, const C2DFVector& range, PSplineKernel kernel);
71 
77  double value(const C2DFVectorfield& coefficients) const;
78 
87  double value_and_gradient(const C2DFVectorfield& coefficients, CDoubleVector& gradient) const;
88 
89 protected:
90 
91 
92 
93  const C2DBounds& get_size() const;
94 
95  const C2DFVector& get_range() const;
96 
97  PSplineKernel get_kernel() const;
98 
99  double get_weight() const;
100 
101 private:
102 
103  virtual void do_initialize() = 0;
104 
105  virtual double do_value(const C2DFVectorfield& coefficients) const = 0;
106 
107  virtual double do_value_and_gradient(const C2DFVectorfield& coefficients, CDoubleVector& gradient) const = 0;
108 
109  virtual C2DSplineTransformPenalty *do_clone() const = 0;
110 
111  double m_weight;
112 
113  C2DBounds m_size;
114  C2DFVector m_range;
115  PSplineKernel m_kernel;
116 };
117 
118 typedef std::shared_ptr<C2DSplineTransformPenalty> P2DSplineTransformPenalty;
119 
120 
121 class EXPORT_2D C2DSplineTransformPenaltyPlugin: public TFactory<C2DSplineTransformPenalty> {
122 public:
123  C2DSplineTransformPenaltyPlugin(char const * const name);
124 private:
125  virtual Product *do_create() const __attribute__((warn_unused_result));
126  virtual Product *do_create(float weight) const __attribute__((warn_unused_result)) = 0 ;
127 
128  float m_weight;
129 };
130 
131 
135 
136 
138 
139 
140 EXPORT_2D C2DSplineTransformPenaltyPluginHandler::ProductPtr produce_2d_spline_transform_penalty(const std::string& descr);
141 
142 
144 public:
146 };
147 
149 
150 #endif