convsurg - tpcclib 0.8.0 © 2023 by Turku PET Centre
Convolution of PET time-activity curve (TAC) with response function h(t)
consisting of surge function
h(t) = a*t*exp(-b*t)
Surge function integral from 0 to infinity is a/(b*b).
Output TAC, Co(t), is calculated from input TAC, Ci(t), as
Co(T) = Ci(T) (x) h(T)
, where (x) denotes the operation of convolution.
Usage: convsurg [Options] tacfile a b outputfile
Sample time interval in convolution; by default the shortest interval
in the plasma data; too long interval as compared to b's leads to bias.
Scale response function to have integral from 0 to infinity to unity;
that will lead to similar AUC for the input and output TACs.
Display usage information on standard output and exit.
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.
Suppress displaying normal results on standard output.
Suppress displaying anything except errors.
The units of a and b must be compatible with units of the TAC, and
the optional sample interval.
For accurate results, input TAC should have very short sampling intervals.
Frame durations are not used, even if available in the TAC file.
convsurg -auc=1 plasma.tac 1 0.05 simulated.tac
See also: fit_xsur, convexpf, sim_av, fit2dat, tacadd, taccalc, interpol
Keywords: TAC, simulation, modelling, convolution