JBO Letters

Ultrahigh-speed non-invasive widefield angiography

[+] Author Affiliations
Cedric Blatter, Branislav Grajciar, Tilman Schmoll, Rainer A. Leitgeb

Medical University Vienna, Center of Medical Physics and Biomedical Engineering, Waehringerguertel 18-20, 1090 Vienna, Austria

Thomas Klein, Wolfgang Wieser, Raphael Andre, Robert Huber

Ludwig-Maximilians-Universität München, Lehrstuhl für BioMolekulare Optik, Oettingenstraße 67, 80538 Munich, Germany

J. Biomed. Opt. 17(7), 070505 (Jun 29, 2012). doi:10.1117/1.JBO.17.7.070505
History: Received April 7, 2012; Revised June 6, 2012; Accepted June 11, 2012
Text Size: A A A

Open Access Open Access

Abstract.  Retinal and choroidal vascular imaging is an important diagnostic benefit for ocular diseases such as age-related macular degeneration. The current gold standard for vessel visualization is fluorescence angiography. We present a potential non-invasive alternative to image blood vessels based on functional Fourier domain optical coherence tomography (OCT). For OCT to compete with the field of view and resolution of angiography while maintaining motion artifacts to a minimum, ultrahigh-speed imaging has to be introduced. We employ Fourier domain mode locking swept source technology that offers high quality imaging at an A-scan rate of up to 1.68 MHz. We present retinal angiogram over 48deg acquired in a few seconds in a single recording without the need of image stitching. OCT at 1060 nm allows for high penetration in the choroid and efficient separate characterization of the retinal and choroidal vascularization.

Figures in this Article

Assessment of the retinal and choroidal vascularization yields important diagnostic benefits for major ocular diseases that affect the vascular network already at an early state. The gold standards for their visualization are fluorescein angiography (FA) and indocyanine green angiography. They are commonly used in clinical practice for the diagnosis of vascular occlusions, diabetic retinopathy, and choroidal neovascularization, which is usually a cause of age-related macula degeneration. The invasiveness of these techniques together with undesirable side effects, through the injection of a fluorescent dye, limits the screening capabilities for large populations. Moreover, no depth resolution is available.

Fourier domain optical coherence tomography (FDOCT) is a non-invasive three-dimensional imaging technique with high sensitivity.1 Functional extensions of FDOCT and in particular microcirculation imaging have profited from recent technological developments for high-speed FDOCT.2 Firstly, images are virtually free of motion distortion because of smaller acquisition times. Therefore, the patient comfort is improved by keeping measurement duration to a few seconds. Secondly, the lateral sampling can be increased while keeping the optimal B-scan rate for vascular contrast imaging. If the B-scan rate is low, motion artifacts are more likely to cause unwanted signal decorrelation, reducing the contrast between static tissue and flow. The increase of B-scan rate should, however, be limited so that high sensitivity for small capillary flow is preserved. Recent demonstrations applied B–scan rates of several hundred Hz.3 Last but not least, a large field of view (FOV) can be covered in a single recording. High-speed FDOCT is thus optimal for non-invasive imaging of retinal microvasculature comprehensively in vivo.

The visualization of the microvasculature yields an easily accessible and intuitive way to assess its integrity. During the last years, a number of strategies for contrasting microvasculature based on FDOCT have been introduced. Contrast between flow and static tissue can be obtained by determining changes either between A-scans or successive tomograms in phase, a method known as phase variance,4 or in amplitude, known as speckle variance,5 or by filtering the complex spectrum, known as OMAG.6 Approaches based on phase difference between two spatially separated beams on the retina yielded good contrast of the retinal and choroidal vessel beds over a large range of velocities.7 Other authors introduced a pure intensity-based method to filter out capillaries from high resolution data thus avoiding axial decorrelation artifacts below vessels usually present in variance-based contrasting techniques.8 Generally, the demonstration of these techniques was restricted to small FOV because of limited acquisition speed. It was partially solved by stitching small volumes together.3 A critical point, however, concerning the clinical acceptance of this technique is certainly the associated long total measurement time, because of fixation change and the recording of redundant overlap areas required for registration. Ultrahigh-speed FDOCT is therefore a promising candidate to compete with the FOV and resolution of FA, since a large patch can be covered by a single recording in a few seconds. A rough estimation of the required imaging speed, taking into account a spatial sampling of 2000×1000pixels, 5-fold redundancy for variation calculation, and a total acquisition time shorter than 10 s, shows that an A-scan rate over 1 MHz is required. Recent developments in this imaging speed range are dominated by swept source (SS) systems. Over the last years, multi-MHz imaging was introduced based on the concept of Fourier domain mode locked (FDML) swept sources. Retinal and choroidal imaging with this technology at ultrahigh-speed was demonstrated at a center wavelength of 1060 nm.9 Posterior segment OCT imaging in that water window has the advantage of providing increased penetration into the choroid compared to the common 850 nm region because of reduced scattering. It allows for a better assessment of choroidal vasculature that is particularly important for ocular diagnosis; its network being the main oxygen and nourishment supplier of the retina.

Angiography using intensity variance analysis is especially suitable for high-speed SS-based systems, because it is independent from any acquisition timing problems, such as trigger jitter.

In the present work we therefore chose to combine ultrahigh-speed FDOCT with intensity variance contrasting for non-invasive depth-resolved angiography of the retina and choroid.

The schematic of the optical setup is the same as previously described for MHz imaging.9 The light source is a FDML laser with its Fabry-Perot filter driven at 419kHz. The so-called buffering technique of time-multiplexing is later used to increase the sweep rate by a factor of 4 leading to 1.68 MHz A-scan rate. The spectrum is centered at 1060 nm with a 72 nm sweep range. It produces an 14μm axial resolution in air. Shot noise limited sensitivity of 91dB with 1.7 mW power at the cornea is achieved. Widefield imaging is obtained via a sequence of several commercially available, 2 in. diameter planoconvex lenses (Thorlabs), arranged into 2 groups with effective focal lengths of 78 mm and 41 mm. A 27-μm spot size was produced on the retina. The interference signal was measured with a dual-balanced detector (DBD, WieserLabs, WL-BPD1GA) and digitalized with an 8-bit card at 1.5GSamples/s (Signatec, PX1500-4).

Microcirculation is contrasted against static tissue by calculating first the squared intensity difference between successive tomograms: Display Formula

P(x,yi,z)=[I(x,yi,z)I(x,yi+1,z)]2,
where I(x,y,z)=20log[|FFT(x,y,k)|], and x, y, z are the spatial pixel coordinates corresponding to fast and slow scanning, and the depth coordinate, respectively. The intensity difference tomograms P exhibiting high decorrelation are rejected by comparing their respective mean value to a manually set threshold T. The number of pictures averaged at each position can be formally written as M(y)=i=0N1{[Σx,zP(x,y,z)]<T},with N as the number of B-scans measured at the same position y. The value of T is chosen such as to obtain optimal vessel contrast and to suppress residual bulk motion signatures. The logic operation in the brackets yields 1 or 0 for TRUE or FALSE, respectively. Finally, the motion contrast volume V is obtained by averaging only over the remaining M intensity squared difference tomograms P: Display Formula
V(x,y,z)=1M(y)Σi=0N1{[Σx,zP(x,yi,z)]<T}P(x,yi,z).

The slow axis scanner is driven by a multiple-step function. It permits measuring successive B-scans at the same vertical position y, giving an almost perfect correlation for static tissue. We measure N=5B-scans at 800 vertical positions. The fast axis driving function is a ramp of 70% duty cycle and a frequency of 560 Hz. Each B-scan consists of 2060 A-scans. This leads to an effective volume sampling of 2060×800×460 (xyz). The scan amplitude is set to produce an 48deg FOV on the retina. The total acquisition time for the full FOV is only 7 s.

Figure 1(a) shows the en-face mean projection of the intensity data set of the retina of a healthy volunteer acquired over a wide FOV of 48deg. Large retinal and choroidal vessels are already visible, however, smaller vessels lack contrast. They are, on the other hand, well appreciated in the high contrast en-face mean projection of the calculated 3-D intensity variance set [Fig. 1(c)]. The FOV of our label-free and non-invasive widefield angiography can be well compared to that of standard FA. Figure 1(b) shows a standard FA of a patient for qualitative comparison. An important advantage of our technique as compared to previous non-invasive methods based on OCT is the small acquisition time for such image obtained in a single recording.

Grahic Jump LocationF1 :

(a) Pseudo-SLO fundus obtained by en-face mean projection of the intensity data set. ONH: optic nerve head. Black arrow: low signal region. (b) Representative standard fluorescein angiogram of a different subject for comparison (c) ∼48 deg widefield angiogram, en-face mean projection of the intensity difference 3-D data set calculated from (a). (d) Five-fold averaged tomogram after flattening to the retinal pigment epithelium (RPE). Segmentation borders for Fig. 2 are the inner limiting membrane (ILM) and the RPE.

The depth resolution of OCT allows for differentiating the retinal and choroidal vasculature network for further investigation. For this task, the intensity tomograms were first flattened by detection of the retinal pigment epithelium (RPE) layer [Fig. 1(d)]. In a second step, the position of the inner limiting membrane (ILM) was determined. The corresponding coordinates were used for segmented en-face projection of the intensity difference 3-D set. The first segment consists of the region from ILM to the RPE layer [red bar in Fig. 1(d)], the second segment comprises the vascular structures down to about 50 µm below the RPE (blue bar). Large choroidal vessels are finally segmented in a third layer (green bar). For the retinal layer, the intensity difference 3-D data set was multiplied by a manually thresholded copy of the intensity data set. It reduces background noise, so that high back-scattered blood vessels are better visible. The segmentation of the RPE is not possible in the optic nerve head (ONH) region, hence the ONH vasculature is attributed to the upper layer.

Figure 2 shows a color-coded combined representation of the segmented en-face mean projections for the retina and choroid according to Fig. 1(d). It is possible to resolve, simultaneously, retinal and choroidal vessels. The high lateral sampling allows revealing also small vessels and fine network as shown in Fig. 2(c). The choroid exhibits different kinds of vessels along the depth. Close below the RPE, a dense network of small choroidal vessels is visible [Fig. 2(c)]. In deeper regions, the vessel size increases [Fig. 2(b)]. The low contrast region in the picture’s lower part is the result of low signal in the original intensity tomograms [cf. Fig. 1(a), black arrow].

Grahic Jump LocationF2 :

Color-coded en-face mean projection of the retinal and choroidal vasculature. (a) An 48deg depth-resolved angiogram. (b) and (c) zoom showing large choroidal vessels and fine vasculature.

Small retinal capillaries, particularly in the parafoveal region are only weakly visible. This is due to the coarse sampling. Nevertheless, the current protocol would provide a first widefield overview that can then be used as guide to focus on a suspicious area or a region of interest. Figure 3 shows an 12deg FOV centered on the fovea acquired with the same parameters as for the large FOV. The lateral resolution is obviously sufficient to contrast small parafoveal capillaries efficiently [Fig. 3(b) and 3(c)] in the en-face mean projection of the intensity difference data set, while absent of the intensity projection [Fig. 3(a)]. Furthermore, choroidal vessels are well appreciated [Fig. 3(d)]. It is, however, difficult to resolve single choroidal vessels in the fovea region due to the dense vascular network directly underneath the RPE. The latter also shades the signatures of deeper choroidal vessels in the en-face view. Faster B-scan rates could decrease the difference signal of the choriocapillary layer and enhance the contrast of larger vessels underneath. This has recently been demonstrated employing a dual-beam technique.7

Grahic Jump LocationF3 :

An 12deg FOV centered at the fovea. (a) Pseudo-SLO fundus obtained by en-face mean projection of the intensity data set. (b) Color-coded en-face mean projection of the retinal and choroidal vasculature. (c) Retinal vasculature. (d) Choroidal vasculature.

In conclusion, single recording with short acquisition time is an important step towards large field label-free and non-invasive angiography. Compared to previous work in vascular imaging using volume stitching, the introduction of ultrahigh-speed allows for small acquisition time, as well as short total measurement duration. It is promising, therefore, as a viable alternative to retinal FA. Still, the latter remains advantageous for the detection of blood leakage, as well as for studying the time evolution of vessel filling.

Acknowledgments

Funding from the EU FP7 program is acknowledged (FUN OCT, Grant No. 201880). BG is funded by the Austrian Christian Doppler Association.

Leitgeb  R., Hitzenberger  C., Fercher  A., “Performance of Fourier domain vs. time domain optical coherence tomography,” Opt. Express. 11, (8 ), 889 –894 (2003). 1094-4087 CrossRef
Schmoll  T., Kolbitsch  C., Leitgeb  R., “Ultra-high-speed volumetric tomography of human retinal blood flow,” Opt. Express. 17, (5 ), 4166 –4176 (2009). 1094-4087 CrossRef
Kim  D. et al., “In vivo volumetric imaging of human retinal circulation with phase-variance optical coherence tomography,” Biomed. Opt. Express. 2, (6 ), 1504 –1513 (2011). 2156-7085 CrossRef
Fingler  J. et al., “Volumetric microvascular imaging of human retina using optical coherence tomography with a novel motion contrast technique,” Opt. Express. 17, (24 ), 22190 –22200 (2009). 1094-4087 CrossRef
Mariampillai  A. et al., “Speckle variance detection of microvasculature using swept-source optical coherence tomography,” Opt. Lett.. 33, (13 ), 1530 –1532 (2008). 0146-9592 CrossRef
Wang  R. et al., “Depth-resolved imaging of capillary networks in retina and choroid using ultrahigh sensitive optical microangiography,” Opt. Lett.. 35, (9 ), 1467 –1469 (2010). 0146-9592 CrossRef
Jaillon  F., Makita  S., Yasuno  Y., “Variable velocity range imaging of the choroid with dual-beam optical coherence angiography,” Opt. Express. 20, (1 ), 385 –396 (2012). 1094-4087 CrossRef
Schmoll  T. et al., “Imaging of the parafoveal capillary network and its integrity analysis using fractal dimension,” Biomed. Opt. Express. 2, (5 ), 1159 –1168 (2011). 2156-7085 CrossRef
Klein  T. et al., “Megahertz OCT for ultrawide-field retinal imaging with a 1050 nm Fourier domain mode-locked laser,” Opt. Express. 19, (4 ), 3044 –3062 (2011). 1094-4087 CrossRef
© 2012 Society of Photo-Optical Instrumentation Engineers

Citation

Cedric Blatter ; Thomas Klein ; Branislav Grajciar ; Tilman Schmoll ; Wolfgang Wieser, et al.
"Ultrahigh-speed non-invasive widefield angiography", J. Biomed. Opt. 17(7), 070505 (Jun 29, 2012). ; http://dx.doi.org/10.1117/1.JBO.17.7.070505


Figures

Grahic Jump LocationF1 :

(a) Pseudo-SLO fundus obtained by en-face mean projection of the intensity data set. ONH: optic nerve head. Black arrow: low signal region. (b) Representative standard fluorescein angiogram of a different subject for comparison (c) ∼48 deg widefield angiogram, en-face mean projection of the intensity difference 3-D data set calculated from (a). (d) Five-fold averaged tomogram after flattening to the retinal pigment epithelium (RPE). Segmentation borders for Fig. 2 are the inner limiting membrane (ILM) and the RPE.

Grahic Jump LocationF2 :

Color-coded en-face mean projection of the retinal and choroidal vasculature. (a) An 48deg depth-resolved angiogram. (b) and (c) zoom showing large choroidal vessels and fine vasculature.

Grahic Jump LocationF3 :

An 12deg FOV centered at the fovea. (a) Pseudo-SLO fundus obtained by en-face mean projection of the intensity data set. (b) Color-coded en-face mean projection of the retinal and choroidal vasculature. (c) Retinal vasculature. (d) Choroidal vasculature.

Tables

References

Leitgeb  R., Hitzenberger  C., Fercher  A., “Performance of Fourier domain vs. time domain optical coherence tomography,” Opt. Express. 11, (8 ), 889 –894 (2003). 1094-4087 CrossRef
Schmoll  T., Kolbitsch  C., Leitgeb  R., “Ultra-high-speed volumetric tomography of human retinal blood flow,” Opt. Express. 17, (5 ), 4166 –4176 (2009). 1094-4087 CrossRef
Kim  D. et al., “In vivo volumetric imaging of human retinal circulation with phase-variance optical coherence tomography,” Biomed. Opt. Express. 2, (6 ), 1504 –1513 (2011). 2156-7085 CrossRef
Fingler  J. et al., “Volumetric microvascular imaging of human retina using optical coherence tomography with a novel motion contrast technique,” Opt. Express. 17, (24 ), 22190 –22200 (2009). 1094-4087 CrossRef
Mariampillai  A. et al., “Speckle variance detection of microvasculature using swept-source optical coherence tomography,” Opt. Lett.. 33, (13 ), 1530 –1532 (2008). 0146-9592 CrossRef
Wang  R. et al., “Depth-resolved imaging of capillary networks in retina and choroid using ultrahigh sensitive optical microangiography,” Opt. Lett.. 35, (9 ), 1467 –1469 (2010). 0146-9592 CrossRef
Jaillon  F., Makita  S., Yasuno  Y., “Variable velocity range imaging of the choroid with dual-beam optical coherence angiography,” Opt. Express. 20, (1 ), 385 –396 (2012). 1094-4087 CrossRef
Schmoll  T. et al., “Imaging of the parafoveal capillary network and its integrity analysis using fractal dimension,” Biomed. Opt. Express. 2, (5 ), 1159 –1168 (2011). 2156-7085 CrossRef
Klein  T. et al., “Megahertz OCT for ultrawide-field retinal imaging with a 1050 nm Fourier domain mode-locked laser,” Opt. Express. 19, (4 ), 3044 –3062 (2011). 1094-4087 CrossRef

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging & repositioning the boxes below.

Related Book Chapters

Topic Collections

PubMed Articles
Bilateral perifoveal macular ischemia in sarcoidosis. Retin Cases Brief Rep 2014;8(3):212-4.
Advertisement
  • Don't have an account?
  • Subscribe to the SPIE Digital Library
  • Create a FREE account to sign up for Digital Library content alerts and gain access to institutional subscriptions remotely.
Access This Article
Sign in or Create a personal account to Buy this article ($20 for members, $25 for non-members).
Access This Proceeding
Sign in or Create a personal account to Buy this article ($15 for members, $18 for non-members).
Access This Chapter

Access to SPIE eBooks is limited to subscribing institutions and is not available as part of a personal subscription. Print or electronic versions of individual SPIE books may be purchased via SPIE.org.