oxygen-15 positron emission tomography ( o pet ) can be used to measure cerebral blood flow ( cbf ) , cerebral blood volume ( cbv ) , cerebral metabolic rate of oxygen ( cmro2 ) , and oxygen extraction fraction ( oef ) [ 18 ] . several methods have been proposed for determining cbf , cbv , cmro2 , and oef using o pet . administration of radiotracers using a steady - state technique , which is recognized for its capacity to provide oef values with very high accuracy , is one of those methods [ 7 , 1012 ] . this method is based on the continuous administration of o - labeled tracers ( usually by inhalation ) until constant activity concentrations are reached in the tissue and blood ( dynamic equilibrium ) . quantitative measures of cbf , cmro2 , and oef can then be obtained from the observed activity concentrations using very simple equations . disadvantages of the steady - state technique are the need for stable arterial concentrations [ 10 , 14 ] , the relatively high radiation dose , and the need to assume a value for the volume of distribution of water . in addition , not all cyclotrons have the capacity to produce o - labeled tracers continuously . an alternative method to calculate cbf , cmro2 , and oef is the use of dynamic scanning protocols together with bolus injections and/or brief inhalations of the o - labeled tracers . generation of parametric cbf data using a basis function approach has been discussed previously [ 15 , 16 ] . in order to generate parametric cmro2 and oef images with a dynamic scanning protocol , several methods have been described [ 1720 ] . day - to - day variability of cbf , cmro2 , and oef measurements using dynamico pet scans in addition to absolute values , left - to - right ratios of cbf , cmro2 , and especially oef can be important for identifying patients with impaired cerebral perfusion caused by occlusive cerebral artery disease , who might be treated with extracranial to intracranial bypass surgery [ 21 , 22 ] . again , day - to - day variability of these ratios using dynamic o pet scans has not been determined . the purpose of the present study was to generate reference data of cbf , cmro2 , and oef using dynamic o pet scans and to evaluate associated day - to - day test  retest variability . test  retest pet studies were performed in ten healthy volunteers ( six men , four women ) with a mean age of 69  years ( range , 57 to 80  years ) who were scanned on two separate occasions ( with an interval of 354  days ) . in addition , all patients underwent magnetic resonance imaging ( mri ) , which was needed for pet image analysis shortly before the first pet study . the study was approved by the medical ethics committee of the university medical center utrecht and all volunteers provided informed consent . pet scans were acquired using an ecat exact hr+ scanner ( cti / siemens , knoxville , tn , usa ) which has been characterized elsewhere . subjects were scanned under baseline conditions , i.e. , lights in the room were dimmed and music was turned off . heads were positioned in an immobilization device to minimize patient motion during the study and in between scans . each study consisted of [ o]h2o , [ o]o2 , and [ o]co emission scans . prior to the emission scans , for the subsequent cbf study , a bolus ( ~5  s ) of 1,100  mbq [ o]h2o was administered intravenously , simultaneously starting a 3d dynamic emission scan ( 25 frames over a period of 600  s ) . the [ o]o2 study was performed 5  min after the end of the cbf study to allow for decay of residual o activity to background levels . a 3d dynamic emission scan ( 20 frames over a period of 600  s ) was then started simultaneously with the net inhalation of approximately 300  mbq [ o]o2 gas during a period of ~30  s. during administration of the [ o]o2 gas , the subjects continued breathing normally to avoid changes in cerebral hemodynamics , but they were instructed to inhale through the nasal administration cannula ( oxygen catheter dk3460 ; diameter , 4.7  mm ; unomedical , birkerd , denmark ) and exhale through the mouth . note that , in this study , scan durations of 600  s were used for both the water and oxygen studies . initially , it was attempted to use these longer scans to derive both perfusion and oxygen extraction fraction from a single [ o]o2 study . however , it was found that this was not feasible , and data were therefore analyzed by reusing perfusion and volume of distribution from the water scan ( see  image analysis  ) . although , by reusing parameters from the water scan , shorter scan durations may have been applied , it would not have shortened the overall study duration as some time interval between scans is needed to allow for radioactive decay . for the cbv study , a bolus of approximately 200  mbq [ o]co gas was administrated in a manner identical to the [ o]o2 procedure . one minute after the end of inhalation , a 3d dynamic emission scan ( three frames over a period of 360  s ) was started . all scans were reconstructed using fore  +  2d fbp reconstruction with a hanning filter at nyquist frequency , including all usual corrections for quantification . a matrix size of 256    256 and a zoom of 2.1 were applied , resulting in a voxel size of 1.2    1.2    2.4  mm and a final image resolution of ~7  mm full width at half maximum . for all dynamic scans , the arterial input function was measured using an online continuous blood - sampling device ; details on this device are given in . in short , the device consists of a bismuth germanate ( bgo)-based detector , shielded by a 3-cm lead , a withdrawal pump , and a computer for reading out the detector and controlling the pump . during blood withdrawal , a 1-mm - inner - diameter tubing that is connected to an arterial cannula is positioned in a slit that runs through the bgo detector . on the other end , the tubing is connected to a 3-mm - inner - diameter tubing that runs through the pump , and finally blood is collected in a 3-cm lead shielded waste bag . using scripts , both detector readout sampling and pump speed ( i.e. , withdrawal rates ) can set and automatically change during the studies . when necessary , pump speed can be overruled manually . during the first 4  min of each study , blood is withdrawn at a rate of 450  ml / h , and readout of the detector is performed continuously using 1-s readout samples . after 4  min onwards , pump speed is set to 150  ml / h to minimize the amount of blood withdrawn , and the detector readout sampling time is set to 10  s. manual samples were taken at discrete times for calibration purposes and for assessment of plasma - to - whole blood ratios : three samples at 5.5 , 8 , and 10  min during the [ o]h2o scan , six samples at 1 , 2 , 3 , 4.5 , 7 , and 9  min during the [ o]o2 scan , and three samples at 1 , 3 , and 5  min during the [ o]co scan . activity concentrations in the blood and plasma were determined for all samples . for the samples taken during the oxygen scan , the plasma - to - whole blood ratios were determined in order to separate the whole blood input function into water and oxygen input functions needed for kinetic analysis . for the latter purpose , a multi - exponential function was fitted through the six plasma - to - whole blood ratios . parametric cbf and the volume of distribution of water ( vt ) images were generated using a basis function implementation of the blood flow model , including corrections for dispersion , delay , and arterial blood volume as described previously . next , parametric cmro2 and oef images were generated using a slightly modified implementation of the method published by hattori et al .   with the following differences : ( 1 ) as scans were acquired over 600  s , all frames during this interval were included in fitting the data ; ( 2 ) the parametric vt data derived from the cbf scans were reused voxelwise rather than using the ( same ) fixed value for all voxels [ 16 , 25 ] ; and ( 3 ) the blood volume fraction was included as fit parameter rather than reusing it from the co study as it has been shown that this procedure provides more accurate fits near large blood vessels . more details about the implementation of the basis function method are given in the appendix . parametric cbv data were derived from the co scan as described by martin et al .   and hattori et al . . subsequently , an mri template with volumes of interest ( vois ) of the 3d flow territories of middle cerebral artery ( mca ) , anterior cerebral artery ( aca ) , and posterior cerebral artery ( pca ) ( fig . warping was performed using the automated image registration ( air ) software developed by woods et al . . nonrigid 12-parameter perspective warping was applied to fit the template onto the subject 's mri scan . statistical parametric mapping ( spm02 , wellcome department of cognitive neurology , london , uk , application written in matlab 5.2 ; mathworks , inc . , natick , md , usa ) was used for gray  white matter segmentation of the coregistered mri scan . these gray and white matter segmentations were then used in combination with the warped voi template to calculate mean cbf , cbv , cmro2 , and oef in the different flow territories , separating data / voi into gray and white matter subregions . the generation of cbf , cbv , cmro2 , and oef data and voi analyses were performed using in - house developed software ( idl , 6.2 , itt , boulder co , usa ) . individual cbf , cbv , cmro2 , and oef test  retest data will be displayed for all subjects . to assess intersubject variability of gray matter cbf , cbv , cmro2 , and oef measurements , standard deviation and coefficient of variation \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \left ( { { \hbox{cv } } = { \hbox{sd}}/{\hbox{mean } } } \right ) $ $ \end{document } were calculated across subjects for the defined flow territories . day - to - day reproducibility of gray matter cbf , cbv , cmro2 , and oef measurements as stated above , left - to - right ratios of hemodynamic parameters , in particular oef for the mca territory , may be important for selecting patients with compromised cerebral hemodynamics . therefore , left - to - right ratios of cbf , cbv , cmro2 , and oef were also calculated as well as their intersubject variability . test  retest pet studies were performed in ten healthy volunteers ( six men , four women ) with a mean age of 69  years ( range , 57 to 80  years ) who were scanned on two separate occasions ( with an interval of 354  days ) . in addition , all patients underwent magnetic resonance imaging ( mri ) , which was needed for pet image analysis shortly before the first pet study . the study was approved by the medical ethics committee of the university medical center utrecht and all volunteers provided informed consent . pet scans were acquired using an ecat exact hr+ scanner ( cti / siemens , knoxville , tn , usa ) which has been characterized elsewhere . subjects were scanned under baseline conditions , i.e. , lights in the room were dimmed and music was turned off . heads were positioned in an immobilization device to minimize patient motion during the study and in between scans . each study consisted of [ o]h2o , [ o]o2 , and [ o]co emission scans . prior to the emission scans , for the subsequent cbf study , a bolus ( ~5  s ) of 1,100  mbq [ o]h2o was administered intravenously , simultaneously starting a 3d dynamic emission scan ( 25 frames over a period of 600  s ) . the [ o]o2 study was performed 5  min after the end of the cbf study to allow for decay of residual o activity to background levels . a 3d dynamic emission scan ( 20 frames over a period of 600  s ) was then started simultaneously with the net inhalation of approximately 300  mbq [ o]o2 gas during a period of ~30  s. during administration of the [ o]o2 gas , the subjects continued breathing normally to avoid changes in cerebral hemodynamics , but they were instructed to inhale through the nasal administration cannula ( oxygen catheter dk3460 ; diameter , 4.7  mm ; unomedical , birkerd , denmark ) and exhale through the mouth . note that , in this study , scan durations of 600  s were used for both the water and oxygen studies . initially , it was attempted to use these longer scans to derive both perfusion and oxygen extraction fraction from a single [ o]o2 study . however , it was found that this was not feasible , and data were therefore analyzed by reusing perfusion and volume of distribution from the water scan ( see  image analysis  ) . although , by reusing parameters from the water scan , shorter scan durations may have been applied , it would not have shortened the overall study duration as some time interval between scans is needed to allow for radioactive decay . for the cbv study , a bolus of approximately 200  mbq [ o]co gas was administrated in a manner identical to the [ o]o2 procedure . one minute after the end of inhalation , a 3d dynamic emission scan ( three frames over a period of 360  s ) was started . all scans were reconstructed using fore  +  2d fbp reconstruction with a hanning filter at nyquist frequency , including all usual corrections for quantification . a matrix size of 256    256 and a zoom of 2.1 were applied , resulting in a voxel size of 1.2    1.2    2.4  mm and a final image resolution of ~7  mm full width at half maximum . for all dynamic scans , the arterial input function was measured using an online continuous blood - sampling device ; details on this device are given in . in short , the device consists of a bismuth germanate ( bgo)-based detector , shielded by a 3-cm lead , a withdrawal pump , and a computer for reading out the detector and controlling the pump . during blood withdrawal , a 1-mm - inner - diameter tubing that is connected to an arterial cannula is positioned in a slit that runs through the bgo detector . on the other end , the tubing is connected to a 3-mm - inner - diameter tubing that runs through the pump , and finally blood is collected in a 3-cm lead shielded waste bag . using scripts , both detector readout sampling and pump speed ( i.e. , withdrawal rates ) can set and automatically change during the studies . when necessary , pump speed can be overruled manually . during the first 4  min of each study , blood is withdrawn at a rate of 450  ml / h , and readout of the detector is performed continuously using 1-s readout samples . after 4  min onwards , pump speed is set to 150  ml / h to minimize the amount of blood withdrawn , and the detector readout sampling time is set to 10  s. manual samples were taken at discrete times for calibration purposes and for assessment of plasma - to - whole blood ratios : three samples at 5.5 , 8 , and 10  min during the [ o]h2o scan , six samples at 1 , 2 , 3 , 4.5 , 7 , and 9  min during the [ o]o2 scan , and three samples at 1 , 3 , and 5  min during the [ o]co scan . activity concentrations in the blood and plasma were determined for all samples . for the samples taken during the oxygen scan , the plasma - to - whole blood ratios were determined in order to separate the whole blood input function into water and oxygen input functions needed for kinetic analysis . for the latter purpose , a multi - exponential function was fitted through the six plasma - to - whole blood ratios . parametric cbf and the volume of distribution of water ( vt ) images were generated using a basis function implementation of the blood flow model , including corrections for dispersion , delay , and arterial blood volume as described previously . next , parametric cmro2 and oef images were generated using a slightly modified implementation of the method published by hattori et al .   with the following differences : ( 1 ) as scans were acquired over 600  s , all frames during this interval were included in fitting the data ; ( 2 ) the parametric vt data derived from the cbf scans were reused voxelwise rather than using the ( same ) fixed value for all voxels [ 16 , 25 ] ; and ( 3 ) the blood volume fraction was included as fit parameter rather than reusing it from the co study as it has been shown that this procedure provides more accurate fits near large blood vessels . more details about the implementation of the basis function method are given in the appendix . parametric cbv data were derived from the co scan as described by martin et al .   and hattori et al . . subsequently , an mri template with volumes of interest ( vois ) of the 3d flow territories of middle cerebral artery ( mca ) , anterior cerebral artery ( aca ) , and posterior cerebral artery ( pca ) ( fig . warping was performed using the automated image registration ( air ) software developed by woods et al . . nonrigid 12-parameter perspective warping was applied to fit the template onto the subject 's mri scan . afterwards , the warped , 3d vois were inspected and manually corrected when needed . statistical parametric mapping ( spm02 , wellcome department of cognitive neurology , london , natick , md , usa ) was used for gray  white matter segmentation of the coregistered mri scan . these gray and white matter segmentations were then used in combination with the warped voi template to calculate mean cbf , cbv , cmro2 , and oef in the different flow territories , separating data / voi into gray and white matter subregions . the generation of cbf , cbv , cmro2 , and oef data and voi analyses were performed using in - house developed software ( idl , 6.2 , itt , boulder co , usa ) . individual cbf , cbv , cmro2 , and oef test  retest data will be displayed for all subjects . to assess intersubject variability of gray matter cbf , cbv , cmro2 , and oef measurements , standard deviation and coefficient of variation \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \left ( { { \hbox{cv } } = { \hbox{sd}}/{\hbox{mean } } } \right ) $ $ \end{document } were calculated across subjects for the defined flow territories . day - to - day reproducibility of gray matter cbf , cbv , cmro2 , and oef measurements was calculated in a similar fashion . as stated above , left - to - right ratios of hemodynamic parameters , in particular oef for the mca territory , may be important for selecting patients with compromised cerebral hemodynamics . therefore , left - to - right ratios of cbf , cbv , cmro2 , and oef were also calculated as well as their intersubject variability . complete test and retest parametric cbf , cbv , cmro2 , and oef images could be generated reliably in seven of the ten subjects . in two subjects , the input curve of one of the scans ( i.e. , either test or retest scan ) showed irregularities due to radiation spillover , possibly due to poor positioning of the tubing within the detector and poor inhalation of o2 . in one subject , due to technical reasons and human mistakes , the time duration between the two scans was 54  days in one subject , 28  days in another subject , and 3 to 8  days in the other five subjects . figure  2 shows an example of the parametric images of a healthy subject in comparison with a patient with recurrent transient ischemic attacks ( tias ) associated with a left internal carotid artery occlusion . the cbf , cmro2 , and oef results for the seven successful test  retest studies are shown in fig .  3 and , together with cbv data , summarized in table  1 . left - to - right ratios of test and retest cbf , cmro2 , and oef measurements are shown in fig .  4 and , together with cbv data , summarized in table  2 . the range of the left - to - right oef ratio in the mca gray matter territory was 0.981.03 ( n  =  14 ) . 2typical cbf , oef , and cmro2 parametric maps for a healthy volunteer in comparison with a patient who had recurrent tias associated with a left internal carotid artery occlusion . the pet study of the patient shows in the ipsilateral hemisphere a decreased cbf and increased oef compared with the contralateral hemisphere.fig . 3relationship between test and retest measurements of gray matter cbf , cmro2 , and oef in seven healthy volunteers . black circle test right hemisphere ; black triangle test left hemisphere , black square test both hemispheres . gray circle retest right hemisphere , gray triangle retest left hemisphere , gray square retest both hemispheres . total brain included all flow territories.table  1intersubject variability and day - to - day reproducibility of dynamic gray matter pet measurementsintersubject variabilityday - to - day variabilitymean  sdsdcv ( % ) sdcv ( % ) cbf ( ml / min/100  ml )  mca  left37.4    3.35.514.63.38.9  right37.3    3.35.815.23.49.1  pca  left41.2    3.86.314.64.19.7  right41.3    3.66.615.53.37.8  aca  left37.1    2.85.915.43.910.6  right37.9    3.05.614.43.28.3  total brain37.6    2.95.113.13.38.8cbv ( ml/100 ml )  mca  left2.96    0.50.413.60.413.8  right2.99    0.50.412.00.414.1 pca  left3.35    0.50.515.10.515.9  right3.47    0.60.615.80.516.1  aca  left2.79    0.40.27.90.415.9  right2.90    0.50.618.20.414.9  total brain4.06    0.60.512.10.513.8cmro2 ( ml / min/100  ml )  mca  left3.04    0.20.26.90.14.1  right3.01    0.20.27.20.14.5  pca  left3.39    0.30.37.70.25.3  right3.38    0.30.39.00.25.4  aca  left2.92    0.20.38.80.26.9  right2.96    0.30.39.20.13.8  total brain3.05    0.20.27.10.25.3oef ( % )  mca  left43.4    6.59.522.04.09.8  right43.1    6.59.221.64.19.5  pca  left43.6    5.68.419.64.19.5  right43.4    6.39.021.03.99.3  aca  left41.8    6.69.623.14.310.5  right41.3    6.08.821.54.19.9  total brain43.0    6.39.121.44.09.3sd test  retest standard deviation , cv coefficient of variation , mca middle cerebral artery , pca , posterior cerebral artery , aca anterior cerebral arterybaseline datamean test and retest data . aca anterior cerebral artery , mca middle cerebral artery , pca posterior cerebral artery.table  2intersubject variability and day - to - day reproducibility of left right ratios of dynamic gray matter pet measurementsintersubject variabilityday - to - day reproducibilitymean  sdsdcv ( % ) sdcv ( % ) cbf ( ml / min/100  ml )  l / r ratio  mca1.00    0.020.022.20.011.1  pca1.00    0.030.032.80.032.9  aca0.98    0.040.087.90.044.5cbv ( ml/100  ml )  l / r ratio  mca0.99    0.010.033.40.032.6  pca0.97    0.080.099.40.033.5  aca0.97    0.120.1312.70.077.4cmro2 ( ml / min/100  ml )  l / r ratio  mca1.01    0.020.021.80.021.6  pca1.00    0.040.032.90.033.0  aca0.99    0.040.087.80.044.0oef  l / r ratio  mca1.01    0.010.021.70.011.5  pca1.01    0.020.022.10.010.8  aca1.01    0.020.022.20.011.4sd test  retest standard deviation , cv coefficient of variation , mca middle cerebral artery , pca posterior cerebral artery , aca anterior cerebral arterybaseline datamean test and retest data . mean left right ratios given for comparison typical cbf , oef , and cmro2 parametric maps for a healthy volunteer in comparison with a patient who had recurrent tias associated with a left internal carotid artery occlusion . the pet study of the patient shows in the ipsilateral hemisphere a decreased cbf and increased oef compared with the contralateral hemisphere . relationship between test and retest measurements of gray matter cbf , cmro2 , and oef in seven healthy volunteers . black circle test right hemisphere ; black triangle test left hemisphere , black square test both hemispheres . gray circle retest right hemisphere , gray triangle retest left hemisphere , gray square retest both hemispheres . intersubject variability and day - to - day reproducibility of dynamic gray matter pet measurements sd test retest standard deviation , cv coefficient of variation , mca middle cerebral artery , pca , posterior cerebral artery , aca anterior cerebral artery mean test and retest data . retest standard deviation , cv coefficient of variation , mca middle cerebral artery , pca posterior cerebral artery , aca anterior cerebral artery mean test and retest data . this study provides reference data for cbf , cbv , cmro2 , and oef values derived from dynamic o pet scans . obtained mean values of cbf , cbv , cmro2 , and oef are comparable to those reported previously [ 9 , 29 , 30 ] . compared with the within - session reproducibility of steady - state measurement of cbf and oef , we found a similar day - to - day test  retest variability of cbf and a higher day - to - day test  retest variability of oef measurements . day - to - day variability of cbf measurements in the present study was comparable to that reported previously for healthy volunteers [ 2 , 31 , 32 ] . in addition to measurement errors that affect day - to - day variability of cbf , cbv , cmro2 , and oef measurements , the day - to - day test  retest variability of cbf measurements in the present study is most likely also affected by physiological day - to - day variability . with the present study design , however , it is not possible to distinguish between measurement errors and physiological day - to - day variance . although the exact contribution of physiological and measurement errors could be studied further using a within - session test retest reproducibility study , the present design provides a clinically relevant assessment of the variability in cbf , cbv , cmro2 , and oef measurements for follow - up studies ( e.g. , when monitoring the effects of treatment ) . interestingly , the variability of cmro2 measurements was low and comparable to within - session reproducibility of cmro2 derived from steady - state studies . the relatively high day - to - day variability in cbf and low day - to - day variability in cmro2 also implies a high variability in oef . this finding is in line with previous observations where artificial changes in cbf were compensated by opposite changes in oef   to maintain constant cmro2 needed to preserve cell functions of brain neurons . retest variability of cmro2 validates the use of a dynamic scanning procedure for o pet studies . other o pet techniques have been described in literature , such as autoradiographic methods and dual autoradiographic methods [ 18 , 36 ] and steady - state measurements [ 7 , 1012 ] . a possible advantage of the method used in our study is that it is not required to assume a fixed value for the volume of distribution ( or partition coefficient ) of water and that the blood volume fraction is used as fit parameter . in theory , this may provide more accurate results but at the cost of using two additional fit parameters in the calculations . the test  retest variability and cbf , cmro2 , and oef values observed in our studies , however , are very comparable to those found elsewhere using different scanning techniques as discussed above . moreover , the very low variability seen in cmro2 suggests that the higher day - to - day variabilities observed for cbf and oef are likely determined by physiological variation rather than technical issues or the applied o pet technique . retest variabilities are independent of the applied o pet technique or , at least , that the dynamic scanning technique performs equally as well as others . clearly , critical issues in the dynamic scanning technique are related to the collection of the arterial input function as illustrated by the three study failures . also , shorter scan durations than 10-min acquisitions may have been applied , yet it would not have shortened the overall study duration as some time interval between scans is needed to allow for radioactive decay . moreover , time is needed for quality control of the collected input function data and for repeat production of the o - labeled tracers when needed . when studying left - to - right ratios of cbf , cmro2 , and oef ( table  2 ) , excellent test  retest variabilities were observed , i.e. , in the order of a few percent ( cv ) . this good test retest variability was less prominent for the cbv left - to - right ratio , probably due to higher noise levels in the [ o]co scan and asymmetry in the presence of larger blood vessels in the 3d voi . in the kinetic analysis of cbf , cmro2 , and oef data , the vascular component within a voi was used as fit parameter rather than reusing cbv derived from the [ o]co study . therefore , differences or inaccuracies in vascular fraction seen in the [ o]co scans did not affect regional test  retest variability of cbf , cmro2 , and oef measurements . apart from studying absolute changes in cbf , cmro2 , and oef in case of monitoring response to treatment or following pathophysiological processes over time , dynamic o pet studies may be used clinically as a sensitive means of detecting ipsilateral versus contralateral asymmetries in cbf , cmro2 , and oef . the present values of left - to - right ratios were comparable with previous measurements . taken together , the present study provides additional reference data for absolute and relative ( i.e. , left - to - right ratios ) cbf , cmro2 , and oef measurements using dynamic o pet scans . when asymmetry in cbf or oef is expected such as in patients with unilateral ica occlusion who are at high risk of stroke , left - to - right ratios may be useful because of their excellent test  retest variability ( table  2 ) . it should be emphasized that in our study results were obtained at a baseline condition . in physiologically stressed conditions , retest variability may be similar to those seen under baseline conditions , i.e. , effects of image noise and input function collection and processing on cbf , oef , and cmro2 test yet day - to - day test  retest variability may be more strongly affected by the quality of obtaining reproducible stressed conditions . from the present data , it is not possible to determine the test  retest variability under physiologically stressed or pathological conditions , and this warrants further investigations , yet information on test  retest variability under baseline conditions provide valuable information of the normal range of cbf , oef , and cmro2 , and this information is needed in order to reliably identify subjects or brain regions with abnormally high or low values . in the present study , warped 3d vois were used to measure cbf , cmro2 , and oef in entire mca , pca , and aca flow territories . it prevents possible selection errors that can occur when only a few smaller regions of interest are used for analysis . moreover , as carotid occlusions are likely to affect cbf , cmro2 , and oef in entire flow territories , use of 3d voi will provide data with better statistics and higher sensitivity . however , use of larger 3d voi may underestimate and even miss small focal changes in cbf , cmro2 , and oef because of the relatively large volumes of the flow territories . it should be noted , however , that in the present study first parametric maps were generated . clearly , these parametric images could also be used to analyze the data with smaller voi or at the voxel level . as mentioned above , three out of 20 pet studies could not be analyzed primarily because of errors in the acquired input function data . correct administration of the tracers and collection of the input functions seem to be the most critical issues in obtaining reliable results . because of the short half - life of o , however , acquisition errors could easily be corrected provided that they are detected immediately after or during the scanning procedure . in those cases , a scan could be repeated , adding only 15  min to the overall study duration . therefore , for routine clinical use , it is required to have an online data quality assurance procedure in place . at present , we have implemented an online quality control system based on the number of counts ( activity ) seen in the manual samples to verify sufficient and correct administration of the tracers . furthermore , online review of the input functions allows verification of correct collection of the continuous online blood input curves . day - to - day test  retest variability of cmro2 measurements using dynamic o pet scans is comparable to within - session test  retest variability using steady - state [ o ] pet scans . day - to - day test  retest variability of cbf and oef measurements probably is determined primarily by day - to - day physiological variability in cbf and corresponding adaptation of oef . test retest variability in left - to - right ratios of cbf , cmro2 , and oef is excellent ( cv  <  4.5% ) and may even be used to detect impaired brain perfusion and/or oxygen consumption in individual patients . parametric cbf and volume of distribution of water ( vt ) images were generated using a basis function implementation of the blood flow model , including corrections for dispersion , delay , and arterial blood volume as described previously ( eqs . the model is given by a single tissue compartment model with ( arterial ) blood volume fraction ( va ) correction : 
 1\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ { c_{\rm{t } } } = f \bullet { c_{\rm{a } } } \otimes { e^ { - ( f/{v_{\rm{t } } } ) \bullet t } } $ $ \end{document}2\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ { c_{{\rm{roi / vox } } } } = ( 1 - { v_{\rm{a } } } ) \bullet f \bullet { c_{\rm{a } } } \otimes { e^ { - ( f/{v_{\rm{t } } } ) \bullet t } } + { v_{\rm{a } } } \bullet { c_{\rm{a } } } $ $ \end{document}ct represents the activity concentration in tissue and croi / vox the activity concentration in a region of interest or , in this case , in a voxel . f is the flow or perfusion , vt is the volume of distribution , and ca is the arterial input function ( in this case of arterial whole - blood activity concentration over time ) . in order to reduce computation time , the convolution term in these equations can be precomputed for a range of \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \left ( { f/{v_{\rm{t } } } = { k_{\rm{i } } } } \right ) $ $ \end{document } values and then eq .  1 can be rewritten as a linear equation : 
 3\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ { \hbox{b}}{{\hbox{f}}_{\rm{i}}}(t ) = { c_{\rm{b } } } \otimes { e^ { - f/{v_{\rm{t } } } \bullet t } } = { c_{\rm{b } } } \otimes { e^ { - { k_{\rm{i } } } \bullet t } } $ $ \end{document}and eq .  2 is then rewritten as ( eq .  4 ) : 
 4\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ { c_{{\rm{roi / vox } } } } = ( 1 - { v_{\rm{a } } } ) \bullet f \bullet { \hbox{b}}{{\hbox{f}}_{\rm{i } } } + { v_{\rm{a } } } \bullet { c_{\rm{a } } } $ $ \end{document}finally , for each voxel time  activity curve , the combination of bfi(t ) ( and thus ki ) , f , and va which provides the minimal weighted least squares discrepancy between measured data and eq . next , parametric cmro2 and oef images were generated using a slightly modified implementation of the method published by hattori et al .   with the following differences : ( 1 ) all frames over the interval of 0600  s postinjection were included , ( 2 ) the parametric vt data derived from the cbf scans were reused voxelwise rather than using the ( same ) fixed value for all voxels [ 16 , 25 ] , and ( 3 ) arterial blood volume was used as fit parameter rather than reusing it from the co study as it has been shown that this procedure provides more accurate fits near large blood vessels . the basis function method for analyzing the oxygen studies is thus given by ( eq .  5 ) : 
 5\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ { c_{\rm{t } } } = { v_{\rm{a } } } \bullet c_{\rm{a}}^{{{\rm{o}}_2 } } + ( 1 - { v_{\rm{a } } } ) \bullet [ f \bullet c_{\rm{a}}^{{{\rm{h}}_2}{\rm{o } } } \otimes { e^ { - kt } } + f \bullet e \bullet c_{\rm{a}}^{{{\rm{o}}_2 } } \otimes { e^ { - kt } } ] $ $ \end{document}croi / vox is the pet tissue time  activity curve , corrected for contribution of activity from venous blood using method of holden et al . . vais the arterial blood volume fraction or spillover correction term ( needed in parametric application of the model ) . \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ { c_{\rm{a}}}^{{{\rm{h}}_2}{\rm{o } } } $ $ \end{document}is the arterial input function of water . \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ { c_{\rm{a}}}^{{{\rm{o}}_2 } } $ $ \end{document}is the arterial input function of oxygen . f is the flow or perfusion , e is the oxygen extraction fraction ( oef ) , and \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ k = f/{v_{\rm{t}}},\;{v_{\rm{t } } } = { \hbox{volume}}\;{\hbox{of}}\;{\hbox{distribution}}\;{\hbox{of}}\;{\hbox{water } } $ $ \end{document}. the total arterial input function ca is split into \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ { c_{\rm{a}}}^{{{\rm{h}}_2}{\rm{o } } } $ $ \end{document } and \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ { c_{\rm{a}}}^{{{\rm{o}}_2 } } $ $ \end{document } , based on the activity concentration ratios between whole blood \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \left ( { = \left [ { ^{{15}}{\hbox{o } } } \right]{{\hbox{h}}_2}{\hbox{o } } + \left [ { ^{{15}}{\hbox{o } } } \right]{{\hbox{o}}_2 } } \right ) $ $ \end{document } and plasma \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \left ( { = \left [ { ^{{15}}{\hbox{o } } } \right]{{\hbox{h}}_2}{\hbox{o } } } \right ) $ $ \end{document } sample data , thus assuming that oxygen is bound in red blood cell , i.e. , [ o]o2 is not present in plasma . in our implementation , f and vt and thus the k values are reused voxelwise from the parametric data derived from the o water study . consequently , the convolution terms are precomputed / known for each voxel , and only a linear equation remains . by weighted linear least squares fitting of the voxel tac with eq .  5 , the oef ( = e ) and thus also the product of oef and f are obtained . cmro2 is then derived simply by multiplying the oef by the f image with the oxygen content in blood ( derived from the manual blood samples ) .
purposewe assessed test  retest variability of cerebral blood flow ( cbf ) , cerebral blood volume ( cbv ) , cerebral metabolic rate of oxygen ( cmro2 ) , and oxygen extraction fraction ( oef ) measurements derived from dynamic 15o positron emission tomography ( pet ) scans.proceduresin seven healthy volunteers , complete test  retest 15o pet studies were obtained ; test  retest variability and left - to - right ratios of cbf , cbv , oef , and cmro2 in arterial flow territories were calculated.resultswhole-brain test  retest coefficients of variation for cbf , cbv , cmro2 , and oef were 8.8% , 13.8% , 5.3% , and 9.3% , respectively . test 
 retest variability of cbv left - to - right ratios was < 7.4% across all territories . 
 corresponding values for cbf , cmro2 , and oef were better , i.e. , < 4.5% , < 4.0% , and < 1.4% , respectively.conclusionsthe test  retest variability of cmro2 measurements derived from dynamic 15o pet scans is comparable to within - session test 
 retest variability derived from steady - state 15o pet scans . 
 excellent regional test 
 retest variability was observed for cbf , cmro2 , and oef . 
 variability of absolute cbf and oef measurements is probably affected by physiological day - to - day variability of cbf .