Helper Class for NumericTwoTissueCompartment Model: Defines the differential equations (Mass Balance Equations) in the 2-tissue-compartment model for dynamic PET data modeling. The 2-Tissue Compartment model is defined via the mass balance equations dC1(t)/dt = K1*Ca(t) - (k2 + k3)*C1(t) + k4*C2(t) dC2(t)/dt = k3*C1(t) - k4*C2(t) CT(t) = C_a(t)*VB + (1-VB)*(C1(t)+C2(t) where Ca(t) is the plasma concentration(aterial input function) Boost ODEINT performs a stepwise numeric integration (e.g. via Runge-Kutta method) of the initial value problem x' = dx/dt = f(x,t) It needs an operator () (a functor) that calculates dx/dt = dxdt for a given x and t. Parameters are K1,k2,k3,k4, VB and the time dependent Ca(t) =AIF, that is interpolated to the current step t. More...

#include <mitkTwoTissueCompartmentModelDifferentialEquations.h>

Public Types
typedef std::vector< double >	AIFType

Public Member Functions
void	operator() (const mitk::NumericTwoTissueCompartmentModel::state_type &x, mitk::NumericTwoTissueCompartmentModel::state_type &dxdt, const double t)
	Functor for differential equation of Two Tissue Compartment Model Takes current state x = x(t) and time t and calculates the corresponding dxdt = dx/dt. More...

	TwoTissueCompartmentModelDifferentialEquations ()

void	initialize (double k_1, double k_2, double k_3, double k_4)
	Initialize class with parameters K1, k2, k3 and k4 that are free fit parameters. More...

void	setAIF (AIFType &aif)

void	setAIFTimeGrid (AIFType &grid)

Detailed Description

Helper Class for NumericTwoTissueCompartment Model: Defines the differential equations (Mass Balance Equations) in the 2-tissue-compartment model for dynamic PET data modeling. The 2-Tissue Compartment model is defined via the mass balance equations dC1(t)/dt = K1*Ca(t) - (k2 + k3)*C1(t) + k4*C2(t) dC2(t)/dt = k3*C1(t) - k4*C2(t) CT(t) = C_a(t)*VB + (1-VB)*(C1(t)+C2(t) where Ca(t) is the plasma concentration(aterial input function) Boost ODEINT performs a stepwise numeric integration (e.g. via Runge-Kutta method) of the initial value problem x' = dx/dt = f(x,t) It needs an operator () (a functor) that calculates dx/dt = dxdt for a given x and t. Parameters are K1,k2,k3,k4, VB and the time dependent Ca(t) =AIF, that is interpolated to the current step t.

Definition at line 32 of file mitkTwoTissueCompartmentModelDifferentialEquations.h.

Member Typedef Documentation

◆ AIFType

typedef std::vector< double > mitk::TwoTissueCompartmentModelDifferentialEquations::AIFType

Definition at line 36 of file mitkTwoTissueCompartmentModelDifferentialEquations.h.

Constructor & Destructor Documentation

◆ TwoTissueCompartmentModelDifferentialEquations()

mitk::TwoTissueCompartmentModelDifferentialEquations::TwoTissueCompartmentModelDifferentialEquations ( )

inline

Definition at line 49 of file mitkTwoTissueCompartmentModelDifferentialEquations.h.

Member Function Documentation

◆ initialize()

void mitk::TwoTissueCompartmentModelDifferentialEquations::initialize	(	double	k_1,
		double	k_2,
		double	k_3,
		double	k_4
	)

inline

Initialize class with parameters K1, k2, k3 and k4 that are free fit parameters.

Definition at line 54 of file mitkTwoTissueCompartmentModelDifferentialEquations.h.

◆ operator()()

void mitk::TwoTissueCompartmentModelDifferentialEquations::operator()	(	const mitk::NumericTwoTissueCompartmentModel::state_type &	x,
		mitk::NumericTwoTissueCompartmentModel::state_type &	dxdt,
		const double	t
	)

inline

Functor for differential equation of Two Tissue Compartment Model Takes current state x = x(t) and time t and calculates the corresponding dxdt = dx/dt.

Definition at line 41 of file mitkTwoTissueCompartmentModelDifferentialEquations.h.