imgmtrap - tpcclib 0.8.0 © 2023 by Turku PET Centre

Computation of parametric images of K1 and Vb from dynamic PET image in
ECAT, NIfTI, or Analyze format, applying NNLS.
The multilinear model is either
  Cpet(t) = Vb*Cb(t) + ((1-Vb)*K1+Vb*k2)*Integral[Cb(t)]
            - k2*Integral[Cpet(t)]
or if plasma data is available
  Cpet(t) = Vb*Cb(t) + (1-Vb)*K1*Integral[(1-HCT)*Cp(t)]
 
The program needs the blood curve and dynamic PET image, and optionally
plasma curve, multiplied by (1-HCT).
The name for the resulting parametric K1 image must be given.
 
Usage: imgmtrap [Options] btacfile imgfile k1file [vbfile]
 
Options:
 -k2=<value> | -k2=median
     Parameter k2 is fixed to a known value, or to a median k2 from
     initial fitting of unconstrained model.
 -ptac=<filename>
     Use plasma curve, multiplied by (1-HCT), as input function;
     the -k2 options are not available with this.
 -thr=<threshold%>
     Pixels with last frame concentration less than (threshold/100 x BTAC)
     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.
 -w1, -wf, -wfa
     By default, all weights are set to 1.0 (no weighting, option -w1);
     option -wf sets weights based on frame lengths, and option -wfa based
     on both frame lengths and mean activity during each frame.
 -p2=<filename>
     Image with pixel values of the second fitted parameter, (1-Vb)*K1.
 -p3=<filename>
     Image with pixel values of the third fitted parameter, HCT*kc;
     only available with the first model setting (with unknown PTAC).
 -cbv=<filename>
     Image with pixel values corrected for estimated blood volume,
     Cpet(t) - (1-Vb)*Cb(t).
 -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 parametric images are ml/(min*ml) for K1
and ml/ml for Vb.
 
See also: fitmtrap, b2ptrap, img2dft
 
Keywords: image, modelling, irreversible uptake, perfusion