Support:Documents:Examples:Estimate Input Parameters Simultaneously with Flow

Estimating Input Function Parameters Simultaneously with Flow

This example demonstrates an approach to simultaneously estimating input functions parameters along with blood flow in the context of the PET blood flow model.

Tissue Uptake

The flow model has one tissue compartment and (radioactive) water in the blood is assumed to rapidly exchange with the water in the tissue. The tissue model has two rate constants: K₁ and k₂ and can be depicted in a diagram:

It can also be described by the differential equation:

dC_T/dt = K₁ C_p - k₂ C_T

where C_T is the tissue concentration and C_p is the plasma concentration of radioactive water.

Instead of estimating K₁ and k₂, it sometimes is advantageous estimate a different set of parameters defined as

lambda = K₁/k₂ and F = K₁.

F is blood flow (perfusion) per unit volume of tissue (we're assuming 100% extraction) and lambda is the ratio of tissue to blood concentration at equilibrium.

To use this parameterization in COMKAT, one needs to express the rate constants in terms of lambda and F: K₁ = F, k₂ = K₁/lambda

This example demonstrates how to estimate lambda and F.

Input Function

This example also demonstrates how to use a parameterized input function. Specifically, the plasma concentration vs. time is modeled as

C_p = (p₁t - p₂ - p₃) exp(p₄ t) + p₂ exp(p₅ t) + p₃ exp(p₆t) [aka Feng input]

This input function was proposed in Feng Int. J. Biomed. Comput. 32(2):95-110, 1993 and has been implemented in fin.m in the COMKAT examples folder.

Pixel Values

In a PET image, pixel values represent the radioactivity concentration. There are contributions from radioactivity in blood and tissue. The arterial blood concentration includes radioactive water in plasma and the blood cells and is designated C_a. Here we assume C_a can be described using the Feng input equation but with different values for p₁, p₂, ... p₆ than were used for C_p. The tissue concentration is C_T and is given above. The pixel activity concentration is calculated as the concentration in the tissue plus the pixel's fractional blood space F_v times the arterial blood concentration: C_T + F_vC_a.

COMKAT Implementation

T

Text

Write some introduction to your example.

Code

You can put code in a different format using <pre> tags like:

For text that you want to appear as part of a MATLAB function or program, encase it inside a "pre" block.
Line 2 of pre block
    Line 3 (indented) of pre block.

Edit the current page to see how this is done.

Figures

Remember to include figures and plots in your example to make it look pretty and convey information. To insert a figure, create a figure file in png (preferred) or other graphic format.

figure
plot(1:10)
title('Example Figure')
ylabel('y-axis label')
xlabel('x-axis label')
set(1,'PaperPosition',[0.25 2.5 3 2])
print -dpng ExampleFigure.png

Next, upload the image file ("upload file" at the bottom of any wiki page) giving it a unique name.

Use double-square brackets to insert the figure into the page

Edit the example index page to include your example in the list

Add your new example to the list of examples on the example index page by adding a line like this to the index page:

[http://comkat.case.edu/comkat/comkat_wiki/index.php?title=Support:Documents:Examples:Instructions_for_Adding_Your_Own_Example The title of you example]

Within the brackets, the first part is the URL address to your example page. Then you type your example title following that address. On the index page it will look like:

The title of you example

That is it. It is so simple that even a caveman (or professor) can do it.