Slicer Extension providing pharmacokinetic modeling
C++ Python CMake
Fetching latest commit…
Cannot retrieve the latest commit at this time.
Failed to load latest commit information.


PkModeling is a 3D Slicer Version 4 Extension that provides pharmacokinetic modeling for dynamic contrast enhanced MRI (DCE MRI)[1][2].

PkModeling accepts volumetric timecourse data of signal intensities and computes parametric maps using either a two or three parameter Tofts model. The estimated parameters include

  • Ktrans - volume transfer contrast constant between plasma and the extracellular-extravascular space at each voxel
  • Ve - fractional volume for extracellular space at each voxel
  • fpv - fractional plasma volume at each voxel
  • MaxSlope - maximum slope of the of the time sequence at each voxel
  • AUC - area under the curve of each voxel, measured from the bolus arrival time to the end of the time interval, normalized by the AUC of the arterial input function (AIF)
  • R^2 - goodness of fit value. Since the parametric model is non-linear, R^2 is not strictly bounded by [-1,1]. But larger values still correspond to better fits.

PkModeling can also output a concentration curve view of the original volumetric timecourse as well as the "fitted" concentration curves resulting from the parametric model.

Estimation of the parametric model is controlled through a series of inputs including

  • T1 Blood Value
  • T1 Tissue Value - defaults to published value for prostate in de Bazelaire et al.[3]
  • Relaxivity Value - default 0.0039 corresponds to the Gd-DPTA (Magnevist) at 3T, see Pintaske et al.[4]. This value needs to be adjusted for magnet strength and contrast agent.
  • Hematocrit Value - volume percentage of red blood cells in blood
  • AUC Time Interval Value - time interval for AUC calculation

Furthermore, an arterial input function (AIF) must be specified either by designating a mask corresponding to the voxels on which to base a patient specific estimate of the AIF, or by specifying a population derived AIF curve directly.

Finally, the estimation of the parametric maps can be restricted to a specified mask defining a region of interest.

Acquisition parameters relevent to the parametric model fitting are embedded in the input volumetric timecourse data, either as attributes on a NRRD file or extracted directly from the underlying DICOM structures

  • TR Value - repetition time (in milliseconds)
  • TE Value - echo time (in milliseconds)
  • FA Value - flip angle (in degrees)
  • Timestamps for the timecourses (in milliseconds)


See the MultiVolumeExplorer module in the 3D Slicer.


@millerjv, @fedorov, @zhuy


[1]: Knopp MV, Giesel FL, Marcos H et al. "Dynamic contrast-enhanced magnetic resonance imaging in oncology." Top Magn Reson Imaging, 2001; 12:301-308.

[2]: Rijpkema M, Kaanders JHAM, Joosten FBM et al. "Method for quantitative mapping of dynamic MRI contrast agent uptake in human tumors." J Magn Reson Imaging 2001; 14:457-463.

[3]: de Bazelaire, C.M., et al. "MR imaging relaxation times of abdominal and pelvic tissues measured in vivo at 3.0 T: preliminary results." Radiology, 2004. 230(3): p. 652-9.

[4]: Pintaske J, Martirosian P, Graf H, Erb G, Lodemann K-P, Claussen CD, Schick F. "Relaxivity of Gadopentetate Dimeglumine (Magnevist), Gadobutrol (Gadovist), and Gadobenate Dimeglumine (MultiHance) in human blood plasma at 0.2, 1.5, and 3 Tesla." Investigative radiology. 2006 March;41(3):213–21.