imgflowd - tpcclib 0.8.0 © 2023 by Turku PET Centre
Estimation of rate constants K1 (perfusion), k2 and Va from dynamic
radiowater PET image in ECAT 6.3, ECAT 7.x, NIfTI-1, or Analyze 7.5 file
format using. The two-compartment model equations
dCt(t)/dt = K1*Ca(t) - k2*Ct(t)
Croi(t) = Va*Ca(t) + Ct(t)
are linearized (1) into equation
Croi(t) = Va*Ca(t) + (K1+Va*k2)*Integral[Ca(t)] - k2*Integral[Croi(t)]
from which the model parameters are solved using Lawson-Hanson non-negative
least squares (NNLS) method (2).
Arterial blood curve (BTAC) and dynamic PET image must be corrected for physical decay.
BTAC sample times must be in seconds.
Since time delay varies inside PET image, this program requires a delay map
which contains the time delay for each image pixel.
Usage: imgflowd [Options] btacfile imgfile delayfile flowfile
Options:
-end=<Fit end time (sec)>
Use data from 0 to end time; by default, all of it.
-k2=<filename>
Parametric k2 image is saved; in some situations perfusion calculation
from k2 can be more accurate than the default assumption of f=K1.
Perfusion can be calculated from k2 using equation f=k2*pH2O, where
pH2O is the physiological partition coefficient of water in tissue.
-Va=<filename> | -Va=0
Parametric Va image is saved, or set Va to 0 if image is pre-corrected
for arterial blood volume; by default Va is fitted.
-max=<Max value>
Upper limit for blood flow (K1) values.
-noneg
Pixels with negative K1 estimates are set to zero.
-h, --help
Display usage information on standard output and exit.
-v, --version
Display version and compile information on standard output and exit.
-d[n], --debug[=n], --verbose[=n]
Set the level (n) of debugging messages and listings.
-q, --quiet
Suppress displaying normal results on standard output.
-s, --silent
Suppress displaying anything except errors.
The units of pixel values in the blood flow (K1) image is
(mL blood)/((mL tissue) * min), in k2 image 1/min, and
in Va image (mL blood/mL tissue).
References:
1. Blomqvist G. On the construction of functional maps in positron
emission tomography. J Cereb Blood Flow Metab. 1984;4:629-632.
2. Lawson CL & Hanson RJ. Solving least squares problems.
Prentice-Hall, 1974.
See also: imgdelay, imgflowm, imgbfh2o, imgcbv
Keywords: image, modelling, perfusion, blood flow, radiowater, NNLS