imglhk1 - tpcclib 0.8.0 © 2023 by Turku PET Centre

Computation of parametric images of K1 and Vb, and optionally k2, from
dynamic PET image in ECAT, NIfTI, or Analyze format applying
(ir)reversible on1-tissue compartmental model with arterial plasma input.
The compartmental models are transformed to general linear least squares
functions (1), which are solved using Lawson-Hanson non-negative least
squares (NNLS) algorithm (2).
Dynamic PET image and plasma time-activity curve (PTAC) must be corrected
for decay to the tracer injection time.
Usage: imglhk1 [Options] ptacfile imgfile k1file vbfile [k2file]
     Pixels with AUC less than (threshold/100 x PTAC AUC) are set to zero;
     default is 1%.
 -end=<Fit end time (min)>
     Use data from 0 to end time; by default, model is fitted to all frames.
 -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 unidirectional back-transport rate k2 is considered in the model setting
only if file name for k2 image is given.
Note that this model can correctly estimate Vb only if
1) plasma does not contain any labelled metabolites, and
2) plasma and blood curves are similar in shape.
Note also that Cpet is modelled as Cpet=Vb*Cb + Ct, thus K1 may need to
be corrected by factor 1/(1-Vb).
The units of pixel values in the parametric images are ml/(min*ml) for K1
1/min for k2, and ml/ml for Vb.
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, ISBN 0-89871-356-0.
See also: imglhk3, imgcbv, imgki, fitk3
Keywords: image, modelling, irreversible uptake, NNLS