9 September 2014 Visualization of the ocular pulse in the anterior chamber of the mouse eye <italic<in vivo</italic< using phase-sensitive optical coherence tomography
Author Affiliations +
J. of Biomedical Optics, 19(9), 090502 (2014). doi:10.1117/1.JBO.19.9.090502
We report on a phase-based method for accurately measuring the ocular pulse in the anterior chamber in vivo. Using phase-sensitive optical coherence tomography with optimized scanning protocols and equations for compensating bulk motion and environmental vibrations, a high sensitivity of 0.9  μm/s minimal velocity is demonstrated at a wide detection band of 0 to 380 Hz. The pulsatile relative motion between cornea and crystalline lens in rodents is visualized and quantified. The relative motion is most likely caused by respiration (1.6 Hz) and heartbeat (6.6 Hz). The velocity amplitude of the relative motion is 10.3±2.4  μm/s . The displacement amplitudes at the respiratory and cardiac frequencies are 202.5±64.9 and 179.9±49.4  nm , respectively. The potential applications of the measurement technique can be found in the evaluation of intraocular pressure and the measurement of biomechanical properties of the ocular tissue, which are important in several ocular diseases.
Li, Ding, Ni, Xu, Zhao, Shen, Du, and Jiang: Visualization of the ocular pulse in the anterior chamber of the mouse eye in vivo using phase-sensitive optical coherence tomography

Due to the ocular hemodynamics in the elastic eyeball, particularly the vascular volume changes of the choroid, several ocular elements experience a pulsatile movement, i.e., ocular pulse (OP).1,2 OP is highly correlated with the biomechanical properties of the eye and may also act as a surrogate for changes in choroidal blood flow. The ability to noninvasively monitor the OP in vivo may play an important role in a number of ocular diseases such as glaucoma, age-related macular degeneration, and diabetic retinopathy. Fourier-domain optical coherence tomography (FDOCT) is a promising noninvasive method for high-resolution and dynamic eye imaging. Making use of the high-axial resolution of the FDOCT, the axial displacements of cornea and retina and the fundus pulse were measured by tracking the tissue position in sequential structural images.34.5 Generally, the precision of the intensity-based measurement is limited and is highly dependent on the accuracy of image segmentation. In order to suppress the segment-error-induced noise, the movement was processed using a bandpass filter at cardiac frequency.4,5 However, such a bandpass filter at certain frequencies is not desirable for the measurement of the minute OP in the anterior chamber (termed as aOP and defined as the relative motion between cornea and crystalline lens). As we know, the anterior chamber is filled with aqueous humor, and the aOP would be closely correlated to the intraocular pressure (IOP). Accordingly, similar to the IOP, the aOP may also be affected by several factors of different frequencies such as heartbeat and respiration.6 Therefore, it would be particularly attractive to capture the aOP at a high sensitivity and a wide frequency-band to retain the valuable information which would otherwise be lost.

In contrast, the phase shift ϕ during a time interval τ is inherently sensitive to the subwavelength displacement D of scattering particles in the phase-sensitive OCT (PhS-OCT),7,8, according to the relation expressed by ϕ=4πn/λ0×D, where n is the refractive index of the ocular tissue which is typically 1.36, and λ0 is the center wavelength used in OCT. The detectable displacement range of this phase-based method is bounded by the system signal-to-noise ratio (SNR) (lower limit, λ0/4πnSNR) and the π ambiguity (upper limit, λ0/4n).8 Depending on the choice of the time interval τ, different velocity ranges can be achieved. However, there are several factors that make it challenging to visualize and quantify the dynamic aOP by the phase-based method. First, due to involuntary head movements and microsaccades, the small relative cornea-lens motion is accompanied with a large and time-variant bulk motion. Accordingly, the time interval (determined by the B-frame rate in PhS-OCT) should be carefully selected to make the velocity range cover both the slow relative motion and the fast bulk motion. Second, since the phase shift is calculated between successive B-frames in PhS-OCT, the accumulated phase drift ϕdrift(x,i) caused by the random environmental vibrations is no longer negligible during such a long-time interval and requires reliable compensation. Here, we report for the first time the application of PhS-OCT to visualize the dynamic aOP in the rodent eye in vivo, along with optimized scanning protocols and equations for compensating bulk motion and random vibrations.

Mice (strain C57BL/6) of 10 to 12 weeks of age were used in this experiment. All mice were raised and treated in compliance with the Association for Research in Vision and Ophthalmology statement. Mice were anesthetized by intraperitoneal injection. A total of six eyes were scanned by PhS-OCT. The PhS-OCT system is based on a spectral-domain configuration, similar to the system in Ref. 7. A broadband superluminescent light source with a central wavelength of 850 nm and a spectral bandwidth of 100 nm was used to achieve an axial resolution of 3.2μm in air. A high-speed spectrometer is used to record the OCT spectrum, providing a 70-kHz line scan rate and a measured imaging range of 6.19mm in air at half Fourier space.

Based on the structural image, the cornea and lens components are generally 30dB above the noise floor [refer to Fig. 2(a)], i.e., the system phase noise can be considered as 0.03 rad. By tracking the cornea/lens position in the M-mode structural image [refer to Fig. 2(a)], the maximal velocity of the bulk motion can reach 107μm/s, whereas the velocity of the relative motion is in the range of several micrometers per second [refer to Fig. 2(c)]. To capture both the bulk and the relative motions, a B-frame rate of 780 fps was used to provide a theoretical detectable velocity range of 1.2 to 122μm/s at a wide detection band of 0 to 390 Hz (Nyquist theorem). Moreover, the slow scanning was stopped for repeated B-scanning in the horizontal meridians. In total, 3000 frames were collected covering a dynamic optical interval of 3.85s.

Figure 1(a) shows a wide-field OCT image of the full mouse eye using a full-range complex method, which has been detailed in Ref. 9 and is not elaborated here. In this study, the field of view is focused on the central part of anterior chamber [as marked by the rectangle in Fig. 1(a)] for dynamic imaging. As shown in Fig. 1(b), the raw phase cross-section ϕ0(x,z,i) can be generated based on the equation


where r(x,z,i) is the complex OCT cross-sectional signal at the i’th acquisition of successive frames. The acquisition number i can be related to the time t by t=iτ, and τ is the time interval between successive frames. * denotes the complex conjugation of a complex signal. arg represents the operation of phase extraction. According to the depth information in the phase images ϕ0(x,z,i), the raw cornea phase ϕc0(x,i) and the raw lens phase ϕl0(x,i) can be separated and, respectively, expressed as




where Δϕlc(x,i) is the phase difference corresponding to the relative motion between the cornea and the lens. ϕB(x,i) and ϕdrift(x,i) are the phase changes as a result of the bulk motion and the random vibrations, respectively. Therefore, as shown in Fig. 1(c), the relative phase difference Δϕlc(x,i) can be extracted based on the following equation:


where the raw cornea phase ϕc0(x,i) is reliably determined by a histogram-based phase selecting process.7,10

Fig. 1

(a) Wide-field structural cross-section in the horizontal meridian using a full-range complex method (refer to Ref. 9) to illustrate the architecture of the full mouse eye. Phase cross-section of the regional anterior chamber (as marked by the rectangle in a) before (b) and after (c) phase compensation.


In order to visualize the time-variant aOP, M-mode structural and phase images are generated along the center A-line of the B-frame (denoted as x0) as shown in Figs. 2(a) and 2(b). The anterior and posterior corneal boundaries as well as the anterior lens boundary are segmented based on the derivative of the intensity signal in the z-direction. As indicated in Fig. 2(b), the residual phase of the cornea component is around zero (green color), and the lens component presents a periodical variation versus time after compensation. Therefore, the compensated lens motion can be regarded as the relative cornea-lens motion, i.e., the aOP. To further quantify the aOP, the relative cornea-lens velocity (ΔVcl) is plotted versus time along the depth z0 [around the center part of the anterior lens, marked by the red bold arrows in Fig. 2(b)], based on the following equation:


as reported by the red curve in Fig. 2(c). The residual cornea velocity is also plotted with a root-mean-square level of 0.9μm/s (i.e., minimal detectable velocity), as shown by the blue dashed line in Fig. 2(c). According to the power spectrum analysis, the relative velocity ΔVcl is mainly composed of two frequency components at 1.6 and 6.6 Hz, respectively. The low-frequency component 1.6 Hz is most likely due to the mouse respiration, whereas the 6.6-Hz frequency component is induced by the heartbeat. The relative displacement (ΔDcl) pulse is readily obtained by integrating the velocity ΔVcl versus the time t, as shown in Fig. 3(a). The waveform of the pulsatile ΔDcl is quite similar to the IOP variations due to the heartbeat and respiration, as reported in Ref. 6. According to the frequency analysis, the relative displacement is further decomposed into two components: the heartbeat-induced and the respiration-induced, as respectively shown in Figs. 3(b) and 3(c). According to the statistics generated from a total of six mouse eyes, the velocity amplitude is 10.3±2.4μm/s. The displacement amplitudes at the respiratory and cardiac frequencies are, respectively, 202.5±64.9 and 179.9±49.4nm.

Fig. 2

Dynamic ocular pulse in the anterior chamber (aOP) image. (a) M-mode structural image generated along the center A-line of the B-frame. (b) M-mode structural image with phase encoded as hue. (c) The relative cornea-lens velocity (red line) and the residual cornea velocity (blue dashed line) versus time. The inset is the power spectrum of the relative velocity.


Fig. 3

Relative cornea-lens displacement. (a) Original displacement deduced from the velocity curve of Fig. 2(c) can be decomposed into two components: the heartbeat-induced (b) and the respiration-induced (c).


This preliminary study demonstrates the feasibility of using a phase-based method for aOP imaging/measurement in vivo. The method provides a high sensitivity of 0.9μm/s minimal velocity at a wide detection band of 0 to 380 Hz. Further research could be performed to correlate the aOP to the fundus pulse, the choroidal blood flow, the IOP, and the biomechanical properties of the ocular tissue.


We acknowledge financial supports from the Chinese Natural Science Foundation (61335003, 61275196, and 61327007), Zhejiang Province Science and Technology Grant (2012C33031), Zhejiang Provincial Natural Science Foundation of China (LY14F050007), and the Fundamental Research Funds for the Central Universities (2014QNA5017).



D. M. SilverR. A. Farrell, “Validity of pulsatile ocular blood flow measurements,” Surv. Ophthalmol. 38(Suppl), S72–S80 (1994).SUOPAD0039-6257http://dx.doi.org/10.1016/0039-6257(94)90049-3Google Scholar


J. C. Liet al., “The dynamic response of intraocular pressure and ocular pulse amplitude to acute hemodynamic changes in normal and glaucomatous eyes,” Invest. Ophthalmol. Visual Sci. 54(10), 6960–6967 (2013).IOVSDA0146-0404http://dx.doi.org/10.1167/iovs.13-12405Google Scholar


R. de Kinkelderet al., “Heartbeat-induced axial motion artifacts in optical coherence tomography measurements of the retina,” Invest. Ophthalmol. Visual Sci. 52(6), 3908–3913 (2011).IOVSDA0146-0404http://dx.doi.org/10.1167/iovs.10-6738Google Scholar


K. Singhet al., “Development of a novel instrument to measure the pulsatile movement of ocular tissues,” Exp. Eye Res. 91(1), 63–68 (2010).EXERA60014-4835http://dx.doi.org/10.1016/j.exer.2010.03.022Google Scholar


K. Singhet al., “Measurement of ocular fundus pulsation in healthy subjects using a novel Fourier-domain optical coherence tomography,” Invest. Ophthalmol. Visual Sci. 52(12), 8927–8932 (2011).IOVSDA0146-0404http://dx.doi.org/10.1167/iovs.11-7854Google Scholar


R. L. Cooperet al., “Continual monitoring of intraocular pressure: effect of central venous pressure, respiration, and eye movements on continual recordings of intraocular pressure in the rabbit, dog, and man,” Br. J. Ophthalmol. 63(12), 799–804 (1979).BJOPAL0007-1161http://dx.doi.org/10.1136/bjo.63.12.799Google Scholar


P. Liet al., “Pulsatile motion of the trabecular meshwork in healthy human subjects quantified by phase-sensitive optical coherence tomography,” Biomed. Opt. Express 4(10), 2051–2065 (2013).BOEICL2156-7085http://dx.doi.org/10.1364/BOE.4.002051Google Scholar


K. E. O’Haraet al., “Measuring pulse-induced natural relative motions within human ocular tissue in vivo using phase-sensitive optical coherence tomography,” J. Biomed. Opt. 18(12), 121506 (2013).JBOPFO1083-3668http://dx.doi.org/10.1117/1.JBO.18.12.121506Google Scholar


P. Liet al., “Extended imaging depth to 12 mm for 1050-nm spectral domain optical coherence tomography for imaging the whole anterior segment of the human eye at 120-kHz A-scan rate,” J. Biomed. Opt. 18(1), 016012 (2013).JBOPFO1083-3668http://dx.doi.org/10.1117/1.JBO.18.1.016012Google Scholar


V. X. D. Yanget al., “Improved phase-resolved optical Doppler tomography using the Kasai velocity estimator and histogram segmentation,” Opt. Commun. 208(4–6), 209–214 (2002).OPCOB80030-4018http://dx.doi.org/10.1016/S0030-4018(02)01501-8Google Scholar

Peng Li, Zhihua Ding, Yang Ni, Baishen Xu, Chen Zhao, Yi Shen, Chixin Du, Bo Jiang, "Visualization of the ocular pulse in the anterior chamber of the mouse eye <italic<in vivo</italic< using phase-sensitive optical coherence tomography," Journal of Biomedical Optics 19(9), 090502 (9 September 2014). http://dx.doi.org/10.1117/1.JBO.19.9.090502

Optical coherence tomography


In vivo imaging



Image segmentation


Back to Top