/usr/include/mia-2.2/mia/2d/timestep.hh is in libmia-2.2-dev 2.2.7-3.
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 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 | /* -*- mia-c++ -*-
*
* This file is part of MIA - a toolbox for medical image analysis
* Copyright (c) Leipzig, Madrid 1999-2015 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_timestep_hh
#define mia_2d_timestep_hh
#include <mia/core/factory.hh>
#include <mia/core/typedescr.hh>
#include <mia/2d/image.hh>
#include <mia/2d/transform.hh>
NS_MIA_BEGIN
/**
\ingroup registration
\brief The time step class for time-marching registration algorithms
Registration time step for a time-marching registration algorithm like
those based on a variational model.
Generally the time marhcing model consists of two parts, a PDE describing
the regularization model of the registration and the time step that
advances the registration transformation towards a minimum.
*/
class EXPORT_2D C2DRegTimeStep : public CProductBase {
public:
/// plugin search path data component helper
typedef C2DImage plugin_data;
/// plugin search path plugin type component helper
typedef timestep_type plugin_type;
/**
Constructor
\param min minimum time step
\param max maximum time step
*/
C2DRegTimeStep(float min, float max);
virtual ~C2DRegTimeStep();
/**
Evaluate the pertuberation of the vector field in combination with the next
transformation to be applied. What actually happens here depends on the
time step model.
\param[in,out] io vector field resulting from the solution of the PDE,
may be overwritted by its pertuberated version
\param shift current transformation
\returns the norm of the maxium transformation over the transformation domain
*/
float calculate_pertuberation(C2DFVectorfield& io,
const C2DTransformation& shift) const;
/**
Depending on the time step model, a regridding may be used - specifically this
is the case for the fluid dynamic model
\param b current transformation
\param v velocity field in case of a fluid dynamic model
\param delta time step
\returns true if regridding is necessary, false if not
*/
bool regrid_requested(const C2DTransformation& b, const C2DFVectorfield& v,
float delta) const;
/**
Decrease the time step by dividing by 2.0. - if the time step falls
below the appointed minimum it will be adjusted accordingly
\returns true if the time-step was decreased, and false if the time
step was already at the minimum
*/
bool decrease();
/// increase thetime step by multiplying with 1.5
void increase();
/** evaluate the time step based on the maximum shift resulting from
\a calculate_pertuberation
\param maxshift maximum shift allowed for all pixles
\returns the time step delta to be used
*/
float get_delta(float maxshift) const;
/** \returns true if the time step model supports regridding and false if not */
bool has_regrid () const;
private:
virtual float do_calculate_pertuberation(C2DFVectorfield& io,
const C2DTransformation& shift) const = 0;
virtual bool do_regrid_requested (const C2DTransformation& b,
const C2DFVectorfield& v, float delta) const = 0;
virtual bool do_has_regrid () const = 0;
float m_min;
float m_max;
float m_current;
};
/// pointer type for the 2D registration time step
typedef std::shared_ptr<C2DRegTimeStep > P2DRegTimeStep;
/**
\ingroup infrastructur
\brief Plug-in to create the time step evaluation
Plug-in to create the time step evaluation in time-marching registration
algorithms.
*/
class EXPORT_2D C2DRegTimeStepPlugin : public TFactory<C2DRegTimeStep>
{
public:
/**
Initialize the plug in with its given name
@param name
*/
C2DRegTimeStepPlugin(const char *name);
protected:
/// @returns minium time step
float get_min_timestep() const;
/// @returns maximum time setp
float get_max_timestep() const;
private:
float m_min;
float m_max;
};
/**
\ingroup infrastructur
Time step model plugin handler
*/
typedef THandlerSingleton<TFactoryPluginHandler<C2DRegTimeStepPlugin> >
C2DRegTimeStepPluginHandler;
NS_MIA_END
#endif
|