Skip to content

Section M: Perfusion Models

General forward model

Code OSIPI name Alternative names Notation Description Reference
M.GF1.001 Forward model -- -- A forward model to be inverted (select from Section M). --

MR signal models

This section covers models that describe how the measured MR signal S depends on electromagnetic properties, such as the relaxation rates R1, R1 and R2* or the magnetic susceptibility \(\chi\) , and on MR sequence parameters such as TR and TE.
The exception is section "Magnitude models: DCE - R1 in the fast water exchange limit, direct relationship with indicator concentration", in which the model describes how the MR signal depends directly on the indicator concentration.

Magnitude models: DSC

Code OSIPI name Alternative names Notation Description Reference
M.SM1.001 Gradient echo model -- GE model This forward model is given by the following equation:
\(S=S_0\cdot e^{-TE\cdot R_2^*}\)
with
\(TE\) (Q.MS1.005),
\(S_0\) (Q.MS1.010),
\(R_2^*\) (Q.EL1.007),
\(S\) (Q.MS1.001).
Jackson et al. 2005
M.SM1.002 Spin echo model -- SE model This forward model is given by the following equation:
\(S=S_0\cdot e^{-TE\cdot R_2}\)
with
TE (Q.MS1.005),
S0 (Q.MS1.010),
R2 (Q.EL1.004),
S (Q.MS1.001).
Jackson et al. 2005
M.SM1.999 Model not listed -- -- This is a custom free-text item, which can be used if a model of interest is not listed. Please state a literature reference and request the item to be added to the lexicon for future usage. --

Magnitude models: DCE - R1 in the fast water exchange limit

Code OSIPI name Alternative names Notation Description Reference
M.SM2.001 Linear model -- Linear This forward model is given by the following equation:
\(S=k \cdot R_1\)
with
k (Q.GE1.009),
R1 (Q.EL1.001),
S (Q.MS1.001).
--
M.SM2.002 Spoiled gradient recalled echo model FLASH model SPGR model This forward model is given by the following equation: \(S=S_0 \cdot \frac{sin(\alpha)[1-e^{TR\cdot R_1}]}{1-cos(\alpha)\cdot e^{-TR\cdot R_1}}\)
with
S0 (Q.MS1.010),
R1 (Q.EL1.001),
TR (Q.MS1.006),
\(\alpha\) (Q.MS1.007),
S (Q.MS1.001)
--
M.SM2.003 Single-shot saturation recovery SPGR with centric encoding model SS-SR-FLASH-c model SS-SR-SPGR-c model This forward model is given by the following equation: \(S = S_0\cdot(1-e^{-PD\cdot R_1})\)
with
S0 (Q.MS1.010),
R1 (Q.EL1.001),
PD (Q.MS1.008),
S (Q.MS1.001)
Parker et al. 2000
M.SM2.004 Saturation-recovery SPGR with linear encoding model SR-turboFLASH-lin model SR-turboSPGR-lin model This forward model is given by the following equation: \(S=S_0\cdot sin(\alpha)\cdot[(1-e^{-PD\cdot R_1})a^{n-1}+b\frac{(1-a^{n-1})}{(1-a)}]\)
with
\(a=cos(\alpha)\cdot e^{-TR\cdot R_1}\),
\(b=1-e^{-TR\cdot R_1}\),
S0 (Q.MS1.010),
R1 (Q.EL1.001),
PD (Q.MS1.008),
TR (Q.MS1.006),
\(\alpha\) (Q.MS1.006)= 90°,
n (Q.MS1.011),
S (Q.MS1.001)
Larson 2001
M.SM2.005 Single-shot inversion recovery SPGR with centric encoding model SS-IR-FLASH-c model SS-IR-SPGR-c model This forward model is given by the following equation: \(S = S_0\cdot(1-2e^{-PD\cdot R_1})\)
with
S0 (Q.MS1.010),
R1 (Q.EL1.001),
PD (Q.MS1.008),
S (Q.MS1.001)
Ordidge et al. 1990
M.SM2.006 Inversion-recovery spoiled gradient recalled echo with linear encoding model IR-turboFLASH-lin model IR-turboSPGR-lin model This forward model is given by the following equation: \(S=S_0\cdot sin(\alpha)\cdot [\frac{(C+bA-\frac{1}{cos(\alpha)}D+1)}{1+BD}\cdot e^{-PD\cdot R_1}a^{n-1}\)
  \(+(1-e^{-PD\cdot R_1})a^{n-1}+b\frac{1-a^{n-1}}{1-a}]\)
with
\(a=cos(\alpha)e^{-TR\cdot R_1}\),
\(b=1-e^{-TR\cdot R_1}\),
\(C=a^{N-1}(1-e^{-PD\cdot R_1})\),
\(A=\frac{1-a^{N-1}}{1-a}\),
\(D=cos(\alpha)e^{-PD\cdot R_1}\),
\(B=a^{N-1}e^{-PD\cdot R_1}\),
S0 (Q.MS1.010),
R1 (Q.EL1.001),
PD (Q.MS1.008),
\(\alpha\) (Q.MS1.007) ,
TR (Q.MS1.006),
n (Q.MS1.011),
N (Q.MS1.012),
S (Q.MS1.001)
Larson 2001
M.SM2.999 Model not listed -- -- This is a custom free-text item, which can be used if a model of interest is not listed. Please state a literature reference and request the item to be added to the lexicon for future usage. --

Magnitude models: DCE - R1 in the fast water exchange limit, direct relationship with indicator concentration

In this section models are described which assume a direct relationship between the measured MR signal S and the indicator concentration.

Code OSIPI name Alternative names Notation Description Reference
M.MS3.001 Linear model -- Linear This forward model is given by the following equation:
\(S=k\cdot C\)
with
C (Q.IC1.001),
k (Q.GE1.009),
S (Q.MS1.001)
--
M.MS3.002 Absolute signal enhancement model -- ASE This forward model is given by the following equation:
\(\left\| S-S_{BL} \right\|=k\cdot C\) ,
with
C (Q.IC1.001),
SBL (Q.MS1.002),
k (Q.GE1.009),
S (Q.MS1.001)
Ingrisch and Sourbron 2013
M.MS3.003 Relative signal enhancement model -- RSE This forward model is given by the following equation:
\(\left\| \frac{S}{S_{BL}}-1 \right\|=k\cdot C\)
with
C (Q.IC1.001),
SBL (Q.MS1.002),
k (Q.GE1.009),
S (Q.MS1.001)
Ingrisch and Sourbron 2013
M.SM3.999 Model not listed -- -- This is a custom free-text item, which can be used if a model of interest is not listed. Please state a literature reference and request the item to be added to the lexicon for future usage. --

Magnitude models: DCE - R1 in the presence of restricted water exchange

Notation: \(S_{DCE,FXL}\) is one of the DCE signal models from the table above (DCE - R1 in fast water exchange limit). In all the models below, water exchange across red blood cell membranes is considered to be in the fast exchange limit. We also assume that the tissue-blood partition coefficient for water is equal to 1. Models with restricted water-exchange between one or more compartments are also referred to as shutter-speed models.

Magnitude models: DCE - R1 in the presence of no water exchange

Code OSIPI name Alternative names Notation Description Reference
M.SM4.001 Fast transendothelial water exchange, no transcytolemmal water exchange -- -- This forward model is given by the following equation:
\(S = \left( p_{b} + p_{e}\right)S_{DCE,FXL}\left( R_{1,be}\right)+p_iS_{DCE,FXL}\left( R_{1,i}\right)\),
where \(S_{DCE,FXL}\) is a forward model selected from Magnitude models: DCE- R1 in fast water exchange limit,
\(p_b\) (Q.PH1.017.b),
\(p_e\) (Q.PH1.017.e),
\(p_i\) (Q.PH1.017.i),
\(R_{1,be}\) (Q.EL1.001.be),
\(R_{1,i}\) (Q.EL1.001.i),
S (Q.MS1.001)
--
M.SM4.002 No transendothelial water exchange, fast transcytolemmal water exchange -- -- This forward model is given by the following equation:
\(S = p_bS_{DCE,FXL}\left( R_{1,b}\right)+\left( p_e + p_i\right)S_{DCE,FXL}\left( R_{1,ei}\right)\),
where \(S_{DCE,FXL}\) is a forward model selected from Magnitude models: DCE- R1 in fast water exchange limit,
\(p_b\) (Q.PH1.017.b),
\(p_e\) (Q.PH1.017.e),
\(p_i\) (Q.PH1.017.i),
\(R_{1,b}\) (Q.EL1.001.b),
\(R_{1,ei}\) (Q.EL1.001.ei),
S (Q.MS1.001)
--
M.SM4.003 No water exchange -- -- This forward model is given by the following equation:
\(S = p_bS_{DCE,FXL}\left( R_{1,b}\right)+p_eS_{DCE,FXL}\left( R_{1,e}\right)+p_iS_{DCE,FXL}\left( R_{1,i}\right)\),
where \(S_{DCE,FXL}\) is a forward model selected from Magnitude models: DCE- R1 in fast water exchange limit,
\(p_b\) (Q.PH1.017.b),
\(p_e\) (Q.PH1.017.e),
\(p_i\) (Q.PH1.017.i),
\(R_{1,b}\) (Q.EL1.001.b),
\(R_{1,e}\) (Q.EL1.001.e),
\(R_{1,i}\) (Q.EL1.001.i),
S (Q.MS1.001)
Bains et al. 2010
M.SM4.999 Model not listed -- -- This is a custom free-text item, which can be used if a model of interest is not listed. Please state the doi of a literature reference and request the item to be added to the lexicon for future usage. --

Magnitude models: DCE - R1 in the presence of finite water exchange (Shutter speed models)

Non-zero exchange for the restricted ones

Code OSIPI name Alternative names Notation Description Reference
M.SM5.001 Fast transendothelial water exchange, finite transcytolemmal water exchange -- -- This forward model is given by the following equation:
\(S=p_+S_{DCE,FXL}\left( R_{1,+}\right)+\left( 1-p_+\right)S_{DCE,FXL}\left( R_{1,-}\right)\),
with
\(R_{1\pm}=0.5\left( R_{1,i}+\tau_i^{-1}+R_{1,be}+\tau_{be}^{-1} \right)\)
   \(\pm0.5\sqrt{\left( R_{1,i}+\tau_i^{-1}-R_{1,be}-\tau_{be}^{-1} \right)^2+4\tau_i^{-1}\tau_{be}^{-1}}\)

\(p_+=0.5-0.5\frac{\left( R_{i,i}-R_{1,be}\right)\left( 2(p_{be}+p_i)-1\right)+\tau_i^{-1}+\tau_{be}^{-1}}{\sqrt{\left( R_{1,i}+\tau_i^{-1}-R_{1,be}-\tau_{be}^{-1} \right)^2+4\tau_i^{-1}\tau_{be}^{-1}}}\),
where \(S_{DCE,FXL}\) is a forward model selected from Magnitude models: DCE- R1 in fast water exchange limit,
\(p_{be}\) (Q.PH1.017.be),
\(p_i\) (Q.PH1.017.i),
\(R_{1,be}\) (Q.EL1.001.be),
\(R_{1,i}\) (Q.EL1.001.i),
\(\tau_{be}\) (Q.PH1.016.be),
\(\tau_i\) (Q.PH1.016.i),
S (Q.MS1.001)
Buckley 2018
M.SM5.002 Finite transendothelial exchange, fast transcytolemmal water exchange -- -- This forward model is given by the following equation:
\(S=p_+S_{DCE,FXL}\left( R_{1,+}\right)+\left( 1-p_+\right)S_{DCE,FXL}\left( R_{1,-}\right)\),
with
\(R_{1\pm}=0.5\left( R_{1,ei}+\tau_{ei}^{-1}+R_{1,b}+\tau_{b}^{-1} \right)\)
   \(\pm0.5\sqrt{\left( R_{1,ei}+\tau_{ei}^{-1}-R_{1,b}-\tau_b^{-1} \right)^2+4\tau_{ei}^{-1}\tau_b^{-1}}\)

\(p_+=0.5-0.5\frac{\left( R_{1,ei}-R_{1,b}\right)\left( 2(p_{ei}+p_b)-1\right)+\tau_{ei}^{-1}+\tau_b^{-1}}{\sqrt{\left( R_{1,ei}+\tau_{ei}^{-1}-R_{1,b}-\tau_b^{-1} \right)^2+4\tau_{ei}^{-1}\tau_b^{-1}}}\),
where \(S_{DCE,FXL}\) is a forward model selected from Magnitude models: DCE- R1 in fast water exchange limit,
\(p_{b}\) (Q.PH1.017.b),
\(p_{ei}\) (Q.PH1.017.ei),
\(R_{1,b}\) (Q.EL1.001.b),
\(R_{1,ei}\) (Q.EL1.001.ei),
\(\tau_{b}\) (Q.PH1.016.b),
\(\tau_{ei}\) (Q.PH1.016.ei),
S (Q.MS1.001)
Schwarzbauer et al. 1997
Dickie 2019
M.SM5.003 Finite transendothelial exchange, finite transcytolemmal water exchange -- -- This forward model is given by the following equation:
\(S=S_b+S_e+S_i\),
where
\(\begin{pmatrix} S_b\\S_e\\S_i\end{pmatrix}=S_{DCE,FXL}(\textbf{X})\begin{pmatrix} p_b\\p_e\\p_i\end{pmatrix}\)
where
\(\textbf{X}=\begin{pmatrix} -R_{1,b}-\tau_b^{-1}&\tau_e^{-1}-\frac{p_i\tau_i^{-1}}{p_e}&0\\\tau_b^{-1}&-R_{1,e}-\tau_e^{-1}&\tau_i^{-1}\\0&\tau_e^{-1}-\frac{p_b\tau_e^{-1}}{p_e}&-R_{1,i}-\tau_i^{-1}\end{pmatrix}\)
where \(S_{DCE,FXL}\) is a forward model selected from Magnitude models: DCE- R1 in fast water exchange limit with 1’s replaced with the 3 x 3 identity matrix,
\(p_b\) (Q.PH1.017.b),
\(p_e\) (Q.PH1.017.e),
\(p_i\) (Q.PH1.017.i),
\(R_{1,b}\) (Q.EL1.001.b),
\(R_{1,e}\) (Q.EL1.001.e),
\(R_{1,i}\) (Q.EL1.001.i),
\(\tau_{b}\) (Q.PH1.016.b),
\(\tau_e\) (Q.PH1.016.e),
\(\tau_i\) (Q.PH1.016.i),
S (Q.MS1.001)
Ruiliang et al. 2020
Bains et al. 2010
M.SM5.999 Model not listed -- -- This is a custom free-text item, which can be used if a model of interest is not listed. Please state the doi of a literature reference and request the item to be added to the lexicon for future usage. --

Magnitude models: Combined DCE/DSC - R1/R2/R2*

Code OSIPI name Alternative names Notation Description Reference
M.SM6.001 DSC Multi-echo (GE) model -- -- This forward model is given by the following equation:
\(S=S_{DCE,FXL}(R_1)e^{-TE\cdot R_2^*}\)
where \(S_{DCE,FXL}\)(\(R_1\)) is a forward model selected from Magnitude models: DCE- R1 in fast water exchange limit,
TE (Q.MS1.005),
R2* (Q.EL1.007),
S (Q.MS1.001)
--
M.SM6.002 DSC Multi-echo (SE) model -- -- This forward model is given by the following equation:
\(S=S_{DCE,FXL}(R_1)e^{-TE\cdot R_2}\)
where \(S_{DCE,FXL}\)(\(R_1\)) is a forward model selected from Magnitude models: DCE- R1 in fast water exchange limit,
TE (Q.MS1.005),
R2 (Q.EL1.004),
S (Q.MS1.001)
--
M.SM6.999 Model not listed -- -- This is a custom free-text item, which can be used if a model of interest is not listed. Please state a literature reference and request the item to be added to the lexicon for future usage. --

Phase models: DSC

Code OSIPI name Alternative names Notation Description Reference
M.SM7.001 Linear susceptibility signal model -- -- This forward model is given by the following equation:
\(S=k\cdot \chi\),
with
k (Q.GE1.009),
\(\chi\) (Q.EL1.011),
S (Q.MS1.001).
--
M.SM7.999 Model not listed -- -- This is a custom free-text item, which can be used if a model of interest is not listed. Please state a literature reference and request the item to be added to the lexicon for future usage. --

Electromagnetic property models

This section covers models that describe how electromagnetic properties (EP), such as relaxation rates R1, R2 and R2* or the magnetic susceptibility \(\chi\) , are modulated by the indicator concentrations.

Code OSIPI name Alternative names Notation Description Reference
M.EL1.001 Transverse relaxation rate (GE), linear with relaxivity model Effective relaxation rate (GE), linear with relaxivity model -- This forward model is given by the following equation:
\(R_2^*=R_{20}^*+r_2^*\cdot C\)
with
R20* (Q.EL1.008),
r2* (Q.EL1.017),
C (Q.IC1.001),
R2* (Q.EL1.007)
(Rosen et al. 1990)
M.EL1.002 Transverse relaxation rate (SE), linear with relaxivity model Natural relaxation rate (GE), linear with relaxivity model -- This forward model is given by the following equation:
\(R_2=R_{20}+r_2\cdot C\)
with
R20 (Q.EL1.005),
r2 (Q.EL1.016),
C (Q.IC1.001),
R2 (Q.EL1.004)
(Rosen et al. 1990)
M.EL1.003 Longitudinal relaxation rate, linear with relaxivity model -- -- This forward model is given by the following equation:
\(R_1=R_{10}+r_1\cdot C\)
with
R10 (Q.EL1.002),
r1 (Q.EL1.015),
C (Q.IC1.001),
R1 (Q.EL1.001)
(Rosen et al. 1990)
M.EL1.004 Transverse relaxation rate (GE) with gradient leakage correction model -- -- This forward model is given by the following equation:
\(R_2^*=R_{20}^*+r_{2v}^*\left\| C_p-C_e \right\| +r_{2e}^*C_e,\)
with
\(R_{20}^*\) (Q.EL1.008),
\(C_p\) (Q.IC1.001.p),
\(C_e\) (Q.IC1.001.e),
\(r_{2e}^*\) (Q.EL1.017.e),
\(r_{2v}^*\) (Q.EL1.017.v),
\(R_2^*\)(Q.EL1.007)
Sourbron et al 2012
M.EL1.005 Transverse relaxation rate (GE), quadratic model -- -- This forward model is given by the following equation:
\(R_2^*=R_{20}^*+k_1\cdot C_p+k_2C_p^2\)
with
\(R_{20}^*\) (Q.EL1.008),
\(C_p\) (Q.IC1.001.p),
[\(k_1\),\(k_2\)] (Q.EL1.020),
\(R_2^*\) (Q.EL1.007)
Van Osch 2003 (also see Calamante 2013)
M.EL1.006 Linear susceptibility concentration model -- -- This forward model is given by the following equation:
\(\chi=\chi_0+\delta\chi\cdot C\)
with
\(\chi_0\) (Q.EL1.012),
\(\delta\chi\) (Q.EL1.013),
\(C\) (Q.IC1.001),
\(\chi\) (Q.EL1.011)
(Conturo et al. 1992)
M.EL1.999 Model not listed -- -- This is a custom free-text item, which can be used if a model of interest is not listed. Please state a literature reference and request the item to be added to the lexicon for future usage. --

Indicator concentration models

This section covers models that describe how indicator concentrations in tissue and blood vary with time.

Indicator kinetic models

In the current version of the lexicon the list of indicator kinetic models is restricted to linear and stationary tissues, and specific models with two distribution spaces. A summary of common pharmacokinetic models in contrast-agent based perfusion MRI is given in Sourbron and Buckley 2013. We provide the differential equations and impulse response functions using the consistent parameterizations to enable straightforward comparisons between models.

Code OSIPI name Alternative names Notation Description Reference
M.IC1.001 Linear and stationary system model -- LSS model This forward model is given by the following equations:
\(C(t)=I(t)\ast C_{a,p}(t)\)
with
[I (Q.IC1.005), t (Q.GE1.004)],
[\(C_{a,p}\) (Q.IC1.001.a,p), \(t\) (Q.GE1.004)],
[\(C_t\) (Q.IC1.001.t), \(t\) (Q.GE1.004)]
(Rempp et al. 1994)
M.IC1.002 One-compartment, no indicator exchange model -- 1CNEX model The one compartment no indicator exchange model describes an intravascular model with no vascular to EES indicator exchange. This forward model is given by the following differential equation:
\(v_{p}\frac{dC_{t}(t)}{dt} = F_{p}C_{a,p} - F_{p}C_{c,p}(t)\)
The impulse response function is given by:
\(I(t) = F_{p}e^{{-\frac{F_{p}}{v_{p}}t}}\)
with
[\(C_{a,p}\) (Q.IC1.001.a,p), t (Q.GE1.004)],
[\(C_{c,p}\) (Q.IC1.001.c,p), t (Q.GE1.004)],
[\(C_t\) (Q.IC1.001.t), t (Q.GE1.004)],
[I (Q.IC1.005), t (Q.GE1.004)],
\(F_p\) (Q.PH1.002),
\(v_{p}\) (Q.PH1.001.p)
(Tofts et al. 1999)
M.IC1.003 One-compartment, fast indicator exchange model -- 1CFEX model The one compartment fast exchange model describes infinitely fast bi-directional exchange of indicator between vascular to extravascular extracellular spaces. The capillary and EES effectively act as a single compartment. This forward model is given by the following differential equation:
\(\frac{dC_{t}(t)}{dt} = F_{p}C_{a,p} - \frac{F_{p}}{v_{p} + v_{e}}C_{t}(t)\)
The impulse response function is given by:
\(I(t) = F_{p}e^{{-\frac{F_{p}}{v_{p} + v_{e}}t}}\)
with
[\(C_{a,p}\) (Q.IC1.001.a,p), t (Q.GE1.004)],
[\(C_t\) (Q.IC1.001.t), t (Q.GE1.004)],
[I (Q.IC1.005), t (Q.GE1.004)],
\(F_p\) (Q.PH1.002),
\(v_{p}\) (Q.PH1.001.p),
\(v_{e}\) (Q.PH1.001.e)
(Sourbron et al. 2013)
M.IC1.004 Tofts Model Kety model, Generalized Kinetic Model TM The Tofts model describes bi-directional exchange of indicator between vascular to extravascular extracellular spaces. The capillary compartment is assumed to have negligible volume. The EES is modeled as well-mixed compartment. The forward model is given by the following differential equation:
\(\frac{dC_{t}(t)}{dt} = K^{trans}C_{a,p} - \frac{K^{trans}}{v_{e}}C_{t}(t)\)
The impulse response function is given by:
\(I(t) = K^{trans}e^{{-\frac{K^{trans}}{v_{e}}t}}\)
with
[\(C_{a,p}\) (Q.IC1.001.a,p), t (Q.GE1.004)],
[\(C_t\) (Q.IC1.001.t), t (Q.GE1.004)],
[I (Q.IC1.005), t (Q.GE1.004)],
\(K^{trans}\) (Q.PH1.008),
\(v_{e}\) (Q.PH1.001.e)
(Tofts and Kermode 1991)
M.IC1.005 Extended Tofts Model Modified Tofts Model, Extended Generalized Kinetic Model, Modified Kety model ETM The extended Tofts model describes bi-directional exchange of indicator between vascular to extravascular extracellular spaces. The capillary and EES are modeled as well-mixed compartments. It is equivalent to the 2CXM in the highly perfused limit. Dispersion of indicator within the capillary bed is assumed negligible:
\(C_{c,p} = C_{a,p}\)
The forward model is given by the following differential equation:
\(v_{e}\frac{dC_{e}(t)}{dt} = PSC_{c,p} - PSC_{e}(t)\)
The impulse response function is given by:
\(I(t) = v_{p}\delta(t) + K^{trans}e^{{-\frac{K^{trans}}{v_{e}}t}}\)
with
[\(C_{c,p}\) (Q.IC1.001.c,p), t (Q.GE1.004)],
[\(C_e\) (Q.IC1.001.e), t (Q.GE1.004)],
[I (Q.IC1.005), t (Q.GE1.004)],
\(\delta\) (M.DM1.009),
PS (Q.PH1.004),
\(v_{e}\) (Q.PH1.001.e),
\(K^{trans}\) (Q.PH1.008)
(Tofts 1997)
M.IC1.006 Patlak Model -- PM The Patlak model allows uni-directional exchange of indicator from vascular to extravascular extracellular spaces. Indicator exchange from the EES to the intravascular space is considered negligible during the timeframe of the imaging experiment. The capillary and EES are modeled as well-mixed compartments. It is equivalent to the two compartment uptake model in the highly perfused limit. Dispersion of indicator within the capillary bed is assumed negligible:
\(C_{c,p} = C_{a,p}\)
The forward model is given by the following differential equation:
\(v_{e}\frac{dC_{e}(t)}{dt} = PSC_{c,p}\)
The impulse response function is given by:
\(I(t) = v_{p}\delta(t) + PS\)
with
[\(C_{c,p}\) (Q.IC1.001.c,p), t (Q.GE1.004)],
[\(C_e\) (Q.IC1.001.e), t (Q.GE1.004)],
[I (Q.IC1.005), t (Q.GE1.004)],
\(\delta\) (Q.PH1.009),
PS (Q.PH1.004),
\(v_{p}\) (Q.PH1.001.p)
(Patlak et al. 1983)
M.IC1.007 Two compartment uptake model -- 2CUM The 2CU model allows uni-directional exchange of indicator from vascular to extravascular extracellular spaces. Indicator exchange from the EES to the intravascular space is considered negligible during the timeframe of the imaging experiment. The capillary and EES are modeled as well-mixed compartments. The forward model is given by the following differential equations:
\(v_{p}\frac{dC_{c,p}(t)}{dt} = F_{p}C_{a,p} - F_{p}C_{c,p} - PSC_{a,p}\)

\(v_{e}\frac{dC_{e}(t)}{dt} = PSC_{a,p}\)
The impulse response function is given by:
\(I(t) = F_{p}e^{-({\frac{F_{p} + PS}{v_{p}}})t} + E(1 - e^{-({\frac{F_{p} + PS}{v_{p}}})t})\)
with
[\(C_{a,p}\) (Q.IC1.001.a,p), t (Q.GE1.004)],
[\(C_{c,p}\) (Q.IC1.001.c,p), t (Q.GE1.004)],
[\(C_{e}\) (Q.IC1.001.e), t (Q.GE1.004)],
[I (Q.IC1.005), t (Q.GE1.004)],
\(F_p\) (Q.PH1.002),
PS (Q.PH1.004),
E (Q.PH1.005),
\(v_{e}\) (Q.PH1.001.e),
\(v_{p}\) (Q.PH1.001.p)
(Pradel et al. 2003), (Sourbron 2009)
M.IC1.008 Plug flow uptake model -- PFUM The plug flow uptake model allows uni-directional exchange of indicator from vascular to extravascular extracellular spaces. Indicator exchange from the EES to the intravascular space is considered negligible during the timeframe of the imaging experiment. The capillary space is modeled as a plug flow system and the EES as a well-mixed compartment. The forward model is given by the following differential equations:
\(v_{p}\frac{\partial C_{c,p}(x_{ax}, t)}{\partial t} = -F_{p}L_{ax}\frac{\partial C_{a,p}(x_{ax}, t)}{\partial x_{ax}} - PSC_{a,p}(x_{ax}, t)\)

\(v_{e}\frac{dC_{e}(t)}{dt} = PS \int_{0}^{L_{ax}} C_{c,p} (x_{ax},t) dx\)
The impulse response function is ... TO ADD IRF
with
[\(C_{c,p}\) (Q.IC1.001.c,p), t (Q.GE1.004)],
[\(C_{e}\) (Q.IC1.001.e), t (Q.GE1.004)],
[I (Q.IC1.005), t (Q.GE1.004)],
\(F_p\) (Q.PH1.002),
PS (Q.PH1.004),
\(v_{e}\) (Q.PH1.001.e),
\(v_{p}\) (Q.PH1.001.p),
\(L_{ax}\) (Q.GE1.007),
\(x_{ax}\) (Q.GE1.008)
(St. Lawrence and Frank 2000)
M.IC1.009 Two compartment exchange model -- 2CXM The 2CX model allows bi-directional exchange of indicator between vascular and extravascular extracellular compartments. Indicator is assumed to be well mixed within each compartment. The forward model is given by the following differential equations:
\(v_{p}\frac{dC_{c,p}}{dt} = F_{p}C_{a,p}(t) - F_{p}C_{c,p}(t) - PSC_{c,p}(t) + PSC_{e}(t)\)

\(v_{e}\frac{dC_{e}}{dt} = PSC_{c,p}(t) - PSC_{e}(t)\)
The impulse response function is given by
\(I(t) = F_{p}e^{-K_{+}t} + E_{-}(e^{-K_{+}t} - e^{-K_{-}t})\)

\(K_{\pm} = \frac{1}{2}\Biggl(\frac{F_{p} + PS}{v_{p}} + \frac{PS}{v_{p}} \pm \sqrt{(\frac{F_{p} + PS}{v_{p}} + \frac{PS}{v_{e}})^{2} - 4\frac{F_{p}PS}{v_{p}v_{e}}}\Biggl)\)

\(E_{-} = \frac{K_{+} + \frac{F_{p}}{v_{p}}}{K_{+} + K_{-}}\)
with
[\(C_{a,p}\) (Q.IC1.001.a,p), t (Q.GE1.004)],
[\(C_{c,p}\) (Q.IC1.001.c,p), t (Q.GE1.004)],
[\(C_{e}\) (Q.IC1.001.e), t (Q.GE1.004)],
[I (Q.IC1.005), t (Q.GE1.004)],
\(F_p\) (Q.PH1.002),
PS (Q.PH1.004),
\(v_{e}\) (Q.PH1.001.e),
\(v_{p}\) (Q.PH1.001.p)
(Brix et al. 2004), (Sourbron et al. 2009), (Donaldson et al. 2010)
M.IC1.010 Distributed parameter model -- DPM This is a placeholder for the distributed parameter model (Sangren and Sheppard 1953) (Sourbron 2011)
M.IC1.011 Tissue homogeneity model Johnson-Wilson model THM This is a placeholder for the tissue homogeneity model (Johnson and Wilson 1966) (Lawrence and Lee 1998) (Kershaw 2010) (Koh et al. 2003)
M.IC1.012 Adiabatic Approximation to the Tissue homogeneity model AATHM This is a placeholder for the adiabatic approximation to the tissue homogeneity model (Lawrence and Lee 1998) Kershaw et al. 2010) (Sourbron et al. 2012)
M.IC1.013 Two compartment filtration model -- 2CFM The 2CFM models unidirectional flow (filtration) from a vascular compartment into an extravascular compartment. A fraction of the filtrate (1-f) is reabsorbed. This model is appropriate for the kidney cortex or whole kidney. The forward model is given by the following differential equations:
\(v_{p}\frac{dC_{c,p}}{dt} = F_{p}C_{a,p}(t) - F_{p}C_{c,p}(t) - PSC_{c,p}(t)\)

\(v_{e}\frac{dC_{e}}{dt} = PSC_{c,p}(t) - (1 - f)PSC_{e}(t)\)
. The impulse response function is given by:
\(I(t) = v_{p}C_{c,p} + PSe^{-t\frac{(1-f)PS}{v_{e}}} \circledast C_{c,p}\)

where \(C_{c,p} = \frac{F_{p}}{v_{p}}e^{-t\frac{F_{p}}{v_{p}}}\)
with
[\(C_{a,p}\) (Q.IC1.001.a,p), t (Q.GE1.004)],
[\(C_{c,p}\) (Q.IC1.001.c,p), t (Q.GE1.004)],
[\(C_{e}\) (Q.IC1.001.e), t (Q.GE1.004)],
[I (Q.IC1.005), t (Q.GE1.004)],
\(F_p\) (Q.PH1.002),
PS (Q.PH1.004),
\(v_{e}\) (Q.PH1.001.e),
\(v_{p}\) (Q.PH1.001.p),
\(f\) (Q.PH1.018)
Sourbron et al. 2008
M.IC1.999 Model not listed -- -- This is a custom free-text item, which can be used if a model of interest is not listed. Please state a literature reference and request the item to be added to the lexicon for future usage. --

Arterial input function models

Code OSIPI name Alternative names Notation Description Reference
M.IC2.001 Parker AIF model -- Parker AIF This forward model is given by the following equation:
\(C_{a,b}(t)=\sum_{n=1}^{2}\frac{A_n}{\sigma_n\sqrt{2\pi}}e^{-\frac{(t-T_n)^2}{2\sigma_n^2}}+\frac{\alpha e^{-\beta t}}{1+e^{-s(t-\tau)}},\)
where \(A_n\), \(T_n\) and \(\sigma_n\) are the scaling constants, center and widths of the nth Gaussian; \(\alpha\) and \(\beta\) are the amplitude and decay constants of the exponential; and \(s\) and \(\tau\) are the width and center of the sigmoid, and [\(C_{a,b}\) (Q.IC1.001.a,b), t (Q.GE1.004)]. If not specified otherwise, the values from the publication are assumed: [\(A_1\), \(A_2\), \(T_1\), \(T_2\), \(\sigma_1\), \(\sigma_2\), \(\alpha\) , \(\beta\) , s, \(\tau\) ] = [48.54 mmol \(\cdot\) s, 19.8 mmol \(\cdot\) s, 10.2276 s, 21.9 s, 3.378 s, 7.92 s, 1.050 mmol, 0.0028 s-1, 0.6346 s-1, 28.98 s].
(Parker et al. 2006)
M.IC2.002 Georgiou AIF model -- Georgiou AIF The AIF between the start of the nth recirculation and (n+1)th recirculation is given by:
\(C_{a,p}(t)=(\sum_{i=1}^{3}A_ie^{-m_it})\cdot(\sum_{j=0}^{n}\gamma((j+1)\alpha+j,\beta,t-j\tau)),\)
with \(n\tau<t<(n+1)\tau\)
and \(\gamma(\alpha, \beta, \tau)=\frac{t^\alpha e^{-\frac{t}{\beta}}}{\beta^{\alpha+1}\Gamma(\alpha +1)}\) with \(t\geq 0\),
where \(A_1\), \(A_2\), \(A_3\) and \(m_1\), \(m_2\), \(m_3\) are the amplitudes and time constants of the exponential terms, the coefficients and reflect the number of theoretical mixing chambers and the ratio of the volume of the mixing chambers to the volume flow rate, and is the recirculation time, [\(C_{a,p}\) (Q.IC1.001.a,p), t (Q.GE1.004)]. \(\gamma = 0\) for \(t<0\).
If not specified otherwise, the values from the publication are assumed: [ \(A_1\), \(A_2\), \(A_3\), \(m_1\), \(m_2\), \(m_3\), \(\alpha\), \(\beta\) , \(\tau\) ] = [0.37 mM, 0.33 mM, 10.06 mM, 0.002 \(s^{-1}\), 0.02 \(s^{-1}\), 0.267 \(s^{-1}\), 5.26, 1.92 s, 7.74 s].
(Georgiou et al. 2019)
M.IC2.003 Weinmann AIF model -- Weinmann AIF This forward model is given by the following equation:
\(C_{a,p}(t)=D(a_1e^{-m_1t}+a_2e^{-m_2t}),\)
where \(D\) is the dose of contrast agent and \(a_1\), \(a_2\), \(m_1\) and \(m_2\) are the amplitudes and time constants of the exponential terms, and [\(C_{a,p}\) (Q.IC1.001.a,p), t (Q.GE1.004)].
If the model parameters are not specified, the values from the publication are assumed: [\(D\), \(a_1\), \(m_1\), \(a_2\), \(m_2\)] = [0.25 mmol/kg, 3.99 kg/l, 0.0024 \(s^{-1}\), 4.78 kg/l, 0.0002 \(s^{-1}\)].
(Weinmann et al. 1984)
M.IC2.999 Model not listed -- -- This is a custom free-text item, which can be used if a model of interest is not listed. Please state a literature reference and request the item to be added to the lexicon for future usage. --

Descriptive models

Code OSIPI name Alternative names Notation Description Reference
M.DM1.001 Unit step model Heaviside step model u This forward model is given by the following equation:
\(f(x-T)=0,\ \ x\leq T\),
\(f(x-T)=1,\ \ x\gt T\) ,
where T is a defined data grid point, at which the step function changes from 0 to 1 and \([x, f(x)]\)= [Data grid (Q.GE1.001), Data (Q.GE1.002)].
--
M.DM1.002 Linear-quadratic model -- LQM This forward model is given by the following equation:
\(f(x)=f_{BL},\ \ x\leq BAT\),
\(f(x)=f_{BL}+\beta_1(x-BAT)+\beta_2(x-BAT)^2,\)
\(x\gt BAT\),
where BAT is the bolus arrival time (Q.BA1.001), fBL the baseline (Q.BL1.001), \(\beta\)1 the slope after the BAT, \(\beta\)2 a quadratic component and \([x, f(x)]\)= [Data grid (Q.GE1.001), Data (Q.GE1.002)].
(Cheong et al. 2003)
M.DM1.003 Two step linear model -- 2SLM This forward model is given by the following equation:
\(f(x)=f_{BL},\ \ x\leq BAT,\)
\(f(x)=f_{BL}+b_1(x-BAT),\ \ x\gt BAT,\)
where BAT is the bolus arrival time (Q.BA1.001), fBL is the baseline (Q.BL1.001), b1 the slope after the BAT and \([x, f(x)]\)= [Data grid (Q.GE1.001), Data (Q.GE1.002)].
Cheong et al. 2003
M.DM1.004 Three step linear mode -- 3SLM This forward model is given by the following equation:
\(f(x)=f_{BL},\ \ x\leq BAT,\)
\(f(x)=f_{BL}+b_1(x-BAT),\ \ BAT\leq x\leq \beta,\)
\(f(x)=f_{BL}+b_1(x-BAT)+b_2(x-\beta),\ \ x\gt\beta\) ,
where BAT is the bolus arrival time (Q.BA1.001), fBL is the baseline (Q.BL1.001), \(\beta\) is the point of intersection of the 2nd and 3rd line segment. b1 and b2 are the slopes of the 2nd and 3rd line segments and \([x, f(x)]\)= [Data grid (Q.GE1.001), Data (Q.GE1.002)].
(Singh et al. 2009)
M.DM1.005 Multi-exponential model -- -- This forward model is given by the following equation:
\(f(x)=A_1\cdot e^{-x\cdot a_1}+;...+A_n\cdot e^{-x\cdot a_n}\)
where A1, …, An and a1,...,an are arbitrary coefficients and \([x, f(x)]\)= [Data grid (Q.GE1.001), Data (Q.GE1.002)].
--
M.DM1.006 Gamma-variate model -- -- This forward model is given by the following equation:
\(f(x)=\frac{1}{\Gamma(\alpha)\beta^{\alpha}}(x-BAT)^{\alpha-1}e^{-(x-BAT)/\beta}\)
where BAT is the bolus arrival time (Q.BA1.001), \(\alpha\) is a shape parameter, \(\beta\) is a scale parameter, \(\Gamma\)(\(\alpha\)) is the gamma function and \([x, f(x)]\)= [Data grid (Q.GE1.001), Data (Q.GE1.002)].
Mouridsen 2006
M.DM1.007 Fermi model -- -- This forward model is given by the following equation:
\(f(x)=F\cdot\frac{1+b}{1+b\cdot e^{x\cdot a}}\) ,
where \(F\) is the blood flow, \(a\) und \(b\) are arbitrary fit parameters and \([x, f(x)]\)= [Data grid (Q.GE1.001), Data (Q.GE1.002)].
(Brinch et al. 1999)
M.DM1.008 Normal distribution model Gaussian distribution model N This forward model is given by the following equation:
\(f(x)=\frac{1}{\sigma\sqrt{2\pi}}e^{-\frac{(x-\mu)^2}{2\sigma^2}},\)
where \(\mu\) is the mean (population) (Q.US1.007), \(\sigma\) is the standard deviation (Q.US1.010) and \([x, f(x)]\)= [Data grid (Q.GE1.001), Data (Q.GE1.002)].
--
M.DM1.009 Dirac delta model Unit pulse model \(\delta\) This forward model is given by the following equation:
\(\delta(x)=+ \infty\) for \(x = 0\),
\(\delta = 0\), elsewhere
with \([x, \delta(x)]\)= [Data grid (Q.GE1.001), Data (Q.GE1.002)].
--
M.DM1.999 Model not listed -- -- This is a custom free-text item, which can be used if a model of interest is not listed. Please state a literature reference and request the item to be added to the lexicon for future usage. --

Leakage correction models

This section is concerned with models for DSC leakage correction. They are not descriptive models in the sense that they are defined for very specific physical quantities, but at the same time cannot be derived as a composition of kinetic models, electromagnetic tissue property models or MR signal models.

Code OSIPI name Alternative names Notation Description Reference
M.LC1.001 Boxerman-Schmainda-Weisskoff (BSW) leakage correction model -- BSW leakage correction model This forward model is given by the following equation:
\(R_2^*(t)=R_{20}^*+K_1\overline{\Delta R_{2,ref}^*(t)}-K_2\int_0^t \overline{\Delta R_{2,ref}^*(t')dt'},\)
with
[\(\overline{\Delta R_{2,ref}^*}\)(Q.EL1.010), t (Q.GE1.004)]
\(R_{20}^*\) (Q.EL1.008),
\(K_1\) (Q.LC1.001),
\(K_2\) (Q.LC1.002),
[\(R_2^*\) (Q.EL1.007), t (Q.GE1.004)]
--
M.LC1.002 Bidirectional leakage correction model -- -- This forward model is given by the following equation:
\(R_2^*(t)=R_{20}^*+K_1\overline{\Delta R_{2,ref}^*(t)}\)
    \(-K_2\int_0^t \overline{\Delta R_{2,ref}^*(t')}\cdot e^{-k_{e->p}(t-t')}dt',\)
with
[\(\overline{\Delta R_{2,ref}^*}\)(Q.EL1.010), t (Q.GE1.004)]
\(R_{20}^*\) (Q.EL1.008),
\(K_1\) (Q.LC1.001),
\(K_2\) (Q.LC1.002),
\(k_{e->p}\) (Q.PH1.009.e->p)
[\(R_2^*\) (Q.EL1.007), t (Q.GE1.004)]
--
M.LC1.999 Model not listed -- -- This is a custom free-text item, which can be used if a model of interest is not listed. Please state a literature reference and request the item to be added to the lexicon for future usage. --

Perfusion identity models

This group lists relationships between perfusion quantities that can be used to derive one quantity from another under certain assumptions. This group is divided into the derivation of scalar quantities and the derivation of scalar derived from dynamic curves.

Scalar quantities

Code OSIPI name Alternative names Notation Description Reference
M.ID1.001 Central volume theorem -- CVT This forward model is given by the following equation:
\(v_p=MTT\cdot F_p\)
with
MTT (Q.PH1.006),
\(F_p\) (Q.PH1.002),
\(v_p\) (Q.PH1.001.p)
--
M.ID1.002 Total volume of distribution -- -- This forward model is given by the following equation:
\(v=v_p+v_e+v_i\)
with
\(v_p\) (Q.PH1.001.p),
\(v_e\) (Q.PH1.001.e),
\(v_i\) (Q.PH1.001.i),
\(v\) (Q.PH1.001)
--
M.ID1.003 Blood vs plasma volume fraction -- -- This forward model is given by the following equation:
\(v_b=\frac{v_p}{(1-Hct)}\)
with
\(v_p\) (Q.PH1.001.p),
\(Hct\) (Q.PH1.012),
\(v_b\) (Q.PH1.001.b).
--
M.ID1.004 Blood vs plasma flow -- -- This forward model is given by the following equation:
\(F_b=\frac{F_p}{(1-Hct)}\)
with
\(F_p\) (Q.PH1.002),
\(Hct\) (Q.PH1.012),
\(F_b\) (Q.PH1.003)
--
M.ID1.005 Blood vs plasma AIF -- -- This forward model is given by the following equation:
\(C_{a,b}(t)=C_{a,p}(t)\cdot(1-Hct_a)\),
with
[\(C_{a,p}\) (Q.IC1.001.a,p), t (Q.GE1.004)],
[\(C_{a,b}\) (Q.IC1.001.a,b), t (Q.GE1.004)],
\(Hct_a\) (Q.PH1.012.a)
--
M.ID1.006 Small vessel hematocrit correction -- -- This forward model is given by the following equation:
\(Hct_f=\frac{1-Hct_a}{1-Hct_t}\)
with
\(Hct_a\) (Q.PH1.012.a),
\(Hct_t\) (Q.PH1.012.t),
\(Hct_f\) (Q.PH1.013)
--
M.ID1.007 Compartment extraction fraction -- -- This forward model is given by the following equation:
\(E=\frac{PS}{F_p+PS}\)
with
\(PS\) (Q.PH1.004),
\(F_p\) (Q.PH1.002),
\(E\) (Q.PH1.005)
--
M.ID1.008 Plug flow extraction fraction -- -- This forward model is given by the following equation:
\(E=1-e^{-\frac{PS}{F_p}}\)
with
\(PS\) (Q.PH1.004),
\(F_p\) (Q.PH1.002),
\(E\) (Q.PH1.005)
--
M.ID1.009 Plasma MTT identity -- -- This forward model is given by the following equation:
\(MTT_p=\frac{v_p}{F_p+PS}\)
with
\(v_p\) (Q.PH1.001.p),
\(PS\) (Q.PH1.004),
\(F_p\) (Q.PH1.002),
\(MTT_p\) (Q.PH1.006.p)
--
M.ID1.010 Interstitial MTT identity -- -- This forward model is given by the following equation:
\(MTT_e=\frac{v_e}{PS}\)
with
\(v_e\) (Q.PH1.001.e),
\(PS\) (Q.PH1.004),
\(MTT_e\) (Q.PH1.006.e)
--
M.ID1.011 \(K^{trans}\) identity -- -- This forward model is given by the following equation:
\(K^{trans}=E\cdot F_p\),
with
E (Q.PH1.005),
\(F_p\) (Q.PH1.002),
\(K^{trans}\) (Q.PH1.008)
--
M.ID1.012 \(k_{ep}\) identity -- -- This forward model is given by the following equation:
\(k_{ep}=\frac{K^{trans}}{v_e}\),
\(K^{trans}\) (Q.PH1.008),
\(v_e\) (Q.PH1.001.e),
\(k_{e->p}\) (Q.PH1.009.e->p)
--
M.ID1.999 Model not listed -- -- This is a custom free-text item, which can be used if a model of interest is not listed. Please state a literature reference and request the item to be added to the lexicon for future usage. --

Scalars derived from dynamic curves

Code OSIPI name Alternative names Notation Description Reference
M.ID2.001 Bolus delay identity -- -- This forward model is given by the following equation:
\(MTT_a=\int_{0}^{\infty}h_a(t)dt\)
with
[ \(h_a\) (Q.IC1.004), t (Q.GE1.004)],
\(MTT_a\) (Q.PH1.006.a)
--
M.ID2.002 Tissue mean transit time identity -- -- This forward model is given by the following equation:
\(MTT_t=\int_{0}^{\infty}R(t)dt\)
with
[ \(R\) (Q.IC1.002),
\(t\) (Q.GE1.004)],
\(MTT_t\) (Q.PH1.006.t)
--
M.ID2.003 Blood plasma flow from maximum -- -- This forward model is given by the following equation:
\(F_p=max(I(t))\)
with
[ \(I\) (Q.IC1.005), \(t\) (Q.GE1.004)],
\(F_p\) (Q.PH1.002)
--
M.ID2.004 Blood plasma flow from first time frame -- -- This forward model is given by the following equation:
\(F_p=I(0)\)
with
[ \(I\) (Q.IC1.005), \(t\) (Q.GE1.004)],
\(F_p\) (Q.PH1.002)
--
M.ID2.005 Capillary transit time heterogeneity identity -- -- This forward model is given by the following equation:
\(CTTH=\int_{0}^{\infty}\sqrt{(t-MTT)^2}h(t)dt\)
with
[\(h\) (Q.IC1.003), \(t\) (Q.GE1.004)],
\(MTT\) (Q.PH1.006),
\(CTTH\) (Q.PH1.014)
--
M.ID2.999 Model not listed -- -- This is a custom free-text item, which can be used if a model of interest is not listed. Please state a literature reference and request the item to be added to the lexicon for future usage. --