22 December 2010

Data-driven arterial input functions

I haven't devoted any time/effort/code to the extraction and parameter estimation of arterial input functions (AIFs) for DCE-MRI in either the vignette or the submission to JSS.  Frankly, I wanted to communicate the main features of the R package dcemriS4 and data-driven AIFs just haven't made the list... up 'til now.  So let's assume one wants to fit a parametric model to observed contrast-agent concentration time curves (CTCs) obtained from a DCE-MRI acquisition.

The function extract.aif() requires a seed voxel to perform a three-dimensional region-growing algorithm based on correlation.  I use FSLView, but feel free to select the voxel using your favorite piece of software.  The code below uses the RIDER Neuro MRI data from the NBIA (National Biomedical Imaging Archive) that has been (a) converted from DICOM to NIfTI, (b) brain extracted by applying BET on the 15-degree multiple flip-angle acquisition and (c) converted from signal intensity to contrast agent concentration in all brain voxels.  The main steps in estimating a parametric AIF are
  1. Read in the contrast-agent CTCs (a four-dimensional NIfTI file)
  2. Read in the acquisition times (obtained from the DICOM data)
  3. Use extract.aif() with the seed voxel and correlation threshold
  4. Fit a parametric model (Orton et al. 2008) to the average CTC from the seed-growing algorithm



The thin black lines correspond to individual CTCs from the voxels, the red line is the average CTC used to fit the parametric model and the green curve is the fitted model. Note, in order to fit the entire time series a compromise was struck between fitting the first pass of the contrast versus the tail. It may be necessary to modify the range of acquisition times in order to emphasize the beginning or the end of the AIF, depending upon which aspect of the dynamic acquisition is more important.

Contrast-agent concentration time curves from voxel (128,75,9).

The four parameters that characterize the AIF model may be obtained via

2 comments:

  1. Nice blog and great article. This site is very helpful and informative. I am come back to see your website and I am waiting for your next post. Thanks for sharing with us.

    ReplyDelete
  2. I always enjoy visiting your site because I know that you’ll be posting something interesting every time and this article certainly did not disappoint. It’s a great read – thank you for putting in all that work into writing such an informative piece on my favorite topic!Womens Shearling Jackets

    ReplyDelete