p2t_v3c 1.2.1 © 2001-2013 by Turku PET Centre
Simulation of PET tissue time-radioactivity concentration curve (TAC)
from arterial plasma (Ca) and blood (Cb) TACs, based on three-tissue
compartmental model, where the compartments are in series (default):
____ K1 ____ k3 ____ k5 ____
| Ca | ----> | C1 | ----> | C2 | ----> | C3 |
|____| <---- |____| <---- |____| <---- |____|
k2 k4 k6
dC1(t)/dt = K1*Ca(T) - (k2+k3)*C1(T) + k4*C2(T)
dC2(t)/dt = k3*C1(T) - (k4+k5)*C2(T) + k6*C3(T)
dC3(t)/dt = k5*C2(T) - k6*C3(T)
Ct(T) = C1(T) + C2(T) + C3(T)
Cvb(T) = Cab(T) - dCt(t)/dt / f
Cpet(T)= Vb*fA*Cab(T) + Vb*(1-fA)*Cvb(T) + (1-Vb)*Ct(T)
, or, optionally, three-tissue compartmental model, where the 2nd and 3rd
tissue compartments are parallel, often used to represent specific and
non-specific binding:
____ K1 ____ k3 ____
| Ca | ----> | C1 | ----> | C2 |
|____| <---- |____| <---- |____|
k2 k4
| ^
k5 | | k6
v |
____
| C3 |
|____|
dC1(t)/dt = K1*Ca(T) - (k2+k3+k5)*C1(T) + k4*C2(T) + k6*C3(T)
dC2(t)/dt = k3*C1(T) - k4*C2(T)
dC3(t)/dt = k5*C1(T) - k6*C3(T)
Ct(T) = C1(T) + C2(T) + C3(T)
Cvb(T) = Cab(T) - dCt(t)/dt / f
Cpet(T)= Vb*fA*Cab(T) + Vb*(1-fA)*Cvb(T) + (1-Vb)*Ct(T)
Command-line arguments:
1) Plasma TAC filename
2) Blood TAC filename
3-8) K1, k2, k3, k4, k5, k6
9) Blood flow f (ml/(min*ml)); set to 0 to assume that f>>K1 and Cvb=Cab
10) Vascular volume fraction Vb (%)
11) Arterial fraction of vascular volume fA (%)
12) Filename for simulated TACs
Options:
-paral[lel]
Model with parallel compartments C2 and C3 is applied.
-ser[ies]
Model with compartments C1, C2, and C3 in series is applied (default).
-sub | -nosub
TACs of sub-compartments (C1, C2 and C3) are written (-sub, default)
or not written (-nosub) into the output file.
-h or --help
Print this message and exit.
-v, --version or --build
Print software build information and exit.
--silent
Program works silently, printing only error and warning messages.
--verbose
Program prints more information about what it is doing.
If the times in plasma file are in seconds, the units of rate constants
(k's) and blood flow (f) must also be specified as 1/sec.
For accurate results, plasma TAC should have very short sampling intervals.
To reduce the model, k5 or k3 can be set to 0.
Simulated TACs are written in ASCII format with columns:
1) Sample time
2) Total tissue activity concentration (Cpet)
3) Activity concentration in 1st tissue compartment, (1-Vb)*C1
4) Activity concentration in 2nd tissue compartment, (1-Vb)*C2
5) Activity concentration in 3rd tissue compartment, (1-Vb)*C3
6) Arterial concentration to tissue activity, Vb*fA*Cab
7) Venous contribution to tissue activity, Vb*(1-fA)*Cvb
References:
1. TPCMOD0001.
See also:
p2t_3c,
p2t_di,
dftadd,
dftren,
fr4sim,
dft2csv,
dft2img
Keywords: DFT, simulation, modelling, compartmental model
This program comes with ABSOLUTELY NO WARRANTY. This is free software, and
you are welcome to redistribute it under GNU General Public License.