The most common implementation of autoradiography (ARG) in PET is the quantitation of perfusion (blood flow) from [15O]H2O bolus studies.
Autoradiography is also used in ex vivo quantitative phosphor imaging.
Autoradiography in the estimation of tissue perfusion with [15O]H2O
ARG method has been used to quantitate perfusion in the brain, tumours and skeletal muscle. The method has been validated for skeletal muscle in Turku (Ruotsalainen et al., 1997). ARG method is very robust and fast to calculate, and therefore it has been used mainly to compute parametric images of perfusion.
A specific compartmental model for radiowater has been developed, based on the principle of exchange of inert gas between blood and tissues. In the ARG method the apparent partition coefficient (p=K1/k2) is assumed to be known. Contribution of vascular radioactivity is assumed to be negligible, which leads to overestimation of the perfusion estimates:
- Bolus [15O]H2O infusion, or [15O]CO2 inhalation
- Arterial blood sampling, preferably with on-line sampler
- Static imaging starting when PET count rate starts to increase, or dynamic scanning; scan length must be relatively short when high perfusion values are to be measured (for example 90 s for the brain) and longer for low perfusion values (for example 360 s for resting skeletal muscle)
- Both PET and blood data are corrected for decay
- Blood TAC is corrected for dispersion and delay time (Iida et al., 1986)
Implementation of ARG method
1. Tissue TACs are simulated
2. Lookup table is constructed
These simulated tissue curves are integrated over a specified time, equal to the duration of static PET scan or integration time of dynamic PET scan.
The perfusion values used to simulate the TACs, and the corresponding TAC integral (AUC) values are stored in a “look-up table”.
Example of look-up table file:
0.00000e+000 0.00000e+000 3.24911e+001 2.40481e-001 6.48763e+001 4.80962e-001 9.71560e+001 7.21443e-001 1.29331e+002 9.61924e-001 1.61401e+002 1.20240e+000 1.93366e+002 1.44289e+000 2.25228e+002 1.68337e+000 2.56985e+002 1.92385e+000 2.88640e+002 2.16433e+000 3.20192e+002 2.40481e+000 3.51641e+002 2.64529e+000 ...
3. Integral image is computed
In case of dynamic PET scan, the AUC of each image voxel is computed during the integration time. In case of static PET scan, the image data is simply multiplied by the frame length.
4. Computation of perfusion image
For each voxel in the integral image, the nearest AUC value in the look-up table is found, and the corresponding blood flow value from the look-up table is written into the parametric image.
Computation of blood flow images in TPC
If you are still using Windows XP, and your computer is logged into TPC/hospital network, you can use the GUI or VBS script as instructed below. GUI does not work in Windows 7 or later. If your computer is outside TPC network or operating system is not Windows XP, you can do the same steps that the script would do in the command-line interface, preferably writing your own study protocol dependent script.
The perfusion image can then be used in SPM analysis, or regional flow values can be retrieved directly from the ROIs.
Calculation step-by-step with examples
Blood input TAC must be prepared first. In the next example the count-rate curve was not available, and it is made from the dynamic image, and then given to the input processing script. Dispersion in the tubing is corrected by the script automatically, and the physiological dispersion is assumed negligible in the example:
imghead -thr=30,fh P103192.v P103192head.tac water_input -fit:yes -cr:P103192head.tac -disp:0.0 P103192.lis P103192blof.bld tacunit -yconv=Bq/ml P103192blof.bld
Make sure that calibration units are the same in the PET image and blood data. In the example, PET image was calibrated to units Bq/ml, and the units of blood data were converted accordingly.
1. Simulate tissue TACs and compute the look-up tableTo compute the look-up table, the following command-line arguments need to be specified to the program arlkup:
- corrected arterial blood datafile (for Windows script see water_input.bat)
- partition coefficient of water (recommended values are 0.8 for the brain, 1.0 for tumours, 0.99 for muscle, and 0.19 for white adipose tissue)
- maximal blood flow expected to be found in the image, in units (mL blood)×(dL tissue)-1×min-1
- integration start time (sec); usually the time where blood curve starts to rise
- integration duration (sec); the dynamic image or sinogram must be integrated over the same time; in case of static scan, this is the duration of the scan
- filename for the look-up table.
If static PET scan (single frame) was performed, then you should use option
arlkup. If dynamic PET image, corrected for physical decay,
was summed into one frame, this option must not be used, but if dynamic PET image or sinogram, not
corrected for physical decay, was summed, then this option should be used.
Optionally, the size of the look-up table can be given; by default it is set to 5000 but, if maximal blood flow was not set to too high value, 2000 gives an eligible flow map, and is faster to compute.
Perfusion studies with count-based scan start: PET scanning may be started automatically when the count rate starts to rise. After delay correction, blood curve starts at the same time. Thus, the integration time for look-up table and for image (next step) is easy to set in these studies: start time is 0 and duration is the length of the scan in static imaging or the total length of required number of frames.
Continuing the previous example, here the look-up table is calculated from the corrected blood curve assuming that p=0.99 (skeletal muscle), maximum perfusion is 80 (mL blood)×(dL tissue)-1×min-1, PET scan was started when counts started to increase, and integration duration is set to 300 seconds:
arlkup P103192blof.bld 0.99 80 0 300 P103192.lkup
2. Integration of the image or sinogram
If the study was dynamic, the quality of the flow image will be better, if frames of the sinogram are summed before image reconstruction. To sum the selected frames of an ECAT sinogram use ecatssum.
After the dynamic or static image is available, an integral image must be computed using imginteg (if the image is dynamic, the frames should be selected so that their total length is the same as the integration time for the blood curve. The frame number and times can be checked with eframe).
Continuing the previous example, the dynamic image is integrated from 0 to 300 seconds:
imginteg P103192.v 0 300 P103192auc.v
3. Conversion of the integral image to the perfusion image
Program imglkup is used to convert the integral image to the perfusion image using the look-up table.
Continuing the previous example, the integrated image is converted to perfusion image:
imglkup P103192auc.v P103192.lkup P103192flow.v
If your Windows XP computer is logged into TPC/hospital network, you can use the pre-made script with either graphical user interface (GUI) or command-line user interface (CLI):
If you are using Windows XP and old MS Explorer, you can open the GUI by double-clicking file
S:\bin\windows\imgarg.hta, or, click Windows XP Start button, then Run, write
imgarg.hta and click Ok. GUI does not work in Windows 7 or later.
In the GUI, fill the required fields, and click “Execute” to compute the perfusion image.
This GUI uses script
imgarg.vbs to do the actual analysis;
this script can also be used from command-line (see below).
Open Windows command prompt window, type
This script has replaced the old UNIX scripts
blood data only upto ... sec.”
Blood sampling system introduces considerable delay to the blood before it reaches the online detectors. To be on the safe side, blood sampling should be continued about 60 s after the end of PET scanning. However, sometimes blood sampler has been stopped too soon, leading to this error.
Blood curve can then be extrapolated using extrapol. The resulting extrapolated curve must be checked visually, because extrapolation is based on fitting which may fail, or there may be some another reason for this error.
It is possible to apply ARG method to regional TACs, although this may be useful only in simulation studies to study the biases introduced by the ARG method.
For regional calculation, follow the step-by-step instructions above, except that in step 2 integrate the regional TAC data using program dftinteg, and in step 3 convert the TAC integrals to perfusion values using program taclkup.
Ginsberg MD, Howard BE, Hassel WR. Emission tomographic measurement of local cerebral blood flow in humans by an in vivo autoradiographic strategy. Ann Neurol. 1984; 15: S12-S18.
Herscovitch P, Markham J, Raichle ME. Brain blood flow measured with intravenous H215O. I. Theory and error analysis. J Nucl Med. 1983; 24: 782-789.
Herscovitch P, Raichle ME. What is the correct value for the brain-blood partition coefficient for water? J Cereb Blood Flow Metab. 1985; 5: 65-69.
Howard BE, Ginsberg MD, Hassel WR, Lockwood AH, Freed P. On the uniqueness of cerebral blood flow measured by the in vivo autoradiographic strategy and positron emission tomography. J Cereb Blood Flow Metab. 1983; 3: 432-441.
Iida H, Kanno I, Miura S, Murakami M, Takahashi K, Uemura K. Error analysis of a quantitative cerebral blood flow measurement using H215O autoradiography and positron emission tomography, with respect to the dispersion of the input function. J Cereb Blood Flow Metab. 1986; 6: 536-545.
Liukko K. Effects of decay correction in autoradiography method. TPCMOD0037.
Raichle ME. Quantitative in vivo autoradiography with positron emission tomography. Brain Res Rev. 1979; 1: 47-68.
Raichle ME, Martin WRW, Herscovitch P, Mintun MA, Markham J. Brain blood flow measured with intravenous H215O. II. Implementation and validation. J Nucl Med. 1983; 24: 790-798.
Ruotsalainen U, Raitakari M, Nuutila P, Oikonen V, Sipilä H, Teräs M, Knuuti J, Bloomfield PM, Iida H. Quantitative blood flow measurement of skeletal muscle using oxygen-15-water and PET. J Nucl Med. 1997; 38:314-319.
Created at: 2008-05-12
Updated at: 2017-09-05
Written by: Vesa Oikonen