Introduction

The two-compartment cerebral blood flow model can be characterized by the following three equations:

In these equations, is the PET data collected over a set of frames, is the delay and dispersion corrected arterial blood sample data, and is the radioactive tracer activity present in cerebral tissue. We know both and , but cannot know without knowing .

One additional point to consider is that these equations are written for continuous functions. However, the PET data that is actually available is average activity across each frame. Therefore, in order to solve the equations, we must average the non PET data across frames. By combining equations (2) and (3), and averaging the non PET data over frames, we get:

where is the PET data that is actually collected (by definition, averaged over frames), and performs an averaging over frames ( is the frame start time, is the frame stop time, and the integral is evaluated for each frame).

We wish to solve equation (4) for , , and . Of course, one approach would be to try to perform an explicit least squares curve fitting. However, this approach would be quite computationally intensive since the fitting would need to be performed for every pixel of a 128x128 pixel image. Fortunately, there is a method of solution that gives good results with reduced computational difficulty.