Quantification of the oxygen saturation of hemoglobin within the vasculature of the eye would be useful in the study and management of a number of diseases associated with regional ischemia and infarct such as diabetic retinopathy,1 ocular vasculopathies such as anterior ischemic optic neuropathy,2 and central and retinal branch artery and vein occlusion.3 For this reason, there have been numerous attempts to develop a spectral technique for the assessment of retinal blood oxygenation dating back to the early . 4, 5, 6, 7, 8, 9, 10 Preliminary data suggests that spectral techniques may be applicable in recently developed spectroscopic optical coherence tomography technology.11, 12, 13
All spectral reflectance methods are based on assumptions about the relationship between light transmittance and oxygen saturation,14 which is expressed asthrough blood is equal to the original intensity multiplied by raised to the factor of a constant times distance traveled through blood times concentration of blood . Solving for yields is the extinction coefficient, which is a term quantifying the combined effects of scattering and absorption per unit distance per unit concentration in a solution. Optical density (OD) is defined as the natural logarithm of .
Ultra-high-resolution 15, 16, 17, 18, 19, 20, 21, 22 and Fourier domain23, 24, 25, 26, 27 optical coherence tomography (OCT) are relatively new technologies that use a broad-spectrum light source. The devices acquire high-resolution images at high speed, which provides highly detailed information from the scanned area. Unlike other application-specific devices such as pulse oximeters, the retinal vessel analyzer that measures blood vessel diameter, or the Canon laser blood flowmeter that measures blood flow in specific retinal vessels, OCT scans of the retina have the potential to simultaneously assess tissue volume, blood flow within that tissue bed, and the oxygen saturation of the blood flowing within those vessels. Combined, OCT scans of the human retina have the potential to provide the only complete and precise assessment of metabolism in medicine. Our high-speed ultra-high-resolution Fourier domain OCT (hsUHR-OCT) uses a light source that includes an oxyhemoglobin-deoxyhemoglobin isosbestic wavelength. Given the previous success of fundus camera-based systems, in which reflectance was quantified with no tissue depth resolution, an analysis of full axial content a-scans may provide sufficient data to quantify an oxygen sensitive optical density ratio (ODR).6, 10 The purpose of the present study is to determine if the spectral data within hsUHR-OCT scans are suitable for the measurement of blood oxygen saturation measurements within retinal arteries and veins.
All procedures were approved by the institutional review board at the University of Pittsburgh School of Medicine. Informed consent was provided by all subjects before participation in the studies, and the study was conducted in compliance with the health insurance portability and accountability act.
Seventeen subjects were recruited for the study. All subjects were recruited from the staff and faculty of the Eye and Ear Institute of the University of Pittsburgh Medical Center, and were normal healthy individuals with no history of ocular disease or trauma.
The hsUHR-OCT may be used to image the layers of the retina (Fig. 1 ). Principles of the hsUHR-OCT have been described previously.23, 24, 28, 29 Briefly, light from a bandwidth low-coherence broadband light source ( , Broadlighter 840, SuperLum Diodes Ltd., Moscow, Russia) is split into two paths in a 50:50 fiber coupler. One beam is projected onto a stationary reference mirror and the second was focused on and scanned the retina (Fig. 2 ). Light reflected from the mirror and the light backscattered from the retina are combined to form an interference pattern. Depth information describing the location of scattering structures within the retina is frequency-encoded in the spectral fringes.
Three dimensional (3-D) optic-disk-centered raster scans were obtained with the hsUHR-OCT device. Scans encompassed a slice of retinal tissue, with a total axial dimension of . The area was scanned with a -point array of a-scans (Fig. 1). Each a-scan contained 1024 axial pixels, and the axial resolution enabled by the source bandwidth was . Image quality was assessed subjectively and only good quality scans were used for the analysis. Specifically, OCT fundus images were assessed for the appearance of sharp vessel borders, suggesting good focus, and evenly spread vignetting, suggesting that the scan was positioned on the center of the field of view with minimal distortion due to eye movements or black lines due to blinking.
Spectral Oximetry and Data Analysis
Data describing reflectance at specific wavelengths were obtained from structural optic-nerve-head (ONH)-centered hsUHR-OCT image data. The wavelengths measured by the CCD camera (Fig. 2) were calibrated by locating peaks within the light source spectrum provided by the manufacturer (Fig. 3 ). Customized code written in Labview (version 7.1, National Instruments, Inc., Austin, Texas) was used to identify the pixel location and bandwidth of the wavelengths of the light source. The light from our bandwidth source was focused onto the central of a linear CCD camera, providing a spectral resolution with a theoretical limit of approximately . Unlike the numerous spectral OCT systems utilizing a light source with a bandwidth, the device in this study uses a light source with a bandwidth of . This necessitates the sharpest possible focus of the spectral signal onto the linear CCD, and minimal spectral aberration. The utilization of a subset of all available pixels on the CCD camera is due to the limitation in the optics available at the time of device construction. Sharp focus could not be obtained at any greater focal length, so the physically limited width is used. In order to sample light with a range of a wavelength of interest, the surrounding that wavelength were analyzed. Zero-padding provided a-scans of for each narrow band of wavelengths sampled. For example, to measure light at , the pixels from 620 to 715 were sampled.
Determination of hemoglobin oxygen saturation requires quantification of the amount of light absorbed or scattered by blood within the retinal vasculature.11 This was accomplished by comparing a-scans at the edges of blood vessels to adjacent a-scans off of the vessel (Fig. 4 ). Artery and vein pairs were located within the 2-D projection of the OCT fundus image. The unprocessed fringe data corresponding to pixels chosen for analysis within the OCT fundus image was identified. A-scans were recalculated, zeroing the constant term in the power spectrum, and used for analysis. Sampling vessels and tissue in close proximity ensured similar tissue thickness and pigment content between the two, leaving differences in reflections primarily due to interaction with blood, or the lack thereof. A-scans at the edge of the blood vessels, as opposed to the center of the blood vessels, were chosen to reduce the influence of specular reflection, and the use of the OCT fundus images allows identification of specific a-scans of interest. Light from the entire a-scan was summed in the spectral analysis with the assumption that the difference between total reflectance within these two a-scans was due primarily to light being scattered or absorbed by blood.30 This approach mimics previous fundus camera techniques, and helps to remove the effects of local variations in absorption by other photoactive materials and the contributions of local tissue geometry to variations in reflectance.
Thirty sets of fringe data were sampled from each image. Ten fringe patterns were sampled from the edge of an artery, adjacent tissue, and the edge of a vein, respectively. Average light intensities for each of the 10 measurements for each subject were calculated, as were average intensities across all subjects for artery, vein, and tissue. Optical densities (ODs) and optical density ratios (ODRs) were calculated for each. Arteries and veins were identified by visual inspection of the OCT fundus image. It was required that sampled pixels be in an area of vessel or tissue with adequate illumination. As previously mentioned, in the case of vessel samples, pixels were required to be near the edge of the vessel in order to avoid specular reflection. Vessels were sampled in the plane of the fundus, outside of the ONH. OD is a measure of light absorbance relative to tissue surrounding a blood vessel,31 and is calculated asand are the average pixel values inside and immediately adjacent to a blood vessel when illuminated at wavelength . The extinction coefficients of oxyhemoglobin and deoxyhemoglobin vary with wavelength (Fig. 5 ).32 For example, light at is equally scattered by oxyhemoglobin and deoxyhemoglobin, and for this reason is defined as an isosbestic wavelength. Ratios of ODs (ODR) between those measured at oxygen sensitive wavelengths and isosbestic wavelengths are linearly proportional to blood oxygen saturation.30 This is expressed as: is percent hemoglobin oxygen saturation, and are constants, and are optical densities at oxygen sensitive and insensitive wavelengths, respectively, and ODR is the optical density ratio. In the present study, was used as the isosbestic wavelength, and was used as the oxygen sensitive wavelength. Here, was selected because it is a location in the spectrum where the difference in extinction coefficients between oxyhemoglobin and deoxyhemoglobin is relatively large within the available bandwidth.
Mean light intensities for each parameter for each subject were calculated. The coefficients of variation for light intensities for each parameter for each subject were calculated. Ninety-five pixels presenting data were summed to measure light intensity in the bandwidths of interest, representing a potential intensity signal of 389,120. In order to obtain an estimate of signal strength, the average intensity measurements for all parameters as a percentage of CCD light saturation level was also calculated. Arterial and venous ODRs, calculated from the average intensities across all subjects, were compared by the Wilcoxon signed rank test. values less than or equal to 0.05 were considered significant.
Seventeen normal healthy subjects (nine female, eight male, ages ) were recruited for the study. Ophthalmic normalcy had been established by subject participation in prior studies. All subjects were successfully examined by hsUHR-OCT using the 3-D protocol. Images were of sufficient quality to allow sampling of on an artery, on a vein, and on tissue between an artery and vein, while avoiding motion artifacts (Fig. 6 ). Eleven of 17 eyes presented with higher oxygen saturation observed in the artery. The mean values of the averages of the light intensity within the respective bandwidths (isosbestic and oxygen sensitive) for each of the locations are presented in Tables 1, 2 for each subject. Calculated OD values calculated for the perspective bandwidths for each location (tissue, artery, and vein) are presented in Tables 1, 2. Arteriovenous ODR differences for each subject are listed in Tables 1, 2, with positive values representing arterial blood oxygen saturation higher than venous oxygen saturation. The coefficient of variation for each of the intensity measurements for each subject was calculated; values are presented in Tables 1, 2 for each subject. Arterial ODRs were significantly greater than venous ODRs ( and respectively; ) ( deviation), suggesting higher oxygen saturation in arterial blood (Fig. 7 ). The average arteriovenous (AV) ODR differences were ( deviation), ranging from to 13.73.
A continuation of Table 1, with mean and standard deviations presented for each of the variables across all subjects. The mean light intensity measurements as a percentage of CCD saturation are presented as an indication of signal intensity.
|11||12||13||14||15||16||17||All SubjectsMean||All SubjectsStandard Deviation||% TotalPotential Light|
|Mean Artery Isosbestic||41.59||22.54||61.26||37.98||53.80||49.74||45.34||45.30||14.30||0.012%|
|Mean Artery Sensitive||48.75||3.43||73.79||37.92||59.82||46.19||49.75||50.19||20.49||0.013%|
|Mean Vein Isosbestic||40.29||21.26||63.73||36.21||57.87||47.39||45.33||44.47||15.18||0.011%|
|Mean Vein Sensitive||47.47||3.64||72.68||39.59||58.67||50.18||50.07||50.08||20.42||0.013%|
|Mean Tissue Isosbestic||42.05||16.44||64.90||41.44||53.34||47.27||43.72||43.41||15.03||0.011%|
|Mean Tissue Sensitive||51.19||3.74||72.66||45.66||61.21||52.24||50.85||50.94||19.65||0.013%|
|OD Artery Isosbestic||0.01||0.06||0.09||0.11|
|OD Artery Sensitive||0.05||0.09||0.19||0.02||0.12||0.02||0.02||0.09|
|OD Vein Isosbestic||0.04||0.02||0.14||0.00||0.09|
|OD Vein Sensitive||0.08||0.03||0.00||0.14||0.04||0.04||0.02||0.03||0.08|
|ODR Art—ODR Vein||2.68||1.07||13.73||2.44||4.01|
|Coefficient of Variation Artery Isosbestic||19.46%||46.32%||19.37%||23.88%||28.19%||24.82%||28.14%|
|Coefficient of Variation Artery Sensitive||17.94%||111.78%||27.29%||28.69%||16.23%||20.72%||17.29%|
|Coefficient of Variation Vein Isosbestic||13.00%||30.11%||27.39%||34.11%||26.14%||29.65%||27.64%|
|Coefficient of Variation Vein Sensitive||24.92%||182.92%||22.62%||46.06%||25.54%||25.38%||26.28%|
|Coefficient of Variation Tissue Isosbestic||9.27%||62.73%||20.88%||28.05%||28.11%||19.58%||15.79%|
|Coefficient of Variation Tissue Sensitive||16.98%||147.59%||19.27%||30.89%||26.43%||36.29%||12.98%|
The “mean” values are the averages of the light intensity within the respective bandwidths (isosbestic and oxygen sensitive). OD values are calculated as the natural log of the ratio of tissue to vessel light intensity for the perspective bandwidths. Positive ODR Art–ODR Vein differences represent subjects in which arterial blood oxygen saturation was higher than venous oxygen saturation. Coefficients of variation for each of the parameters, calculated from the mean and standard deviation of the 10 measurements, are presented.
|Mean Artery Isosbestic||19.14||48.88||61.74||45.87||52.58||72.41||34.18||56.02||27.71||39.24|
|Mean Artery Sensitive||19.23||60.42||82.47||69.90||45.56||75.14||35.45||63.68||39.80||41.97|
|Mean Vein Isosbestic||15.81||48.69||51.32||50.77||54.18||66.92||32.64||63.39||24.38||35.78|
|Mean Vein Sensitive||14.46||68.68||68.92||71.16||43.68||72.66||37.56||71.51||36.02||44.32|
|Mean Tissue Isosbestic||15.26||44.30||55.51||47.27||50.89||65.33||34.44||59.04||23.98||32.84|
|Mean Tissue Sensitive||18.67||61.00||67.93||73.05||41.79||73.77||37.70||71.93||38.24||44.40|
|OD Artery Isosbestic||0.03||0.01||0.05|
|OD Artery Sensitive||0.01||0.04||0.06||0.12||0.06|
|OD Vein Isosbestic||0.08||0.05|
|OD Vein Sensitive||0.26||0.03||0.02||0.00||0.01||0.06||0.00|
|ODR Art—ODR Vein||7.30||2.01||1.83||1.94||0.81||8.31||2.40||3.85|
|Coefficient of Variation Artery Isosbestic||23.57%||17.38%||16.88%||23.64%||13.92%||19.82%||14.21%||22.52%||25.59%||25.52%|
|Coefficient of Variation Artery Sensitive||32.98%||20.22%||21.68%||22.70%||18.63%||19.29%||15.31%||26.61%||18.31%||18.20%|
|Coefficient of Variation Vein Isosbestic||46.84%||29.46%||15.70%||25.17%||17.90%||7.22%||14.85%||16.32%||35.48%||32.90%|
|Coefficient of Variation Vein Sensitive||46.20%||25.78%||20.31%||22.78%||22.22%||14.03%||16.74%||26.04%||20.29%||24.14%|
|Coefficient of Variation Tissue Isosbestic||31.18%||23.54%||17.22%||16.74%||17.75%||13.56%||10.59%||18.96%||15.82%||22.22%|
|Coefficient of Variation Tissue Sensitive||39.06%||20.28%||20.01%||11.83%||19.19%||14.56%||11.26%||22.84%||23.47%||22.34%|
In the present study, we found that significant differences can be observed in spectral absorption between blood in arteries and veins in the back of the eye. Specifically, arterial ODRs were significantly greater than venous ODRs. One contributing factor to the observed difference may be a detectable higher concentration of oxygenated hemoglobin in retinal arteries than in retinal veins. However, since the entire a-scan was included in the analysis, it is possible that local variations in tissue pigmentation or any number of other factors existing outside of the blood column could have contributed to either the observed significant difference, or the high levels of variability observed. Indeed, until the contribution to noise and possible systematic error induced by the various steps involved in the processing and analysis of the data are quantified, it is impossible to confidently conclude that hemoglobin oxygenation was the driving factor in the observed significant difference. Nevertheless, the present data do agree with the physical reality of higher oxygen saturation in arteries than in veins.
Another source of potential noise in the measurement is differing tissue content between a-scans primarily containing tissue and those containing blood large vessels. Figure 5 suggests that the retinal nerve fiber layer (RNFL) layer in tissue adjacent to large blood vessels may be thicker than the RNFL content in the volume of an a-scan that includes a large blood vessel. This difference in the content of backscattering materials in a-scans that the present study assumes to be similar may be contributing to the level of variability in the data.
A-scans were obtained near the vessel edges to avoid specular reflection that occurs near the center (Fig. 6). The intensity of specular reflection decreases as the OD measurements are obtained farther from the vessel center. However, the distance through the blood column also decreases as the measurement distance from the vessel center increases. It is currently unknown if there is some critical level of specular reflection that must be avoided for meaningful OD measurements. For this reason, it is advised that OD measurements of blood columns be made near the edge of the vessel.33
The utility of the measurements in the present study for oximetry depends on the linear relationship between ODR and oxygen saturation10; however, the applicability of ODR to oximetry measurements obtained from backscattered light has several limitations.34 The Beer-Lambert law applies to transmission of light; and furthermore, the extinction coefficient values for light transmission through a solution of hemoglobin were obtained using hemolyzed blood.34 The cells in whole blood as well as the surrounding tissue all create scatter, which is not accounted for in the Beer-Lambert model. Further work is necessary to provide more accurate estimates of the “extinction coefficient” as measured under conditions of light backscattered from whole blood in vessels embedded in tissue and illuminated with a low-coherence superluminescent diode, as is the case in OCT imaging. Concerning the use of the OCT to estimate the intensity of backscattered light, previous techniques utilized film or digital cameras to integrate light backscattered from all reflective retinal sources regardless of tissue depth. The OCT limits the location of reflective signal sources to those contained within the a-scan. Measurement of backscattered light from nearly identical depths, as obtained when analyzing adjacent a-scans in an OCT frame, may be beneficial.
Coefficients of variation in intensity measurements were high, ranging from 7% to 183%, averaging 26%. For example, subject 12 presents with very large coefficients of variation. In this subject, the light levels in at least some of the measurements were very low (Table 2). While the average coefficient of variation was below 30%, the high variability in the level of variation may be due to the extremely low signal level. The average level of light intensity measurement was barely above of 1% of the CCD saturation level.
The analysis of light backscattered by shadowed tissue contributes to the signal in at least three ways. First, light traveling to shadowed tissue must pass through the blood column. This light will have experienced scatter and absorption as a function of wavelength, altering the spectral intensity distribution before arriving at the tissue. Second, light backscattered from shadowed tissue must pass through the blood column a second time on the way to the spectrometer. This light will also experience wavelength-dependent scatter and absorption. Finally, unlike light backscattered directly from the blood column, it is possible that the optical qualities of the tissue may also alter the spectrum of backscattered light in a wavelength-dependent manner. This would be damaging to attempts to calculate an exact ODR based on the arterial or venous a-scan alone; however, the ODs calculated in the present study used light backscattered from adjacent a-scans containing nonshadowed tissue to represent the expected level of backscattered light in the absence of blood interactions. In this case, the local optical properties of shadowed and nonshadowed tissues are assumed to be equal, and their effect removed in calculation of the OD ratio. In this context, the inclusion of shadowed tissue in the calculation of OD actually amplifies the effect of light/blood interaction by making measurements of backscattered light that has twice passed through the blood column.
The ultimate goal of Fourier domain OCT oximetry is to determine a precise numerical value for the percentage of oxygen saturation of blood within retinal arteries and veins. Unfortunately, these vessels are located adjacent to the photoreceptors and retinal pigment epithelium, each containing various pigments that may affect optical density measurements, and behind the vitreous. Light going to and from the vessels must interact with water within the vitreous (relevant in the near-IR range) and retinal pigments. In a previous study of retinal blood oxygenation measured in a two-wavelength spectral system (569 and ), arterial blood oxygen saturation was 97% and venous blood oxygen saturation was 55%.30 In this study the authors corrected ODRs for individual influences of retinal pigments. The influence of spectral absorption in the near-IR wavelengths used in the present study has yet to be determined. It is possible that interindividual variations in pigments may have contributed to the variability observed in the absolute values in ODRs and differences in magnitude in AV ODR differences., it bears repeating that the extinction coefficient combines both absorption and scattering effects. With , , , and assuming a 55% saturation in venous blood 30 and 100% saturation in the arteries, arterial ODRs should be approximately 1.2 times greater than in veins. Our data demonstrated a wide variation in the difference between arterial and venous ODRs. It is possible that this is due to the effects of pigments and water absorption, which were not corrected for, or it may also be due to a variation in venous blood saturation within this experimental cohort.
Differential absorption by water also contributed to differences in ODs between 805 and . The extinction coefficient of water increases by approximately 300% in the range, the bandwidth of our light source. The effect of water was constant between arteries and veins, since the approximate distance through the vitreous and anterior structures to adjacent vessels and tissue is the same; nevertheless, the contribution of water absorption will need to be taken into account before actual saturation values can be calculated.
In conclusion, this pilot study utilized a two-wavelength ODR approach similar to a previous publication with a fundus camera-based system. Calculating an ODR between 855 and , we detected a statistically significant difference in arterial and venous ODRs, which is consistent with the known difference between arterial and venous blood oxygen saturation, a difference known to exist a priori. Future refinements include spatial localization of measurements to relevant tissues (e.g., intravascular), correcting for the absorption of retinal pigments and water, and utilizing the entire available bandwidth to increase accuracy and precision of hsUHR-OCT oximetry.