imgbfk2 - tpcclib 0.8.0 © 2023 by Turku PET Centre

Estimation of rate constants K1, k2, and Vb from dynamic PET image in
ECAT 6.3, ECAT 7.x, NIfTI-1, or Analyze 7.5 file format using
one-tissue compartment model, solved using the basis function approach.
The radioactivity concentration of each voxel is assumed to consist of
tissue (Ct) and blood (Cb):
  Cpet(t) = Ct(t) + Vb*Cb
and differential equation for Ct is:
  dCt(t)/dt = K1*Cp(t) - k2*Ct(t)
, where Cp(t) is arterial plasma concentration curve.
Arterial plasma and blood data must be corrected for decay and delay.
Dynamic PET image must be corrected for decay.
Fit time must be given in seconds.
 
Usage: imgbfk2 [Options] ptacfile btacfile imgfile fittime K1file [k2file]
 
Options:
 -Vt=<filename>
     Parametric K1/k2 (Vt) image is saved.
 -Vb=<filename>
     Parametric Vb image is saved.
 -k2min=<Min k2> and -k2max=<Max k2>
     Enter the minimum and maximum k2 in units 1/min, applying to decay
     corrected data. Defaults are k2min=0 and k2max=10, respectively.
 -K1min=<Min K1> and -K1max=<Max K1>
     Enter the minimum and maximum K1 value; defaults are
     0 and 10 mL/(mL*min), respectively.
 -Vtmin=<Min Vt> and -Vtmax=<Max Vt>
     Enter the minimum and maximum value for apparent Vt (K1/k2);
     defaults are 0 and 10 mL/mL, respectively.
 -Vbmin=<Min Vb> and -Vbmax=<Max Vb>
     Enter the minimum and maximum value for Vb (volume fraction);
     defaults are 0 and 1 mL/mL, respectively.
 -nr=<value>
     Set number of basis functions; default is 500, minimum 100.
 -wss=<filename>
     Weighted sum-of-squares are written in specified image file.
 -thr=<threshold%>
     Pixels with AUC less than (threshold/100 x input AUC) are set to zero;
     default is 0%
 -bf=<filename>
     Basis function curves are written in specified file.
 -err=<filename>
     Excluded voxels and voxels with any of the result parameters at either
     of its limits is written in the specified image file with value 1,
     otherwise with value 0.
 -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.
 
Example 1. Calculation of K1 and k2 images using 0-20 min of
the dynamic image:
  imgbfk2 P223344apc.bld P223344ab.bld P223344.nii 20 P223344k1.nii P223344k2.nii
 
By default, the units of voxel values in the K1 image are
(mL plasma)/((mL tissue) * min), in k2 image 1/min, in Vb image
mL blood/mL tissue), and in Vt image (mL plasma)/(mL tissue),
 
See also: imgbfh2o, fitk2, imgunit, fitdelay, imgcbv
 
Keywords: image, modelling, compartmental model, basis function method