High-resolution wide-field imaging of retinal and choroidal blood perfusion with optical microangiography

Lin An; Hrebesh M. Subhush; David J. Wilson; Ruikang K. Wang

doi:10.1117/1.3369811

1 March 2010 High-resolution wide-field imaging of retinal and choroidal blood perfusion with optical microangiography

Lin An, Hrebesh M. Subhush, David J. Wilson, Ruikang K. Wang

Author Affiliations +

J. of Biomedical Optics, 15(2), 026011 (2010). https://doi.org/10.1117/1.3369811

Abstract

We present high-resolution wide-field imaging of retinal and choroidal blood perfusion with optical microangiography (OMAG) technology. Based on spatial frequency analysis, OMAG is capable of visualizing the vascular perfusion map down to capillary-level resolution. An OMAG system operating at 840 nm is used with an A-scan rate of 27,000 Hz, axial resolution of 8 µm, and sensitivity of 98 dB. To achieve wide-field imaging, we capture 16 optical coherence tomography (OCT) 3-D datasets in a sequential order, which together provide an area of ~7.4×7.4 mm² at the posterior segment of the human eye. For each of these datasets, the bulk tissue motion artifacts are eliminated by applying a phase compensation method based on histogram estimation of bulk motion phases, while the displacements occurring between adjacent B-frames are compensated for by 2-D cross correlation between two adjacent OMAG flow images. The depth-resolved capability of OMAG imaging also provides volumetric information on the ocular circulations. Finally, we compare the clinical fluorescein angiography and indocyanine green angiography imaging results with the OMAG results of blood perfusion map within the retina and choroid, and show excellent agreement between these modalities.

1. Introduction

The importance of ocular circulation to the function of the eye and in eye diseases has been known for years. There is convincing evidence that vascular pathologies play an important role in the etiology and progression of a number of ophthalmic diseases, including glaucoma, diabetic retinopathy, and age-related macular degeneration (AMD). Clinically, a better visualization of the ocular vasculature could provide important diagnostic information for these diseases. Currently, the most widely used diagnostic tools for imaging ocular blood vessels are fluorescein angiography (FA) and indocyanine green angiography (ICGA). Although these two methods can provide visualization of retinal (using FA) and choroidal (using ICGA) blood vessels, they require injection of a contrast agent (sodium fluorescein dye and indocyanine green dye) into the blood stream. Further, sodium fluorescein may cause nausea, vomiting, itching, dyspnea, and additional adverse allergic reactions,¹ while indocyanine green dye must be used cautiously in patients who are allergic to iodine and shellfish.² Another main disadvantage is that neither of these two methods can provide depth information on ocular blood vessels.

Optical coherence tomography (OCT)^{3, 4} has become an important tool in ophthalmology for the diagnosis, management, and high-resolution study of the anterior and posterior segments of the eye. Based on low-coherence interferometry, this noninvasive cross sectional imaging technique offers resolution in the micron range, both latterly and axially, with an ability to image to a depth of several millimeters depending on the opacity of the specimen. Combining OCT with the Doppler principle, phase-resolved Doppler optical coherence tomography (PRDOCT) has successfully been demonstrated for in-vivo imaging of blood vessels and the extraction of functional flow information from blood vessels buried within tissue,^{5, 6} particularly after the advent of Fourier domain (FD) OCT.^{7, 8, 9, 10, 11} PRDOCT is based on the evaluation of phase difference between adjacent A-lines. However, the accuracy of PRDOCT is deteriorated by two factors: system sensitivity¹² and the heterogeneous texture pattern of tissue.¹³ Due to the degradation of depth-dependent sensitivity, it is hard for PRDOCT to provide blood vessel networks in the choroid, where the signal-to-noise ratio is relatively low compared to the retina. Moreover, the heterogeneous texture pattern artifact reduces the image quality of PRDOCT even in the retina layers, causing difficulty in providing capillary level information. To overcome these problems, several methods have been proposed. Szkulmowski reported a functional Doppler OCT method called joint spectral and time domain OCT.¹² This technology relies on analysis of the amplitude of the OCT signal to extract the flow information buried within the tissue, and is supposed to have higher sensitivity at low signal-to-noise ratio regions than the conventional phase resolved method. A problem with this method is that it needs high A-line density along one B-scan, limiting its suitability for acquiring ocular blood vessel networks for in-vivo applications. Wang and Ma¹³ proposed that by employing a reverse scan pattern, the heterogeneous texture pattern of the tissue can be successfully suppressed. This method showed a great potential to increase the velocity sensitivity of the PRDOCT method to capillary level. However, an additional reverse scan required at the same location restricts its applications for in-vivo experiments. Besides these functional Doppler methods, Yasuno ¹⁴ demonstrated a segmentation method, the scattering optical coherence angiography method (S-OCA), which is based on the analysis of the low intensity region in the OCT structural image of the choroid of the human eye. Though this method can image blood vessel networks in the human choroid, it essentially loses the functional information of blood flow, such as flow velocity and direction. Moreover, the contrast of the S-OCA image is provided by the scattering OCT signal, which makes it difficult to distinguish the flow signals from the structure signals.

Recently, based on the principle of full range complex FDOCT^{15, 16, 17} with constant carrier frequency modulation along the B-scan direction of an interferogram, Wang and An introduced a new label-free optical microangiography (OMAG) technique to separate the static and moving signals emerged from an illuminated tissue sample. This method was successfully demonstrated for in-vivo imaging of blood vessels at capillary-level resolution.^{18, 19} By introducing a phase compensation method to compensate the movement artifact caused by head or eye movements, the OMAG method was also applied successfully for in-vivo imaging of ocular blood vessel networks both in human retina and choroid.²⁰ Although OMAG can deliver superior performance for depth-resolved imaging in both retina and choroid, the first version of the OMAG algorithm encountered certain limitations caused by the introduction of constant modulation frequency. First, the introduction of modulation frequency by the system hardware is difficult to readjust for different experimental situations, such as the density of the A-line scan in one B-scan and the flow direction. Another issue is that this modulation frequency will be coupled into the Doppler frequency shift, which is normally used to evaluate the velocity of the flow, leading to difficulty in retrieving the information regarding the true velocity.¹⁸ To solve the first issue, Wang presented a digital modulation frequency method for OMAG technology.²¹ But this digitally introduced modulation frequency demands computation power to obtain velocity information of OMAG flow images. To solve these two issues, a modified Hilbert transform method was introduced into OMAG.^{22, 23} By eliminating the modulation frequency, the speed and flexibility for imaging blood flow can be greatly increased.^{23, 24} We have extended that work to perform functional imaging of the choroidal circulation and to obtain wide-field imaging of the ocular vasculature, which is important for successful use in the clinical setting.

In this work, we used OMAG with the modified Hilbert transform method to obtain ocular blood vessel network images. We carried out wide-field imaging of ocular perfusion of both the retina and choroid by reconstructing multiple sequentially registered OMAG projection maps of small areas $(2 \times 2 mm)$ . For the accurate reconstruction of wide-field OMAG imaging, there are two main issues that need to be resolved in the process of sequential capturing of multiple OMAG images. First, for each of these sequential OMAG datasets, the bulk motion artifacts in each B-scan image need to be suppressed. We solved this issue by directly implementing the phase compensation method.²⁰ In this process, rather than using the phase averaging approach reported in Ref. 20, we instead used the histogram method to estimate the bulk motion phase. The second issue is the displacement between adjacent B-scans. To solve this problem, we used the cross correlation between adjacent 2-D OMAG flow images to compensate for the displacements. After resolving these artifacts, we presented wide-field images of retinal and choroidal blood flow maps, and showed the superior performance of OMAG technology to obtain depth-resolved volumetric ocular perfusion of both retina and choroid. Finally, we compared the wide-field OMAG perfusion map with the clinical FA and ICGA images from the same individual. The well-correlated results show the great potential of OMAG in future clinical applications.

2. System Setup

The schematic of the system used to achieve the wide-field ocular blood perfusion image is illustrated in Fig. 1 . The light source is a superluminecent diode with a spectral bandwidth of $46 nm$ centered at $842 nm$ that yields an axial resolution of $\sim 8 μ m$ in air. The light emitted from the light source was coupled into a Michelson interferometer, which divides the incident beam into the reference arm and sample arm. The sample light was coupled into a custom-made imaging probe with a pair of galvometer scanners to achieve $x - y$ scanning and to collect the light backscattered from the posterior part of the eye. The optical power of the beam on the cornea was $\sim 700 μ W$ , which is lower than the ANSI exposure limit. Together with the human eye lens, the imaging probe includes a collimating lens, objective lens, and an ocular lens that provided a lateral resolution of $\sim 16 μ m$ . In the reference arm, to compensate the dispersion, light was delivered into the stationary mirror through an $\sim 2 - cm$ water cell and a set of optics that is matched with the sample arm. Unlike our previous work, we do not need to introduce a constant modulation frequency into one OCT B-scan. During imaging, the reference arm mirror was kept stationary. The interference signal between the reference arm and sample arm was then collected and delivered into a home-built high-speed spectrometer. This spectrometer had a spectral resolution of $0.055 nm$ that achieved a maximum imaging depth of $\sim 3 mm$ on each side of the zero delay line with a line camera of $2048 pixels$ .

Fig. 1

Schematic of OMAG system used to achieve wide field of view fundus blood vessel map: PC, polarization controller; and CCD, the charged couple device. The aiming source is a blue light source used for minimizing eye movement during the experiment progress. (Color online only.)

For maintaining high-resolution imaging of blood vessels, the area of each OCT scan was set relatively small. In the lateral direction ( $x$ direction), the $x$ -galvoscanner was driven by a $27 - Hz$ sawtooth waveform with $2 - mm$ amplitude, i.e., the imaging speed at 27 frames per second. For elevational direction ( $y$ direction), the $y$ -galvoscanner was driven by a $0.27 - Hz$ sawtooth waveform with $2 - mm$ amplitude. In this study, the camera integrating time was set to be $36 μ s$ , allowing $1 μ s$ for transferring the spectral data from the CCD camera ( $2048 pixels$ , A-scan) to the host computer via Cameralink^™ and a high-speed frame grabber board. For each B-scan we captured 1000 A-scans, and for C-scan (elevational direction) we captured 100 B-scans. Hence, one subarea image dataset contained 2048 (A-scan) by 1000 (B-scan) by 100 (C-scan) pixels.

To achieve wide-field imaging of fundus vasculature $(7.4 \times 7.4 mm)$ , we sequentially scanned different areas of the posterior part of the human eye by a scanning order illustrated in Fig. 2 . In Fig. 2, the background is an OCT fundus image obtained through the technique described in Ref. 25. To obtain this background image, the scan area of our OCT system was set at $\sim 9 \times 9 mm$ . The small blue squares represent the scan area of each smaller OCT dataset with the number indicated referring to the scanning order in the whole experiment. Between every two subarea scans, we introduced $\sim 0.2 - mm$ overlap area to make sure no blank area occurred when reconstructing the wide-field images. For every subscan, it took around $3.7 s$ to capture the dataset. Between two subscans, the volunteer was instructed to close his eye for $15 s$ for a rest. Hence, it took about $5 \min$ to acquire the wide-field view of the fundus blood vessel maps. During the whole experiment, the volunteer fixated on a blue light spot emitted by the aiming source to minimize eye movement. The operating software for the whole system, including beam scanning, data acquisition, data storage, and hand-shaking between them, was developed on a Visual C++ platform.

Fig. 2

Scanning protocol to achieve wide field of view fundus blood vessel image. Numbers on the subareas refer to the order they are acquired. Background image is the OCT fundus image of one OCT volume $(\sim 9 \times 9 mm)$ .

3. Method

3.1.

Modified Hilbert Transform Optical Microangiography Flow Image

Based on the modified Hilbert transform, OMAG can separate the flow signals from the static tissue signals without introducing any modulation frequency. The detailed description of this method has been reported previously.^{23, 24} Here, we briefly review the principle of this method.

The B-scan interference signal captured by the OCT system can be written as:²³

1

I (t, k_{j}) = S (t, k_{j}) {E_{R}^{2} + 2 E_{R} \int_{- \infty}^{\infty} a (z, t) \cos [2 k_{j} n (t) z] d z + 2 E_{R} a (z_{l}) \cos [2 k_{j} n (t) (z_{l} - v t)]},

where

j = 1, 2, \dots, 1024

is the pixel number of the line scan CCD camera,

k_{j}

is the wavenumber of the light captured by the

j

’th pixel,

I (k_{j})

is the light intensity detected by the

j

’th detector,

t

is the time variable when one A-line

I (k)

was captured during one B-scan,

S (k_{j})

is the spectral density of the light source at

k_{j}

,

n

is the refractive index,

a (z)

is the amplitude of the backscattered light and

z

is the depth where the light backscattered from, and

v

is the velocity of moving particles, for example blood cells, which are buried in the tissue at depth

z_{1}

. In Eq. 1, we ignore the self cross-correlation term occurring between the light backscattered within the sample, because it is relatively weak compared to that reflected from the reference mirror.

In Eq. 1, the first term is the dc signal produced by the light reflected from the reference mirror. The second term is the spatial frequency component of the static tissue sample, which we call the heterogeneous frequency component, and which can provide the static structural information of the sample. Because of the optically heterogeneous nature of the biological tissue sample, this frequency is randomly distributed around zero frequency with a bandwidth of $BW$ .¹⁶ The bandwidth of this heterogeneous frequency is determined by the optical characteristics of the tissue sample and the density of A-line scans in a B-scan image.²⁶ The third term is the Doppler beating signal, which is introduced by the moving particles in the tissue sample. Depending on the magnitude of the moving particle velocity, the Doppler beating signal shifts away from the heterogeneous bandwidth $BW$ . If we first apply high-pass filtering with cutting frequency $BW ∕ 2$ and then apply the Hilbert transformation along the B-scan direction, we can extract the negative $(v < 0)$ , Eq. 2, and positive $(v > 0)$ , Eq. 3, directional flows as follows:¹⁸

2

{\tilde{I}}_{N} (t, k_{j}) = 2 S (k_{j}) E_{R} a (z_{l}) {\cos [2 k_{j} n (t) (z_{l} - v t)] - i \sin [2 k_{j} n (t) (z_{l} - v t)]},

3

{\tilde{I}}_{P} (t, k_{j}) = 2 S (k_{j}) E_{R} a (z_{l}) {\cos [2 k_{j} n (t) (z_{l} - v t)] + i \sin [2 k_{j} n (t) (z_{l} - v t)]},

where

i = \sqrt{- 1}

,

{\tilde{I}}_{N} (k_{j}, t)

, and

{\tilde{I}}_{P} (k_{j}, t)

represent the negative and positive flow signals, respectively. By subsequently taking the Fourier transform along the

k

direction, the depth-resolved 2-D OMAG flow map of functional vessels within the microcirculation tissue beds can be retrieved.

3.2.

Phase Compensation for Two-Dimensional Optical Microangiography Flow Image

The movement artifact is one of the most distracting issues for in-vivo imaging of the posterior part of the human eye. This artifact is caused by the involuntary movement of head or eye and can seriously deteriorate the image quality. In our previous work,²³ we proposed a phase compensation method to effectively resolve this issue in the traditional OMAG system. In our modified Hilbert transform-based OMAG system, we directly implemented this phase compensation method to minimize the movement artifact issues.

When the sample moves during imaging, the velocity of this movement $v_{s}$ will introduce a Doppler frequency shift into the light backscattered from the sample with respect to the reference light. If we analyze the frequency components in one B-scan along the $x$ direction, then this Doppler frequency due to the bulk motion will also represent a Doppler beating frequency shift. When this frequency shift is large enough $(> BW ∕ 2)$ , the bulk motion signal can be detected through the OMAG method that inevitably gives a motion artifact in the final OMAG flow image. Figure 3 shows typical imaging results obtained from the OMAG method. Figure 3 presents a typical cross sectional structure image of the posterior part of a human eye from a volunteer in our laboratory. The structure information of different layers, such as retina, choroid, and sclera, are clearly delineated from this image. Figure 3 is the corresponding cross sectional flow image, directly obtained from the OMAG method without applying the phase compensation method. Though two large blood vessels are highlighted by the OMAG method, the other part of the image is corrupted by the motion artifacts (white arrow). In this case, the detected moving signal can be written as:

4

I (t, k_{j}) = 2 S (k_{j}) E_{R} a (z_{l}) \cos [2 k_{j} n (t) (z_{l} - (v_{P} + v_{S}) t)],

where

v_{p}

is the velocity of the moving particles in the specimen. Equation 4 can be further modified in terms of phase

ϕ_{M}

, as shown below:

5

I (t, k_{j}) = 2 S (k_{j}) E_{R} a (z_{l}) \cos {[2 k_{j} n (t) [z_{l} - ϕ_{M} (t, z))]},

where

ϕ_{M} = ϕ_{S} + ϕ_{P}

is the total phase caused by the movement, i.e., the sum of phases caused by the moving particles

ϕ_{P}

and the phase caused by the sample movements

ϕ_{S}

. To eliminate the bulk motion artifacts, we need to compensate the sample bulk motion phase

ϕ_{S}

before applying the OMAG algorithm.

Fig. 3

OMAG image results with and without phase compensation. (a) Conventional cross section OCT structural image. (b) 2-D OMAG flow image without phase compensation. (c) 2-D OMAG flow image with average bulk motion phase estimation for phase compensation.

As we know that the OMAG image has its origin in FDOCT technology, the phase term $ϕ_{M}$ may be obtained from the FDOCT complex signal ${\tilde{B}}_{OCT} (t, z)$ by applying the fast Fourier transform (FFT) along wave number $(k)$ :

6

{\tilde{B}}_{OCT} (t, z) = A (t, z) \exp [i ϕ_{M} (t, z)],

where

A (t, z)

is the amplitude of the OCT complex signal, from which the structural OCT cross sectional image can be reconstructed. Because the time variable

t

during one OCT B-scan corresponds uniquely to the lateral

x

position as well as the index of A-line numbers

(i = 1, 2, 3 \dots 1000)

, the phase term

ϕ_{M} (t, z)

can also be represented as

ϕ_{M}^{i} (x, z)

. By the PRDOCT method, the adjusted phase difference between two A-scan lines

Δ ϕ_{M}^{i}

could be directly extracted from the complex OCT signal, which can be expressed as:

7

Δ ϕ_{M}^{i} = Δ ϕ_{P}^{i} + Δ ϕ_{S}^{i},

where

Δ ϕ_{P}^{i}

and

Δ ϕ_{S}^{i}

are the phase differences between adjacent A-scan lines caused by the moving particles and sample bulk motions, respectively. If

Δ ϕ_{S}^{i}

information can be retrieved separately, the phase estimation of the bulk motion at every A-line position may be obtained through the following recursive function:

8

ϕ_{S}^{1} = 0, ϕ_{S}^{i} = Δ ϕ_{S}^{i} + ϕ_{S}^{i - j} i = 2, 3, 4 \dots 1000 .

Then we use the bulk motion phase term obtained through Eq. 8 to compensate the phase term of Eq. 6, so that the phase change induced by the sample bulk motion is eliminated and the corrected OCT complex signal is obtained:

9

{\tilde{B}}_{C} (t, z) = A (t, z) \exp {i [ϕ_{M}^{i} (x, z) - ϕ_{S}^{i} (x, z)]} .

Through an inverse FFT along the

k

-space direction, we can obtain the corrected FDOCT interferogram as follows:

10

b_{C} (x, k) = {FFT}^{- 1} A ((t, z) \exp {i [ϕ_{M}^{i} (x, z) - ϕ_{S}^{i} (x, z)]}) .

After applying the OMAG algorithm to this corrected interferogram, the flow image without any bulk motion artifacts can be finally reconstructed.

As we discussed before, a key step for the phase compensation method is how to estimate the phases due to bulk tissue motion between adjacent A—lines $Δ ϕ_{S}^{i}$ . In our previous work,²³ we estimated $Δ ϕ_{S}^{i}$ by averaging the phase differences between adjacent A-lines along the $z$ direction within the tissue-signal region of an OCT A-scan signal. This method was based on the assumption that the blood vessels within the cross section of the posterior eye segment are quite small compared to the static tissue region itself. This assumption is true for most parts of the posterior segment, except the optic disk, where the blood vessels are relatively large and in which the blood flow is also very fast. However, the phase changes induced by the flow in large blood vessels within the optic disk are not negligible when using the phase-averaging approach,²³ which instead could impose an additional phase offset to the bulk motion estimation. This additional phase offset will lead to a correction artifact in the final output flow image. Such an artifact is illustrated in Fig. 3, resulting from the phase compensation based on the phase averaging method. As we can see, in the left part where there are no large blood vessels, the bulk motion artifacts are eliminated quite successfully. However, in the right part, the existence of the large blood vessel (blue arrow) causes a serious correction artifact above it (pointed by the red arrow). To solve this problem, a more accurate method is therefore needed for the phase estimation due to the bulk tissue motion.

We note that the phase changes solely due to the bulk motion are almost the same at all positions within the tissue region along one OCT A-scan. On the other hand, the blood vessels often occupy a small proportion of the tissue region along one A-scan, and also the velocity of blood flow is parabolically distributed across a wide range. These facts determine that the number of possible positions carrying the phase caused by the bulk motion will be much larger than those of the positions with other phase values along one A-scan. For this reason, we adopted a histogram approach to estimate the phases that are caused by the bulk tissue motion. Next, we briefly describe the procedures of this histogram method, which are similar to that developed by Makita ¹¹

After estimating the phase differences $Δ ϕ^{i} (z)$ between the adjacent A-scans within the tissue region by PRDOCT method, the maximum value of phase difference $Δ ϕ_{\max}^{i} (z)$ and the minimum value of phase difference $Δ ϕ_{\min}^{i} (z)$ , as well as the number of phase samples $m$ , were calculated. The bin width $h$ for the histogram was then determined by the Freedman and Diaconis rule:²⁷

11

h = 2 I Q m^{- 1 ∕ 3},

where

I Q

is the interquartile range. The number of bins

N

was then calculated by

N = [Δ ϕ_{\max}^{i} (z) - Δ ϕ_{\min}^{i} (z)] ∕ h

. Based on

h

and

N

, the histogram of the phase difference samples in one A-scan

H [Δ ϕ^{i} (z)]

was obtained. Subsequently, the estimation of phases due to the bulk motion

Δ ϕ_{S}

were determined through the

H [Δ ϕ^{i} (z)]

function by locating the phase range that contains the maximum number of phase values. This process is demonstrated in Fig. 4 . The phase difference map directly obtained from the PRDOCT method is shown in Fig. 4, from which we extracted one segment of A-line phase difference signal (presented by the red line). After obtaining the bin width

h

and bin number

N

of this segment, the resulted histogram is presented in Fig. 4. Usually,

h

can be directly calculated from Eq. 11. However, the

h

so evaluated is often too large to achieve an accurate

Δ ϕ_{S}

. For example, the minimal phase difference

Δ ϕ_{S}

that can be estimated from Fig. 4 is

\sim 0.2 rad

, which corresponds to

360 μ m ∕ s

. One way to solve this problem is to further calculate

M

different histograms with different phase origins

Δ ϕ_{SO}

but with identical bin width and bin number.

12

Δ ϕ_{SO}^{j} = Δ ϕ_{\min}^{i} (z) + (\frac{h}{m}) \times (j - 1) j = 1, 2, 3 \dots M .

By averaging these

M

histograms, more accurate bulk motion estimation can be obtained by locating the phase value, which contains the maximum number of

z

positions along one A-scan. In our experiment, we set

M

to be 8. The result of the averaged histogram is shown in Fig. 4, which gives a bin width of

\sim 0.025 rad

, corresponding to a velocity of

\sim 45 μ m ∕ s

. Following the procedure described before [Eqs. 9, 10], the corrected interference spectrum can be obtained. By applying the OMAG algorithm, we can retrieve the 2-D OMAG flow image with the bulk motion artifacts successfully suppressed, even in the large vessel region. Shown in Fig. 4 is the resulting 2-D OMAG flow image based on the histogram bulk motion phase estimation. By comparing Fig. 4 with Fig. 3, there is clear evidence that the bulk motion artifacts are successfully removed, not only in the left part but also in the right part.

Fig. 4

Phase compensation results based on the histogram method. (a) PRDOCT phase difference map. (b) Histogram of bulk motion phase estimations without averaging. (c) Histogram after averaging. (d) OMAG flow image after phase compensation.

3.3.

Displacement Compensation between Adjacent Optical Coherence Tomography B-Scans

In the case of wide-field imaging, another serious problem caused by the involuntary movement of the human head and eye is the displacement between adjacent B-scans. During imaging, although the volunteer was asked to fixate the blue aiming spot to minimize eye movement, the motion artifacts are still quite dominant and occur occasionally through the whole progress of capturing the 16 OCT datasets, from which the wide-field ocular perfusion map is reconstructed. As previously reported,^{11, 20} the axial displacement can be compensated for with relative ease by the cross correlation method between two A-lines at the same position of two adjacent B-scans. This is because the retina has well-layered structures along the depth direction that produce a sharp peak for the cross correlation function. Thus, this cross correlation is suitable for determining the displacement by locating its peak position. Using the position of peak value, we can determine the lateral displacement between these two B-scans. For example, if the position index of the peak value is zero, there is no displacement between them. If the position index is positive, the second B-scan moved to the right compared to the previous B-scan and vice versa. However, for several reasons, the cross correlation method based on an OCT structural image often does not work effectively well to compensate the lateral displacements for most parts of the posterior part of the human eye. First of all, except for the optic nerve head region, the structures of the retina are very smooth. Particularly when a small area of $\sim 2 \times 2 mm$ is scanned, the variance of structures is very small in the lateral direction. This means that the structural information presented by the same layer will be very similar along the lateral direction. Because of these reasons, the calculated cross correlation function will behave as if having a very broad peak, which leads to difficulty in locating the peak position for estimating the lateral tissue motion. To give an example, Fig. 5 shows a typical FDOCT B-scan structure image. The cross correlation between the same depth position marked by the red line of two adjacent B-scans is shown in Fig. 5. As it can be seen, the cross correlation function has a very broad peak, and thus the subsequent determination of the peak position can be inaccurate. To obtain accurate determination of the peak position, we propose to use the OMAG flow images to evaluate the cross correlation function, rather than using the structural image. This is reasonable because the OMAG flow images possess more high spatial frequency components than the structural images. According to the information theory as applied to the pattern recognition; the more high spatial frequency components, the sharper the correlation peak. Figure 5 gives the OMAG flow images corresponding to Fig. 5. The cross correlation function at the same position as for Fig. 5 is shown in Fig. 5. It is clear that the use of OMAG flow images produces a cross-correlation function with a very sharp peak, thus delivering an accurate estimation of the displacement between adjacent B-scans.

Fig. 5

Cross correlation results of FDOCT structure image and OMAG flow image. (a) Typical FDOCT B-scan structure image. (b) Cross-correlation function between two lines at the same depth position [marked by the red line in (a)] between the adjacent B-scans. (c) Typical OMAG B-scan flow image. (d) Cross correlation function between two lines at the same depth position [marked by the red line in (c)] between the adjacent B-scans. (Color online only.)

Based on the prior analysis, we used the cross correlation between adjacent OMAG flow images to obtain an estimation of axial displacement between adjacent B-scans. We also used the same approach to estimate both the axial and lateral displacements directly by 2-D cross correlation method between two neighboring OMAG frames. Figure 6 shows the comparison between the images without and with motion compensation. Figures 6 and 6 are OMAG fundus vasculature images for one 3-D OCT dataset before and after displacement compensation, respectively. As we can see in Fig. 6, the two blood vessels in this data volume have broken into several parts because of the eye motion, marked by the blue circles. After the motion compensation, these two vessels were reconnected in Fig. 6. This method works quite well for compensating lateral displacement, even in a small area without any structural landmarks. Figure 6 is the elevational cross sectional image before motion compensation, which was extracted from the position marked by the red line in Fig. 6. Figure 6 is the corresponding image after the motion compensation. The position of this image is also marked by the red line in Fig. 6, which is the same position as in Fig. 6. By comparing these two images, we can see that Fig. 6 is smoother than Fig. 6, where a big strike (the white arrows) is successfully compensated to a level similar to other positions. We further discuss the motion compensation results in Sec. 4, based on the 2-D OMAG $x - y$ projection images.

Fig. 6

$x - y$ OMAG fundus flow image obtained from one OCT 3-D dataset (a) before and (b) after lateral displacement compensation. FDOCT elevational cross sectional structural image at the position marked as red line in the upper panel (c) before and (d) after axial displacement between B-scans is compensated. (Color online only.)

3.4.

Vascular Perfusion Maps for Retina and Choroid

An important characteristic of OMAG, which is inherited from FDOCT technology, is the depth-resolved imaging capability. Through the OMAG method, blood vessels located at different depths, for example in the retina and choroid, can be separated relatively easily. To achieve the isolation of blood vessel located in retina and choroid, two masks were created through a segmentation process based on the FDOCT B-scan structure image. The detailed procedure of the creation of the masks for both retina and choroid has already been described in our previous work.²⁰ These two masks were then applied onto the OMAG flow image to separate the blood vessels in retina and choroid. After combining 100 B-scans along the elevational direction (C-scan) and doing the maximum projection, the vascular perfusion maps for both the retina and choroid are finally obtained.

3.5.

Threee-Dimensional Volumetric Image of Vascular Perfusion within Retina and Choroid

After compensating the axial and lateral displacement, we could obtain 3-D vascular images by stacking 100 B-scans together. Through the two masks mentioned before, we could also get the 3-D volumetric image of retinal and choroidal blood vessel networks. The results are presented in the following section.

4. Results

4.1.

Demonstration of Optical Microangiography Image Results and Displacement Compensation Effect

Because of the high sensitivity to the flow, OMAG delivers superior performance for imaging the blood vessel networks both in retina and choroid. Especially after displacement compensation, even the microvascular networks could be effectively presented through this method. In this section, we give a detailed demonstration of the results, obtained from the previous described procedure for in-vivo experimental 3-D datasets captured by the system described in Fig. 1. Each 3-D FDOCT interferogram is composed of $2048 pixels$ in one A-line ( $z$ direction), 1000 A-lines in one B-scan ( $x$ direction) and 100 B-scans along one C-scan ( $y$ direction). This dataset covered $\sim 2 mm$ ( $x$ direction) $\times 2 mm$ ( $y$ direction) area information of the posterior part of the volunteer eye. The data acquisition time for such a dataset was $< 4 \sec$ , with the camera line scan rate at $27 kHz$ .

Figure 7 is a typical FDOCT B-scan structure image. According to the structural image, two masks were created from the procedure described in Sec. 3.4 The first one is the anterior boundary presented by the red line and posterior boundary presented by the yellow line, shown in Fig. 7, which is used for vasculature reconstruction in retina. The second one is the anterior boundary presented by the green line and posterior one presented by the blue line in Fig. 7, which is used for choroidal blood vessel reconstruction. Figure 7 is the flow image obtained by applying the OMAG algorithm on the same cross sectional position of Fig. 7. After mirroring the four boundaries in Fig. 7 onto Fig. 7, we successfully segmented out the blood vessels of retina and choroid, respectively. In the retinal region, between the red and yellow line, not only the large vessels but also the microvessels are demonstrated. The choroidal region, between the green and blue lines, presents an abundance of blood vessels. Combining 100 frames along the elevational direction and performing the maximum projection, the ocular vascular projection image for 3-D scanned tissue volume can be obtained. Such results are shown in Fig. 8 .

Fig. 7

Segmentation of retina and choroid. (a) A typical OCT structure image and (b) the corresponding OMAG flow image. The red, yellow, green, and blue curves correspond to the anterior boundary of retina, posterior boundary of retina, anterior boundary of choroid, and posterior boundary of choroids, respectively. (Color online only.)

Fig. 8

OMAG blood perfusion results for (a) retina and (b) choroid before motion compensation; and corresponding (c) retina and (d) choroid after motion compensation.

Figures 8 and 8 are the ocular perfusion maps of retina and choroid before applying the motion compensation method. Though both of them present abundant blood vessels, lots of them (marked by the ellipse) are not connected well with each other because of eye movements. Figures 8 and 8 are the vascular perfusion maps of retina and choroid after applying the displacement compensation. The broken blood vessels presented in Figs. 8 and 8 due to motion are reconnected again, with similar areas marked by the ellipses. Especially in the region marked by the green ellipse in the retina and choroid, the blood vessel networks are reasonably well presented compared to those in Figs. 8 and 8.

After motion compensation, we could obtain the 3-D volumetric flow image of both retina and choroid by combining 100 B-scan frames in each 3-D subdataset. The results are shown in Fig. 9 . Figures 9 and 9 demonstrate the 3-D volumetric blood vessel networks in the retina and choroid. In Fig. 9, the retina blood vessels (coded with green) and choroidal blood vessels (coded in red) are merged together to present their relative positions. Together with a cross section B-scan structure image, Fig. 9 also demonstrates the positions of blood vessels relative to the structure.

Fig. 9

Volumetric image of blood vessel networks in (a) retina and (b) choroid. (c) Combined volumetric image, together with a cross section structure image.

4.2.

Wide Field of View of Ocular Perfusion Map Both for Retina and Choroid

As we discussed in the first section, the wide field of view ocular perfusion maps of retina and choroid may provide a broad spectrum of information for ophthalmic diagnosis and eye disease studies. To obtain the wide-field imaging of ocular perfusion maps while maintaining the high resolution OMAG imaging of blood perfusion, we captured 16 OCT data cubes in a predesigned order (see Fig. 2). Each of these small datasets covered around a $2 \times 2 - mm$ area. In both lateral and elevational directions, we introduced a $\sim 0.2 - mm$ overlap between adjacent small datasets to eliminate blank space during the whole image stitching process. In this case, all 16 datasets covered an entire area of about $7.4 \times 7.4 mm$ over the posterior segment of the eye. For every OCT raw data, we first obtained all cross sectional structural images, i.e., OCT images, to segment out the tissue signal region from which the retina and choroid layers, through the creation of two masks for retina and choroid, were separated. Subsequently, the motion compensation methods described in Secs. 3.2, 3.3 were applied to remove/minimize the motion artifacts due to the inevitable eye movement during imaging. Next, the maximum projection approach along the $z$ direction was used to obtain the projection images $(x - y)$ of the separate vasculature maps of retinal and choroidal blood vessels. Finally, all the projection images of the 16 areas were manually stitched together to form two wide-field vascular perfusion maps of retina and choroid. The results are illustrated in Figs. 10 and 10 , respectively. Figure 10 is the retina blood perfusion map and Fig. 10 is the choroid blood perfusion map. Both images demonstrate clear, well-connected blood vessel networks, not only the large blood vessels but also the capillary vessels.

Fig. 10

Wide field of view blood perfusion maps for (a) retina and (b) choroids, compared with (c) fluorescein angiography and (d) indocyanine green angiography. (Color online only.)

At last, to compare the performance of OMAG technology with standard clinical ocular perfusion imaging modalities, we performed FA and ICGA procedures for the same volunteer on whom we performed OMAG scans. The results are presented in Figs. 10 and 10 for FA and ICGA, respectively. In Figs. 10 and 10, the light-blue squares correspond to the scan area of OMAG imaging shown in Figs. 10 and 10. By comparison, the OMAG results demonstrated very good agreement with the clinical FA and ICGA results. Through scrutinizing the appearances of the blood vessels, almost all noticeable vessels presented by FA and ICGA images are demonstrated by OMAG imaging. Furthermore, in the retinal perfusion map [Fig. 10], the OMAG results demonstrated more microcirculation information in the fovea region than that of the FA image [marked by the ellipse in Fig. 10]. The main advantages of OMAG compared to FA and ICGA are that the images are obtained without injected contrast agents, and the vessels can be correlated with the depth of location in the retina or choroid. These two advantages of OMAG wide field of view perfusion maps indicate this imaging modality will potentially be a very powerful tool for in-vivo clinical application in the future.

5. Conclusion

We use OMAG technology to achieve wide-field imaging of ocular perfusion maps, both in retina and choroid, by manually stitching 16 small area datasets. The datasets provide high-resolution blood vessel maps down to the capillary level. To achieve this, we develop methods to solve the problems associated with subject movements. We use a histogram-based method to estimate the bulk motion phase by which the bulk motion artifacts are eliminated/minimized through a phase compensation method. To remove the displacement between adjacent frames caused by involuntary head and eye movements, we demonstrate a method that uses the correlation function between two adjacent B-scan OMAG flow frames to estimate the displacement. We show that these motion compensation methods work well to remove/minimize the subject movement artifacts presented in OMAG flow images. Finally, we compare the wide-field OMAG perfusion maps with the currently available FA and ICGA images, and the well-correlated results show the great potential of OMAG in future clinical applications.

Acknowledgments

This study was supported in part by research grants from the National Institutes of Health (R01 HL093140-01, R01 EB009682-01, and R01 DC010201-01), the American Heart Association (0855733G), and an unlimited grant from Research to Prevent Blindness. The content is solely the responsibility of the authors and does not necessarily represent the official views of grant-giving bodies.

References

1.

Kwiterovich K. A., Maguire M. G., Murphy R. P., Schachat A. P., Bressler N. M., Bressler S. B., Fine S. L., “Frequency of adverse systemic reactions after fluorescein angiography. Results of a prospective study,” Ophthalmology, 98 1139 –1142 (1991). 0161-6420 Google Scholar

2.

Hope-Ross M., Yannuzzi L., Gragoudas E., Guyer D., Slakter J., Sorenson J., Krupsky S., Orlock D., Puliafito C., “Adverse reactions due to indocyanine green,” Ophthalmology, 101 529 –533 (1994). 0161-6420 Google Scholar

3.

Fercher A. F., Drexler W., Hitzenberger C. K., Lasser T., “Optical coherence tomography—principles and applications,” Rep. Prog. Phys., 66 239 –303 (2003). 10.1088/0034-4885/66/2/204 0034-4885 Google Scholar

4.

Tomlins P. H., Wang R. K., “Theory, development and applications of optical coherence tomography,” J. Phys. D Appl. Phys., 38 2519 –2535 (2005). 10.1088/0022-3727/38/15/002 Google Scholar

5.

Chen Z. P., Milner T. E., Srinivas S., Wang X., Malekafzali A., van Gemert M. J. C., Nelson J. S., “Noninvasive imaging of in vivo blood flow velocity using optical Doppler tomography,” Opt. Lett., 22 (14), 1119 –1121 (1997). 10.1364/OL.22.001119 0146-9592 Google Scholar

6.

Zhao Y. H., Chen Z. P., Ding Z. H., Ren H., Nelson J. S., “Real-time phase-resolved functional optical coherence tomography by use of optical Hilbert transformation,” Opt. Lett., 27 (2), 98 –100 (2002). 10.1364/OL.27.000098 0146-9592 Google Scholar

7.

Leitgeb R. A., Schmetterer L., Drexler W., Fercher A. F., Zawadzki R. J., Bajraszewski T., “Real-time assessment of retinal blood flow with ultrafast acquisition by color Doppler Fourier domain optical coherence tomography,” Opt. Express, 11 3116 –3121 (2003). 1094-4087 Google Scholar

8.

White B. R., Pierce M. C., Nassif N., Cense B., Park B., Tearney G., Bouma B., Chen T., de Boer J., “In vivo dynamic human retinal blood flow imaging using ultra-high-speed spectral domain optical Doppler tomography,” Opt. Express, 11 3490 –3497 (2003). 1094-4087 Google Scholar

9.

Vakoc B. J., Yun S. H., de Boer J. F., Tearney G. J., Bouma B. E., “Phase-resolved optical frequency domain imaging,” Opt. Express, 13 5483 –5493 (2005). 10.1364/OPEX.13.005483 1094-4087 Google Scholar

10.

Leitgeb R. A., Schmetterer L., Hitzenberger C. K., Fercher A. F., Berisha F., Wojtkowski M., Bajraszewski T., “Real-time measurement of in vitro flow by Fourier-domain color Doppler optical coherence tomography,” Opt. Lett., 29 171 –173 (2004). 10.1364/OL.29.000171 0146-9592 Google Scholar

11.

Makita S., Hong Y., Yatagai M. Y. T., Yasuno Y., “Optical coherence angiography,” Opt. Express, 14 7821 (2006). 10.1364/OE.14.007821 1094-4087 Google Scholar

12.

Szkulmowski M., Szkulmowska A., Bajraszewski T., Kowalczyk A., Wojtkowski M., “Flow velocity estimation using joint spectral and time domain optical coherence tomography,” Opt. Express, 16 6008 –6025 (2008). 10.1364/OE.16.006008 1094-4087 Google Scholar

13.

Wang R. K., Ma Z. H., “Real-time flow imaging by removing texture pattern artifacts in spectral-domain optical Doppler tomography,” Opt. Lett., 31 3001 –3003 (2006). 10.1364/OL.31.003001 0146-9592 Google Scholar

14.

Yasuno Y., Hong Y. J., Makita S., Yamanari M., Akiba M., Miura M., Yatagai T., “In vivo high-contrast imaging of deep posterior eye by

1 - μ m

swept source optical coherence tomography and scattering optical coherence angiography,” Opt. Express, 15 6121 –6139 (2007). 10.1364/OE.15.006121 1094-4087 Google Scholar

15.

Wang R. K., “In vivo full rang complex Fourier domain optical coherence tomography,” Appl. Phys. Lett., 90 054103 (2007). 10.1063/1.2437682 0003-6951 Google Scholar

16.

Wang R. K., “Fourier domain optical coherence tomography achieves full range complex imaging in vivo by introducing a carrier frequency during scanning,” Phys. Med. Biol., 52 5897 –5907 (2007). 10.1088/0031-9155/52/19/011 0031-9155 Google Scholar

17.

An L., Wang R. K., “Use of scanner to modulate spatial interferogram for in vivo full range Fourier domain optical coherence tomography,” Opt. Lett., 32 3423 –3425 (2007). 10.1364/OL.32.003423 0146-9592 Google Scholar

18.

Wang R. K., Jacques S. L., Ma Z. H., Hanson S., Gruber A., “Three dimensional optical angiography,” Opt. Express, 15 4083 (2007). 10.1364/OE.15.004083 1094-4087 Google Scholar

19.

Wang R. K., “Three dimensional optical angiography maps directional blood perfusion deep within microcirculation tissue beds in vivo,” Phys. Med. Biol., 52 N531 –N537 (2007). 10.1088/0031-9155/52/23/N01 0031-9155 Google Scholar

20.

An L., Wang R. K., “In vivo volumetric imaging of vascular perfusion within human retina and choroids with optical micro-angiography,” Opt. Express, 16 (15), 11438 –11452 (2008). 10.1364/OE.16.011438 1094-4087 Google Scholar

21.

Wang R. K., “Directional blood flow imaging in volumetric optical micro-angiography achieved by digital frequency modulation,” Opt. Lett., 33 1878 –1880 (2008). 10.1364/OL.33.001878 0146-9592 Google Scholar

22.

Tao Y. K., Davis A., Izatt J. A., “Single-pass volumetric bidirectional blood flow imaging spectral domain optical coherence tomography using a modified Hilbert transform,” Opt. Express, 16 12350 –12361 (2008). 10.1364/OE.16.012350 1094-4087 Google Scholar

23.

Wang R. K., An L., “Doppler optical micro-angiography for volumetric imaging of vascular perfusion in vivo,” Opt. Express, 17 8926 –8940 (2009). 10.1364/OE.17.008926 1094-4087 Google Scholar

24.

Tao Y. K., Kennedy K. M., Izatt J. A., “Velocity-resolved 3D retinal microvessel imaging using single-pass flow imaging spectral domain optical coherence tomography,” Opt. Express, 17 4177 –4188 (2009). 10.1364/OE.17.004177 1094-4087 Google Scholar

25.

Jiao S., Knighton R., Huang X., Gregori G., Puliafito C. A., “Simultaneous acquisition of sectional and fundus images with spectral-domain optical coherence tomography,” Opt. Express, 13 444 –452 (2005). 10.1364/OPEX.13.000444 1094-4087 Google Scholar

26.

Vergnole S., Lamouche G., Dufour M. L., “Artifact removal in Fourier-domain optical coherence tomography with a piezoelectric fiber stretcher,” Opt. Lett., 33 732 –734 (2008). 10.1364/OL.33.000732 0146-9592 Google Scholar

27.

Freedman D., Diaconis P., “On the histogramas a density estimator: L2 theory,” Z. Wahrscheinlichkeitstheor. Verwandte Geb., 57 453 –476 (1981). 10.1007/BF01025868 0044-3719 Google Scholar

Citation Download Citation

Lin An, Lin An, Hrebesh M. Subhush, Hrebesh M. Subhush, David J. Wilson, David J. Wilson, Ruikang K. Wang, Ruikang K. Wang, "High-resolution wide-field imaging of retinal and choroidal blood perfusion with optical microangiography," Journal of Biomedical Optics 15(2), 026011 (1 March 2010). https://doi.org/10.1117/1.3369811 . Submission:

JOURNAL ARTICLE
9 PAGES

DOWNLOAD PDF + SAVE TO MY LIBRARY