imgdv - tpcclib 0.8.0 © 2023 by Turku PET Centre

Calculates the distribution volume (Vt, DV) or distribution volume ratio (DVR)
from dynamic PET image in ECAT, NIfTI-1, or Analyze 7.5 format,
applying the multiple time graphical analysis (MTGA) approach for reversible
tracer uptake (Logan plot) (1, 2, 3).
Simple non-iterative perpendicular line fitting algorithm (4) is applied
in determination of the slope of the Logan plot.
 
Usage: imgdv [Options] input image starttime resultimage
 
Input can be either a TAC file containing arterial plasma PTAC or reference
region TTAC. Data must be in the same concentration units as the data in
image, and times in min; however, if units are specified inside the input
and image file, the program can automatically convert units.
Start time of the linear fit must be given in minutes; set to zero to
let the program automatically determine the start time (not recommended).
 
Options:
 -k2=<k2 of reference region>
     With reference region input it may be necessary to specify also the
     population average for reference region k2 (2).
 -thr=<threshold%>
     Pixels with AUC less than (threshold/100 x input) are set to zero.
     Default is 0 %.
 -max=<Max value>
     Upper limit for Vt or DVR values; by default max is set pixel-wise
     to 10 times the AUC ratio.
 -min=<Min value>
     Lower limit for Vt or DVR values; 0 by default.
 -filter
     Remove parametric pixel values that are over 4x higher than
     their closest neighbours.
 -end=<Fit end time (min)>
     By default line is fitted to the end of data. Use this option to enter
     the fit end time.
 -v=<filename>
     Y axis intercepts times -1 are written as an image in specified file.
 -n=<filename>
     Numbers of selected plot data points are written as an image.
 -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:
     imgdv ua3818ap.kbq ua3818dy1.v 20 ua3818dv.v
 
The unit of pixel values in the parametric Vt image is
(mL plasma)/(mL tissue). Data in DVR image are unitless.
 
In theory, DVR equals BPnd+1, and therefore a binding potential image can
be computed from DVR image by subtraction of unity, e.g. using program
imgcalc. However, this may lead to a large number of pixels with negative
values that may hamper the further statistical or visual analysis.
 
References:
1. Logan J, Fowler JS, Volkow ND, Wolf AP, Dewey SL, Schlyer DJ,
   MacGregor RR, Hitzemann R, Bendriem B, Gatley SJ, Christman DR.
   Graphical analysis of reversible radioligand binding from time-activity
   measurements applied to [N-11C-methyl]-(-)-cocaine PET studies in human
   subjects. J Cereb Blood Flow Metab 1990; 10: 740-747.
2. Logan J, Fowler JS, Volkow ND, Wang GJ, Ding YS, Alexoff DL.
   Distribution volume ratios without blood sampling from graphical
   analysis of PET data. J Cereb Blood Flow Metab. 1996; 16: 834-840.
3. Logan J. Graphical analysis of PET data applied to reversible and
   irreversible tracers. Nucl Med Biol 2000; 27:661-670.
4. Varga J & Szabo Z. Modified regression model for the Logan plot.
   J Cereb Blood Flow Metab 2002; 22:240-244.
 
See also: imgbfbp, imglhdv, img2tif, logan, imgki, imgcalc
 
Keywords: image, MTGA, Logan plot, modelling, Vt, DVR, distribution volume