The present disclosure relates to exemplary methods and apparatus for obtaining information and providing analysis for biological tissues, and more particularly to exemplary methods and apparatus for achieving non-invasive in vivo biomechanical and biophysical characterization, including but not limited to elasticity and viscosity, of various biological tissues, including human and animal eye via, e.g., fast optical spectroscopy and microscopy based on Brillouin light scattering.
Ectasia is also one of the rare but serious adverse outcomes after LASIK (laser-assisted in situ keratomileusis) surgery. Currently about 1.5 million LASIK operations are performed annually in the U.S. As LASIK becomes increasingly popular, the incidence of post-LASIK ectasia has continued to increase. A promising therapeutic approach to corneal ectasia is increasing the stiffness of the stroma by crosslinking the naturally present collagen fibers in the cornea, a procedure known as corneal collagen crosslinking (CXL). The viscoelastic properties of the cornea are also known to affect the tonometry measurement of intraocular pressure.
As a consequence, the biomechanical properties may be an appropriate target for diagnosis and monitoring of onset and progression of cataract and presbyopia as well as corneal pathologies and treatments. For this reason, there has been a great deal of interest in measuring the mechanical properties of the lens corneal tissues for diagnosis and for monitoring of treatments. However, current techniques cannot detect such biomechanical changes in vivo in patients and animal models, seriously frustrating our efforts to develop understanding and treatment of the prevalent ocular problems.
Conventional techniques, from the traditional slit-lamp microscopy to newer imaging technologies (computer videokeratography, OCT, confocal microscopy, ultrasound, Scheimflug photography) are excellent in imaging the structure of cornea and lens but fail to provide their physiological and biomechanical information. Current clinical instruments, such as pachymetry (measuring thickness) and topography (mapping surface curvature), have been limited in screening patients at high risk of post-LASIK ectasia; patients with normal appearing corneas have developed the complication.
Several techniques have been used to characterize the mechanical properties of the cornea and lens ex vivo and in vivo. For example, comprehensive but destructive analysis has been performed by spinning cup, mechanical stretchers, stress-strain equipment or by inflation tests. Other mechanical testing methods include laser induced optical breakdown based on bubble creation and the ocular response analyze measuring corneal hysteresis on the surface without spatial information. Ultrasound is an attractive tool as it allows noninvasive methods such as elastography. Of particular note is ultrasound pulse-echo techniques and ultrasound spectroscopy, where pulsed or continuous-wave acoustic waves are launched onto the cornea, and the propagation speed and attenuation are measured to compute the viscoelastic moduli of the tissue. However, the ultrasound-based techniques have drawbacks of relatively low spatial resolution and measurement sensitivity.
Corneal ectasia refers to a bulging of the cornea, occurring when it is not strong enough mechanically to withstand the intraocular pressure. Ectasia may result from a degenerative disease called Keratoconus, affecting about 1 in 2000 people with a mean onset age of 15.4 years. (See, J. H. Krachmer et al., “Keratoconus And Related Noninflammatory Corneal Thinning Disorders,” Survey of Ophthalmology, vol. 28, pp. 293-322, 1984; and J. O. Jimenez et al., “Keratoconus: Age of onset and natural history,” Optometry and Vision Science, vol. 74, pp. 147-151, March 1997). Ectasia is also one of the rare serious adverse outcomes after LASIK surgery that results in corneal thinning. As LASIK becomes increasingly popular (about 1.5 million operations in the U.S. each year), the incidence of post-LASIK ectasia has continued to increase. Current clinical instruments, such as Pachymetry (measuring the thickness) and topography (mapping surface curvature), have been limited in screening patients at high risk of post-LASIK ectasia.
Concern about development of post-LASIK ectasia and imperfect screening tests prevents ˜15% of patients seeking laser vision correction from benefiting from this procedure in the US. On the other hand, post-LASIK ectasia is reported in ˜0.7% of the remaining cases in patients with normal appearing pachymetry and topography and no other risk factors. (See I. G. Pallikaris et al., “Corneal Ectasia Induced By Laser In Situ Keratomileusis,” J Cataract Refract Surg, vol. 27, pp. 1796-1802, 2001; P. S. Binder et al., “Keratoconus and corneal ectasia after LASIK,” Journal of Refractive Surgery, vol. 21, pp. 749-752, November-December 2005; and Y. S. Rabinowitz, “Ectasia after laser in situ keratomileusis,” Current Opinion in Ophthalmology, vol. 17, pp. 421-426, October 2006). This situation has led to multimillion dollar lawsuits (see, G. McDermott, “Anatomy of a lawsuit.,” Cataract and Refractive Surgery Today, vol. 2005, pp. 93-118, 2005), and has contributed to the decision by the FDA to reevaluate LASIK surgery safety. (See FDA, “http://www.fda.gov/NewsEvents/Newsroom/PressAnnouncements/2009/ucm186858.htm,” 2009). There is a growing consensus that the mechanical properties of the cornea, e.g. elastic modulus or stiffness, given their connection to the pathophysiology of the condition, are highly suited to early identify patients at risk for developing ectasia. (See Y. Rabinowitz, “Ectasia after laser in situ keratomileusis,” Current opinion in ophthalmology, vol. 17, pp. 421-427, 2006; and J. W. Ruberti et al., “Corneal Biomechanics and Biomaterials,” in Annual Review of Biomedical Engineering, Vol 13. vol. 13, M. L. D. J. S. G. M. L. Yarmush, Ed., 2011, pp. 269-295).
A further therapeutic approach to corneal ectasia can increase the stiffness of the stroma by crosslinking the naturally present collagen fibers in the cornea—a procedure known as corneal collagen crosslinking (CXL). For this procedure, monitoring the mechanical properties of the ocular tissue before, during and after, can assist with optimizing the outcome and success of the intervention.
Measuring the mechanical properties of the cornea can also be beneficial for glaucoma screening. Glaucoma is the second highest cause of blindness in the world. About 2.2 million Americans have glaucoma and half of them are not aware of their condition. Left untreated, glaucoma leads to blindness; if diagnosed early and managed, 90% of the patients avoid blindness. Elevated intraocular pressure (IOP) is a known risk factor for glaucoma and is currently used as major screening parameter. To evaluate IOP, applanation tonometry measures the external pressure required to applanate the cornea. However, applanation tonometry can assume the cornea to be thin, spherical and elastic; therefore, pressure readings can be biased by actual corneal thickness, curvature and stiffness. While pachymetry and tomography can measure conical thickness and curvature, stiffness is not currently accounted for. Modeling studies show that lack of stiffness information introduces an uncertainty in TOP estimation of 2 to 3 mmHg (enough to misclassify the risk category of patients). The effect can be more pronounced in patients who underwent refractive surgery or crosslinking. Without corneal stiffness information, IOP measurements are therefore unreliable.
As a consequence, the biomechanical properties may be an appropriate target for diagnosis and monitoring of onset and progression of cataract and presbyopia as well as corneal pathologies and treatments. For this reason, there has been a great deal of interest in measuring the mechanical properties of the lens corneal tissues for diagnosis and for monitoring of treatments. However, prior techniques may not be able to detect such biomechanical changes in vivo, e.g., in patients and animal models, seriously frustrating efforts to develop understanding and treatment of the prevalent ocular problems.
Several techniques have been developed to characterize the mechanical properties of the cornea and lens ex vivo and in vivo. For example, comprehensive but destructive analysis has been performed by spinning cup, mechanical stretchers, stress-strain equipment or by inflation tests. Other prior mechanical testing methods can include laser induced optical breakdown based on bubble creation and the ocular response analyze measuring corneal hysteresis on the surface without spatial information. Ultrasound can be a possible approach, since it facilitates noninvasive methods such as elastography. For example, ultrasound pulse-echo techniques and ultrasound spectroscopy can be used, where pulsed or continuous-wave acoustic waves are launched onto the cornea, and the propagation speed and attenuation are measured to compute the viscoelastic moduli of the tissue. However, the ultrasound-based techniques have certain drawbacks, such as, e.g., a relatively low spatial resolution and a particular measurement sensitivity.
Brillouin scattering results from the interaction of an incident optical wave with hypersonic acoustic phonons inside the medium under test. In Spontaneous Brillouin scattering the acoustic phonons are inherently present in the material due to thermally-induced density fluctuations. However, the Brillouin process can be enhanced or forced by using multiple optical pump waves with frequencies separated by those of the acoustic phonons in the medium. The technique that analyzes the feature of Brillouin scattering signal is known as Brillouin spectroscopy. Various techniques to detect the Brillouin signal have been known in the art and are widely applied in physics, material science, and mechanical engineering.
The magnitude and frequency (spectrum) of the Brillouin scattered light are determined by hypersonic acoustic phonons that are inside the material, the latter being closely related to the mechanical properties of the medium, such as elastic and viscous modulus. These mechanical properties therefore can be measured by Brillouin spectroscopy. In addition, Brillouin spectrum can also be dependent on temperature, pressure, density and refractive index; hence, under specific conditions the analysis of Brillouin spectrum can be used to measure these quantities.
Brillouin spectroscopy of the cornea and lens has been performed using tissues ex vivo. (See Y. S. Rabinowitz, “Ectasia after laser in situ keratomileusis,” Current Opinion in Ophthalmology, vol. 17, pp. 421-426, October 2006). However, the significant potential of using Brillouin scattering for ocular biomechanics has not been fully exploited and in vivo measurement has not been demonstrated. From a technological point of view, the long acquisition times required by the spectral analysis have limited the technique to point spectroscopy while the limited extinction of spectral instruments has made it difficult to discriminate Brillouin scattering signal from the elastic scattering or the scattering from various optical components.
Although Brillouin spectroscopy has been used for material characterization, structural health monitoring, and environmental sensing, the combined requirements of spectral resolution (10−3 nm) and spectral extinction (to suppress strong elastic scattering) forced using slow sequential spectrometers, scanning Fabry-Perot (FP), which acquire a single Brillouin spectrum likely in about a number of minutes to hours. (See B. Culshaw et al., “Smart structures and applications in civil engineering,” Proc. of the IEEE, vol. 84, pp. 78-86, 1996; G. D. Hickman et al., “Aircraft laser sensing of sound velocity in water: Brillouin scattering,” Rem. Sens. Env., vol. 36, p. 165, 1991; and J. R. Sandercock, “Some recent developments in Brillouin scattering,” Rca Review, vol. 36, pp. 89-107, 1975). This can limit the technology to point-sample analysis.
In the past, a spectrometer has been provided which can use a modified FP etalon called a virtually-imaged-phased-array (VIPA), that facilitated parallel detection. (See G. Scarcelli and S. H. Yun, “Brillouin Confocal Microscopy for three-dimensional mechanical imaging,” Nature Photonics, vol. 2, pp. 39-43, 2008; and M. Shirasaki, “Large angular dispersion by a virtually imaged phased array and its application to a wavelength demultiplexer,” Opt. Lett., vol. 21, pp. 366-368, March 1996). This greatly reduced the acquisition time to about 1-10 seconds per spectrum in transparent materials. (See G. Scarcelli and S. H. Yun, “Brillouin Confocal Microscopy for three-dimensional mechanical imaging,” Nature Photonics, vol. 2, pp. 39-43, 2008) This positive result, however, was still not sufficient to achieve an appropriate Brillouin imaging in tissue in vivo, or rapid monitoring of tissue mechanical properties.
Accordingly, there may be a need to overcome at least some of the issues and/or deficiencies described herein above.
To address and/or overcome the above-described problems and/or deficiencies, exemplary embodiments of systems, arrangements and processes can be provided that are capable of obtaining information and providing analysis for biological tissues. For example, exemplary methods and apparatus can achieve non-invasive in vivo biomechanical and biophysical characterization, including but not limited to elasticity and viscosity, of various biological tissues, including human and animal eye via, e.g., fast optical spectroscopy and microscopy based on Brillouin light scattering. In certain exemplary embodiments of the present disclosure, apparatus and methods can be provided to obtain non-invasively the information about intrinsic biomechanical properties of biological tissue (e.g., including but not limited to human and other animal eyes), which can be utilized in a wide range of applications in clinical diagnosis of, e.g., ocular problems, research, development of treatments using animal models, etc.
According to certain exemplary embodiments of the present disclosure, it is possible to perform a spectral analysis to facilitate an in vivo tissue characterization including, e.g., a) sub-second Brillouin spectrum acquisition, b) higher than 60 dB rejection of elastic scattering, and c) sub-GHz elasticity sensitivity. For example, by combining an exemplary spectrometer having certain exemplary characteristics to achieve such exemplary results with a confocal microscope, it can be possible to characterize ocular tissue, lens and cornea in vivo, e.g., in mice and human patients, as well as in and for skin tissue in animals and human. In addition, using such exemplary embodiments, “elasticity imaging” can be accomplished with a high spatial resolution. According to the exemplary embodiments of the present disclosure, there is a benefit over ultrasound techniques, as significantly higher resolution and sensitivity can be obtained. Indeed, it is possible to probe the elasticity of the tissues non-invasively with a microscopic resolution.
Through elasticity imaging, it is possible to utilized certain exemplary protocols, e.g., to extract quantitative biomechanical parameters that can characterize the biophysical/biomechanical state of the biological tissue under consideration. In this manner, health status, risk to develop conditions and outcome of the exemplary procedures can be assessed beyond structural information. Therefore, according to the exemplary embodiments of the technique which uses Brillouin principles, it is possible to enhance management of a variety of ocular problems by, e.g., providing currently inaccessible, biophysical and bio-structural data of the lens and cornea.
According to an exemplary embodiment of the present disclosure which utilizes Brillouin spectroscopy and imaging, it is possible to overcome many of the current state of the art technology limitations and/or deficiencies, such as slit lamp, pachymetry, and topography. Other imaging devices, such as the Pentacam (rotating Schiemflug imaging), generally provide only static information about the lens and cornea and likely do not provide biophysical dynamic information about the cornea or lens. Likewise, the Ocular Analyzer (corneal hysteresis) generally provides only limited information about the cornea regarding its elasticity, and no imaging information or data about the lens. In certain exemplary embodiments of the present disclosure, it is possible to provide quantitative values, and a high quality image of biophysical properties of lens and cornea. This can facilitate in screening of high at-risk patients of corneal ectasia and diagnosing diseases such as keratoconus and cataracts.
According to further exemplary embodiments of the present disclosure, it is possible to utilize to analyze ocular tissue elasticity behavior versus frequency of induced stress and validating Brillouin microscopy against high-end mechanical tests. In addition, an exemplary in vivo biomechanical analysis can be tested against the known natural stiffening that occurs in aging small animals to demonstrate the ability of this technology of providing information of clinical relevance.
In one exemplary embodiment of the present disclosure, it is possible to utilize a spectroscopic analysis that is performed on a high-speed basis, with high spectral resolution and high extinction. The use of such exemplary technique can be done with an ophthalmic slit lamp to facilitate a mechanical image formation. For example, a pump beam can be scanned over a sample (e.g., a biological tissue, such as an eye) through an objective lens, and Brillouin-shifted optical waves are detected to characterize the Brillouin spectrum. The measured spectral features of the Brillouin signal can represent the contrast mechanism for imaging; an image can be obtained by use of an exemplary lookup table and/or an appropriate processing computational routine. Microscopically-resolved three-dimensional images of biomechanical properties of cornea, crystalline lens, aqueous and vitreous humor of the human eye, as well as other tissues can be obtained by using, e.g., a high numerical-aperture objective lens and confocal detection.
According to another exemplary embodiment of the present disclosure, a biomechanical analysis of the human eye can be combined with Raman spectroscopy, fluorescence and auto-fluorescence microscopy, as well as reflectance microscopy. The exemplary co-registration of multiple contrast mechanisms can facilitate, e.g., a complete or near-complete characterization of the biophysical properties, such as structure, elasticity, chemical composition, functional abilities of the human eye in three dimensions with micron-scale resolution, etc.
According to still further exemplary embodiments of the present disclosure, enhanced Brillouin scattering can be facilitated using multiple pump beams, auxiliary ultrasound, or contrast-enhancing nanoparticles. Various pump and probe techniques, heterodyne or spectrometer-based detection techniques can be used. Since the Brillouin shift in human tissues can range typically from about 10 MHz to 10 GHz, the direct electrical detection of the acoustic wave may also be possible.
According to particular exemplary embodiments of the present disclosure, arrangements and methods can be provided for determining information associated with at least one section of at least one biological tissue. For example, it is possible (e.g., using at least one first arrangement) to provide at least one first electro-magnetic radiation to the section(s) of the biological tissue(s) in vivo so as to interact with at least one acoustic wave in the biological tissue(s). At least one second electro-magnetic radiation can be produced based on the interaction. With at least one second arrangement, at least one portion of the second electro-magnetic radiation(s) can be received, and the information associated with the section(s) of the biological tissue(s) can be determined based on the portion of the second electromagnetic radiation(s).
For example, the section(s) of the biological tissue(s) can include an ocular tissue of an eye, which can be cornea, aqueous humor, crystalline lens, vitreous humor, or retina in the eye in vivo. The first arrangement(s) can include a radiation emitting source arrangement which can be configured to provide the first electro-magnetic radiation(s), the first radiation(s) can have a wavelength in the range of about 450-1350 nm with a spectral width less than 1 GHz. The second arrangement(s) can include a spectrally-resolving arrangement which can be configured to measure at least one spectral characteristic of the portion(s) of the second electromagnetic radiation. The second arrangement(s) can obtain the information associated with the section(s) of the biological tissue based on the spectral characteristic(s).
The spectrally-resolving arrangement can include a spectrometer that is configured to disperse the spectrum of the second electromagnetic radiation(s). The spectrometer can have a spectral extinction efficiency greater than about 60 dB. The spectrometer can include at least one virtually imaged phased array (VIPA) etalon that can be configured to disperse the spectrum of the second electromagnetic radiation(s). The spectrometer can include an apodization filter and/or an apodized etalon that is configured to provide a spectral extinction efficiency greater than about 60 dB. 10. The spectral characteristic(s) can be at least one frequency difference between the first electro-magnetic radiation(s) and the portion(s) of the second electro-magnetic radiation(s) The frequency difference(s) can be associated with a propagation speed of the at least one acoustic wave, and the frequency difference can be in the range of about 2 GHz to 20 GHz. The information associated with the section of the biological tissue can be at least one image or a spatially-resolved map of the at section(s) based on at least one parameter associated with the frequency difference. The parameter(s) can include a visco-elastic modulus. The spectral characteristic can be at least one spectral line width of the portion(s) of the second electro-magnetic radiation. The spectral line width can be associated with a propagation attenuation of the acoustic wave(s), and the frequency line width can be in the range of about 0.3 GHz to 3 GHz.
In another exemplary embodiment of the present disclosure, the information associated with the section(s) can be maximum, average, and/or rate of variation of at least one parameter related to the frequency difference over the at least section. At least one third arrangement can be provided that is configured to facilitate positioning the biological tissue(s) with respect to the first electro-magnetic radiation(s). The third arrangements) can include an imaging arrangement configured to measure at least one position of the first electro-magnetic radiation(s) with respect to the biological tissue(s). The second arrangements) can determine the information based on the position(s). At least one fourth arrangement can be provided which is configured to receive the first electromagnetic radiation(s), and generate a fourth radiation based on an interaction of the first electromagnetic radiation(s) with an acoustic wave in a reference material. A spectrum of the fourth radiation can be used to determine the information by the second arrangement(s).
In still another exemplary embodiment of the present disclosure, the portion(s) of the second electromagnetic radiation(s) can be provided by Brillouin scattering using the first electro-magnetic radiation in the biological tissue. The first arrangement(s) can be further configured to provide the first electromagnetic radiation to at least one segment which is at least one point, at least one line and/or an area on or within the biological tissue. The second arrangement(s) can be configured to determine the information from the segment(s) in less than 0.4 seconds (or in less than 1 second and/or at least 0.1 seconds); where the total optical power of the first electro-magnetic radiation is less than 1 mW. The first arrangement(s) can be further configured to provide the first electromagnetic radiation(s) to a further segment, and the second arrangement(s) can be further configured to generate at least one image for the segment and the further segment. The first arrangement(s) can cause a movement of the first electromagnetic radiation(s) from the segment to the further segment in a transverse direction and/or a longitudinal direction with respect to the section(s).
For example, the information can include (i) a biomechanical property, (ii) a stiffness, and/or (iii) a cross-linking of the biological tissue, (iv) a stiffness, (v) an accommodation power, (vi) a presbyopia, or (vii) a cataract of the crystalline lens, (viii) stiffness, (ix) a keratoconus, or (x) a risk of ectasia for a refractive surgery, (xi) collagen crosslinking of the cornea, and/or (xii) intraocular pressure of the eye. At least fourth arrangement can be provided that is configured to produce at least one image associated with the characteristic of the biological tissue. The point can include a plurality of points, and the second arrangement(s) can determine the information from the plurality of points. The second arrangement(s) can determine the information from the plurality of points in less than 0.4 seconds.
In still another exemplary embodiment, apparatus and method can be provided according to the present disclosure. For example, using at least one first arrangement, it is possible to provide at least one first electro-magnetic radiation to at least one portion of at least one sample so as to generate at least one second electro-magnetic radiation with multiple spectral peaks. Further, using at least one second arrangement, it is possible to receive at least one portion of the second electro-magnetic radiation(s) so as to simultaneously acquire spectrum of the portion of the second electromagnetic radiation(s)n having the multiple spectral peaks in a range of about 2 GHz to 200 THz.
The portion(s) of the second electro-magnetic radiation(s) having at least one of the spectral peaks can be associated with Brillouin scattering, Raman scattering, fluorescence and/or luminescence, and the portion(s) of the second electro-magnetic radiation(s) having another one of the spectral peaks can be associated with at least different one of Brillouin scattering, Raman scattering, fluorescence and/or luminescence. The sample can be a biological tissue. The first electro-magnetic radiation(s) can be directed to the biological tissue in vivo.
These and other objects, features and advantages of the present invention will become apparent upon reading the following detailed description of embodiments of the disclosure, when taken in conjunction with the appended claims.
Further objects, features and advantages of the present disclosure will become apparent from the following detailed description taken in conjunction with the accompanying drawings showing illustrative embodiments of the present disclosure, in which:
is a diagram of an exemplary optical coherence tomography (OCT) system which can perform scan of a healthy eye
Throughout the figures, the same reference numerals and characters, unless otherwise stated, are used to denote like features, elements, components or portions of the illustrated embodiments. Moreover, while the subject invention will now be described in detail with reference to the figures, it is done so in connection with the illustrative embodiments. It is intended that changes and modifications can be made to the described embodiments without departing from the true scope and spirit of the subject disclosure as defined by the appended claims.
According to an exemplary embodiment of the system, method and apparatus according to the present disclosure, it is possible to perform a Brillouin microscopy in ocular tissue in vivo, which can be valuable in ocular biomechanical characterization in diagnosing and treating ocular problems, as well as developing novel drugs or treatments.
There are four anatomical sites in the eye: For example, the cornea is a thin (e.g., less than 1 mm) tissue composed by different layers of varying mechanical strength. The aqueous humor is a liquid with similar properties to water that fills the anterior chamber of the eye. The crystalline lens is a double-convex sphere composed by many layers of different index of refraction, density and stiffness. The vitreous humor is the viscous transparent liquid that fills the posterior chamber of the eye.
Brillouin light scattering in a tissue or any other medium usually arises due to the interaction between an incident light and acoustic waves within the matter. For example, a probe light having a frequency v and a wavelength λ can be used, which may be provided to a sample. In a Spontaneous Brillouin process, the acoustic waves or acoustic phonons are naturally present due to thermal fluctuations. Such fluctuations propagate through the medium in the form of acoustic waves. These acoustic waves can generate periodic modulations of the refractive index. Brillouin scattering can be generated by at least one or many acoustic waves or acoustic phonons, which form phase-matched index modulation.
The electromagnetic radiation 110 can be directed to the eye 120 to probe various portions of ocular tissues, including but not limited to, cornea 122 and a crystalline lens 124. For example, an imaging lens 130 can be used to focus the electromagnetic radiation 110 onto a small spot. The imaging lens 130 can be a spherical convex lens, aspheric lens, objective lens, theta lens, or cylindrical lens for line focusing.
To scan the axial position of the focus within the ocular tissues, the imagining lens 130 can be mounted on a translation stage 134. Alternatively or in addition, a tunable element that can change a divergence of the probe radiation can be employed. To scan the transverse position of the focus, a one- or two-axis beam scanner 140 can be employed. The exemplary scanner 140 can include a galvanometer-mounted mirror, MEMS mirror, translation stages, spatial light modulator, etc.
An acousto-optic interaction in the tissue can give rise to light/radiation scattering, thereby generating at least one second electromagnetic radiation. Several mechanisms for light/radiation scattering are known in the art, including Rayleigh and Mie scattering, Raman scattering, and Brillouin scattering. While biological tissues support these scattering mechanisms, Brillouin scattering is directly associated with the acoustic waves in the medium. A portion of such one or more second electromagnetic radiations can be collected by the imaging lens 130.
The exemplary system of
The beam scanning apparatus can be mounted on an ophthalmic slit lamp instrument, and a very low power auxiliary laser beam can be used so that the operator can view the eye and aim the probe laser at the target location in the eye. For small animal studies, another exemplary embodiment of the present disclosure can be used without the slit lamp. If not mounted on the slit lamp, the exemplary system can further comprise a third arrangement, which can be configured to facilitate positioning the eye 120 with respect to the first electro-magnetic radiation(s), and/or the probe light. Preferably, the third arrangement can include at least one of the following features: a forehead rest, a chin rest, an eye fixation beam, and/or a bite bar to maintain repeatable position over long-term analysis, etc. For example, the human interface 180 can employ a camera to measure at least one position of the at least one first electro-magnetic radiation with respect to the at least one ocular tissue. This type of beam guiding arrangement can facilitate aiming the probe beam and provide the position information of the focus, which can be used in making of Brillouin images or the spatial map of the biomechanical properties of the ocular tissue.
Several exemplary configurations for illumination collection of light from ocular tissue can be employed. In an epi-detection configuration, the interacting probe and Brillouin-scattered lights travel in the nearly opposite directions. Alternatively or in addition, a dual-axis configuration can be employed, where the probe and scattered light can form a finite angle. A theta lens or a microscope objective lens with a long working distance can be used for the exemplary Brillouin detection. For coarse resolutions, a collimated pump beam with a relatively small beam diameter can be used. For three dimensional resolutions, the pump/probe beams can be focused into the eye by the use of objective lenses.
In this respect, cornea, lens, aqueous and vitreous humors do not absorb light in the visible and near infrared wavelengths, while the retina is sensitive to such radiation. Moreover, the eye itself can be considered as an optical system that may tend to focus parallel beams onto the retina. As a consequence, to increase and/or maximize the input power, and thus Brillouin signal, while avoiding eye damage, the focusing strategy should be carefully considered. When possible, it may be preferable to focus the light onto the cornea and on to the lens with high numerical aperture objectives because their fast defocusing will minimize the light intensity delivered to the retina.
The selection of the focusing/collecting strategy and the objective lens generally takes into account other factors. For example, in the epi-detection configuration, the backward-propagating Brillouin light can be detected by the same objective lens used for illumination. In dual axis configuration, e.g., illumination and Brillouin light can be either detected by the same lens (small angles) or detected by two different lenses (larger angles, longer working distance). Epi-detection can achieve the maximum collectable signal at given illumination power. In particular, the advantage in collection efficiency can be calculated to be, e.g., 2*sin(2θ)/3*NA, where 2θ is the dual-axis angle and NA is the numerical aperture of the objective lens. Moreover, epi-detection configuration can facilitate the beam scanning over the sample. On the other hand, the dual axis configuration can increase and/or maximize axial resolution at given NA. In addition, a dual axis configuration greatly reduces back-reflections from optical components, thus relaxing the requirements on extinction from the spectrometer.
In terms of objective lenses, exemplary lenses with low numerical aperture (NA) result in a low transverse resolution, but the longitudinal interaction length is long and well defined. Objective lenses with high NA, generally give better transverse and axial resolution. Since the interaction can be made over a large solid angle, the line width of the scattered light broadens affecting the strength of Brillouin signal and the accuracy of the spectral analysis.
Exemplary confocal techniques can be used to enhance depth sectioning. For the same or similar purpose, a single mode fiber can be used as confocal pinhole. Acting as tight spatial mode filter, the fiber can facilitate a strict confocal imaging, thereby minimizing any stray or spurious unwanted radiation. On the other hand, to increase the signal level, a multi-mode collection can be employed. For example, a double-clad fiber with a large inner cladding can collect light with the collection efficiency about 50-200 times higher than that of a single mode fiber currently used in our prototype. Double-clad fibers generally enhance the signals in fluorescence endoscopy and white-light scattering. This exemplary approach can reduce confocal sectioning, therefore the optimal balance between improved SNR and confocal sectioning can be found.
According to one exemplary embodiment of the present disclosure, the second arrangement 150 can employ at least one spectral analysis unit, such as a spectrometer, a monochromator, fixed or scanning spectral filters, or other devices known in the art. The second arrangement 150 can be configured to measure various properties of the second electromagnetic radiation 144, including but not limited to, the center frequency and width of its spectrum, as well as the intensity and polarization of the electrical field. In particular, the frequency difference between the first electromagnetic radiation(s) 110 entering the tissues and the portion of the second electromagnetic radiation(s) 144, which includes the Brillouin scattered light, can be of importance.
The exemplary system/apparatus can further comprises a fourth arrangement 180 configured to display the information associated with the at least one portion of the ocular tissue in the eye in vivo. The displayed information may include but is not limited to the Brillouin frequency shifts, Brillouin line width, Brillouin images, and the hypersonic viscoelastic moduli, as well as parameters, such as the mean values or slopes, calculated from the Brillouin images or the spatial maps of the viscoelastic properties.
The exemplary system/apparatus can further utilize a fifth arrangement 190 to provide at least one frequency reference. For example, the fifth arrangement 190 can be configured to receive at least one portion of the first electromagnetic radiation through the beam splitter 142, and reemit Brillouin scattered light with at least one, and possibly multiple, spectral peak(s). For example, the frequency reference 190 can include at least one reference material, solids or liquids, with known Brillouin frequency shifts. Alternatively or in addition, the frequency reference 190 can be a light source emitting an electromagnetic radiation at a wavelength locked to the wavelength of the probe light source 100. In both case, the electromagnetic radiation from the frequency reference 190 is directed to the second arrangement 150. An optical switch 192 can be employed to gate the intensity of the electromagnetic radiation. The reference frequency can assist with calibrating the spectral analysis unit in the second arrangement 150, facilitating the spectral analysis.
In another exemplary embodiment, additional arms can be added to the microscope to measure the Brillouin scattering signal from known materials. Having two additional reference materials can be sufficient to have a constantly calibrated instrument that can automatically correct for small variation of inside the spectrometer due to temperature instabilities or mechanical drifts.
The frequency shift νB of the Brillouin scattered light with respect to the probe light 110 can be given by
Where n is the local refractive index in the interrogated tissue, V is the speed of the acoustic wave in the sample, and θ is the scattering angle, i.e. the angle between the incident and the scattered light, such as in the dual-axis geometry. In an epi-backward detection configuration, θ=π can be a reasonably good approximation. In typical soft tissues, the speed of the acoustic wave can range from about 1000 to 3000 m/s, and the Brillouin frequency shifts can typically be between 2 and 20 GHz, depending on the wavelength.
The intrinsic spectral width or line width of the Brillouin scattered light can be given by:
where α is the attenuation coefficient of the acoustic wave in the sample.
The longitudinal complex elastic modulus, M=M′+iM″, where the real part M′ refers to the elastic modulus and the imaginary part M″ is the viscous modulus is given by:
M′=ρV
2; (5)
M″=2ρV3α/νB. (6)
Although the physical mechanism of Brillouin scattering can be known, the meaning of Brillouin signatures has not been studied in the context of its biomechanical significance of complex biological material. The biomechanical information extracted from Brillouin spectroscopy may not be the familiar stiffness that is felt upon palpation and is not the standard elastic modulus that can be obtained from a macroscopic stress-strain test. The Brillouin viscoelastic moduli can be defined in Equations (5) and (6) represent the tissue properties at the hypersonic GHz frequencies. Most soft tissues, including the corneal tissues and crystalline lens, exhibit viscoelastic properties characterized by frequency-dependent moduli. Slower relaxation processes have little time to respond to fast mechanical or acoustic modulation, such as GHz acoustic phonons, and thus hardly contribute to the “softness” of the material. As a consequence, modulus tends to increase with frequency. In addition, the propagation of acoustic phonons can be governed by the longitudinal modulus, which is typically much higher than the Young's or shear modulus owing to the incompressibility (i.e. Poisson's ratio ˜0.5) of water. The effects, finite relaxation time and low compressibility, can provide qualitative explanation for the observed large difference in modulus between the Brillouin and standard mechanical tests.
A correlation experiment has been performed with bovine and porcine lenses. Fresh porcine and bovine lenses have been cut at various ages (from 1 to 18 months) into small pieces of the size our mechanical equipment could handle. The mean Brillouin modulus was calculated from the 3D measurement of Brillouin spectrum and the estimated density and refractive index.
A comparison to Young's modulus measured by a conventional stress-strain test indicates a correlation 230 between Brillouin (M′) and quasi-static modulus (G′) for both porcine and bovine tissues. High correlation (R>0.9) was obtained in curve fit to a log-log linear relationship: log(M′)=a log(G′)+b, where the fitting parameters were a=0.093 and b=9.29 for porcine tissues and a=0.034 and b=9.50 for bovine tissues.
Therefore, the exemplary measurement of the spectral characteristics of the Brillouin scattered light provides the information about the biomechanical properties of the ocular tissue. The useful information obtained by the Brillouin measurement using the exemplary embodiments of the present disclosure includes but is not limited to the acoustic speed, acoustic attenuation coefficient, Brillouin elastic modulus, Brillouin viscous modulus, and electrostriction coefficient. As described further below, by scanning the focus within the tissue different spatial locations can be probed, which can provide the information in a spatially resolved manner. This spatial information can in turn be useful to evaluate for the diagnosis of the mechanical integrity or health of the ocular tissue.
The index of refraction and acoustic speed of a given material are generally dependent on the local temperature and pressure. This dependence may be obtained for the analysis of inflammatory or pathologic states in the eye via the measurement of the temperature or ph-value in the aqueous and vitreous humors. The magnitude of the Brillouin scattered radiation is related to the coupling of acoustic and optical energy inside the sample, which is related to the material properties, such as the electrostriction coefficient.
For example, performing Brillouin microscopy in the eye may not be trivial. Brillouin spectroscopy in the human eye can generally use a low power of illumination and a careful design of the focusing microscope. The potential risks associated with light exposure for Brillouin analysis should be considered. Maximum permissible exposure (MPE) is defined as the highest power or energy density that can be admitted to the eye without causing a biohazard. MPE corresponds to 10% of the dose that has a 50% chance of creating damage in worst-case scenario conditions. In the Brillouin scanner, the laser light can be focused in the cornea and diffused onto the retina. According to the International Commission on Non-Ionizing Radiation Protection (ICNIRP), for cw sources in the wavelength region 400-1050 nm the exposure limit for cornea-lens thermal safety is 4 W/cm2; i.e. MPE=32 mW in a 1-mm-diameter zone (0.79 mm2 in area)[19]. The 1-mm aperture for irradiance averaging can be based on the thermal modeling[20], which shows that the temperature rise is independent of the beam size up to 1 mm due to rapid thermal conduction (−0.5 W/miC) in the cornea and lens.
As for the retina, the dominant mechanism of retinal damage can be thermal for an exposure time (T) longer than 0.25 s. Considering the distance 17 mm between the lens and retina and α=0.2, i.e. NA=0.1, the beam size on the retina is >3 mm. Applying the same thermal limit of 4 W/cm2, we calculate MPE=˜300 mW. According to the guidelines from American National Standard Institute, the exposure limit for the retinal thermal hazard has been expressed as 1.8×10−3 CACET0.25 [W/cm2], where CA=1-1.5 (for λ=400-800 nm), CE=267 (for α=0.2; i.e. NA=0.1), and the aperture size of 7 mm (area of 0.38 cm2)[21]. The MPE is calculated to be 183 mW for T=1 s and 66 mW for T=60 s, which is consistent with the above analysis. Taken together, the Brillouin scan in the cornea and lens does not pose risk to the human eye, if laser light power is below 3 mW.
The low illumination power places a requirement on the sensitivity of the instrument capturing Brillouin scattered light. Moreover, the back-reflected components from optical components inside the microscope as well as back-scattered signals from the eye due to the difference in index of refraction between the various regions of the eye, in most situations and for most instruments would overshadow Brillouin scattering signal. This can place a requirement on the extinction of the instrument capturing Brillouin scattered light. Additionally, Brillouin signatures from ocular tissue will range between 9 and 15 GHz for the lens and 7.5 to 10 GHz for the cornea; thus, the spectral resolution of the instrument can be extremely high (e.g., less than about 1 GHz) in order to detect relevant changes in the biomechanics of the eye.
According to certain exemplary embodiment of the present disclosure, these simultaneous can be utilized for an exemplary system/arrangement to detect weak Brillouin scattering from the eye, spectrally analyze it and reconstruct a two-dimensional or three-dimensional, microscopic, mapping of the mechanical properties of the eye.
For example, a spectrometer with high resolution, high sensitivity and high extinction in the spectral analysis can be used with the exemplary system/apparatus. Previously, Fabry-Perot (FP) interferometry has been used for the spectral analysis of Brillouin signal in both scanning and angle-dispersive configuration. Both of these methods are generally slow, because of their intrinsically limited throughput. In particular, the resolution performance of spectrally dispersive elements are characterized by a parameter known as finesse which can be thought as how many spectral components can be resolved. However, the maximum amount of light that is forwarded to the detector in a Fabry-Perot interferometer is likely inversely proportional to its finesse; if n is the finesse of the Fabry-Perot interferometer, maximum 1/n of the input light can be sent to the detecting device. This tradeoff between achievable spectral resolution and throughput of the instrument is a fundamental limit of prior art which does not include other practical loss mechanisms.
The 1/n throughput limit was overcome in the prior art by a fully parallel-detection spectroscopic approach that uses a diffractive tilted etalon, called virtually-imaged phased array (VIPA) [9] in combination with an array-type detector such as a CCD camera. VIPA's spectral selection is given by the interference of multiple reflections at two optical flats yielding equivalent performances to FP in terms of resolution. However, the coating of the first surface is totally reflective with a narrow AR window to allow all the light to enter the interferometer. Besides minimizing losses, this conventional design avoids all the light being transmitted to the detector. As a consequence, with respect to an equivalent Fabry-Perot spectrometer, the signal strength is improved by a factor n equal to the finesse of the interferometer.
Thus, the spectral unit according to the exemplary embodiment of the present disclosure can include a spectrometer employing at least one virtually imaged phased array (VIPA) etalon 300 as shown in
For example, cascading two or more VIPAs can be done to increase contrast without hurting significantly the sensitivity of the spectrometer. A single VIPA etalon 300 generally produces spectral dispersion along one spatial direction, parallel to its coating direction, while leaving unchanged the optical beam propagation in the direction perpendicular to its coating direction. Multiple single VIPA etalons can be cascaded in such a way that each VIPA's orientation matches the spectral dispersion axis from the previous stage interferometer.
In a two-stage VIPA, the second etalon 310 can be placed orthogonally to the first one 300. The spectral pattern exiting the first stage can enter the second etalon through the input window. Both etalons disperse light, in orthogonal directions, so the overall spectral axis of the two-stage device lies along a diagonal direction, at 135° from the horizontal axis if the etalons have identical dispersive power. The second etalon 310 can separate Brillouin signal from crosstalk because, although spatially overlapped after the first stage, their frequencies at each spatial location are different. Thus, after the second stage, while the Brillouin spectrum lies on a diagonal axis, crosstalk components due to the limited extinctions of the etalons can be separated and mostly confined to the horizontal and vertical axis.
Besides spatial separation of signal and stray light, the exemplary two-stage spectrometer also facilitates the selective spectral filtering. An appropriate aperture mask 320 can be placed at the focal plane of the first VIPA 300, where a highly resolved spectral pattern is formed. Examples of the mask 320 are a slit or a rectangular aperture. This mask allows unwanted spectral components to be blocked and only the desired portion of the spectrum to pass to the second VIPA 310. For a preferred performance, it is possible to maintain only two Brillouin peaks (Stokes and anti-Stokes from two adjacent orders) and have a vertical mask cut off all elastic scattering peaks. This greatly reduces crosstalk in the second-stage VIPA 310 and assist in avoiding a saturation of the pixels in a CCD camera placed afterward, which are illuminated by strong unfiltered elastic scattering light.
This cross-axis cascade can be extended to a third stage. In a three-VIPA spectrometer, a third VIPA 330 is oriented perpendicular to the spectral axis of the preceding two stages, so that the Brillouin spectrum can enter through the input window of the VIPA 330. A second mask 340 can be employed to further reduce crosstalk. Due to the combined dispersion of the three etalons, the overall spectral axis can be further rotated, e.g., at about 170° if all the etalons have the same dispersive power.
Following the same cascading principle, a multiple VIPA interferometer of N stages can be provided. The k-th VIPA can be oriented at an appropriate angle to accept the spectrum dispersed through the preceding k−1 stages. The building block of each stage is modular, comprised of a cylindrical lens Ck, an etalon VIPAk, a spherical Fourier transform lens Skf with focal length fk, a mask and a spherical relay lens Sk,k+1 of focal length fk,k+1.
In the first stage, the VIPA is oriented along the direction v1 at an angle θ1 with respect to the horizontal axis (θ1=90° in our experiment), and its spectral dispersion direction, d1, can be parallel to v1 with ψ1=θ1. In the double-VIPA interferometer, the second etalon can be aligned along v2 at an angle θ2=ψ1±π/2, perpendicular to the spectral direction of the first stage d1 (θ2=180° in our experiment). After the two etalons, the spectrum can emerge at an angle ψ2 along spectral dispersion direction, d2. In a three-stage interferometer, the third VIPA must be oriented perpendicular to d2, at an angle θ3=ψ2±π/2. This exemplary arrangement can result in a final dispersion direction s3 at an angle ψ3.
For each stage, the exemplary dispersion angle, φk, imposed on a beam of wavelength λk by the kth VIPA interferometer was previously derived in both plane-wave and paraxial approximations. The focal length fk of the spherical lens, Skf, after the VIPA determines the linear dispersion power, sk, of the kth stage: sk=φk*fk. Since telescopes are used to link two subsequent VIPA stages, the overall linear dispersion also depends on the magnifications introduced by such optical systems. For example, each kth stage can introduce a magnification Mk=fk/fk−1, on the spectral pattern obtained by the previous k−1 stages, so that the effective linear dispersion, s′k, due to the kth stage is given by: s′k=sk*MN* . . . * Mk+1. Therefore, along the overall spectral axis, the total linear dispersion, SN, of an N-stage multiple VIPA interferometer is calculated to be: SN=sqrt(Σ1Ns′k2), which suggests a theoretical improvement in spectral resolution. When all the spectral dispersions are equal, i.e. s′1=s′2= . . . =s′N≡s, the total dispersion becomes SN=√{square root over (N)}·s.
Knowing the spectral dispersion and the optical magnification introduced by each stage, the overall dispersion axis can be computed and/or determined. It can be shown:
ψk+1−ψk=tan−1(s′k+1/Sk)→tan−1(1/√{square root over (k)}); (7)
θk+1−θk=tan−1(s′k/Sk−1)→tan−1(1/√{square root over (k−1)}). (8)
Here, the expressions marked with arrows apply in the case equal dispersion and unit magnification for all stages.
In terms of extinction, a single VIPA spectrometer has extinction, C proportional to its finesse squared: C˜4F2/π2, for an input beam with Airy profile. After N VIPA etalons of equal finesse F, the spectral extinction or contrast improves, in principle, to: C˜(4F2/π2)N.
It is possible to compare the extinction performance of single-stage, two-stage and three-stage VIPA spectrometers by coupling the single mode laser light directly into the spectrometer and placing a CCD camera in the focal planes of S1f, S2f, and S3f, respectively To overcome the limited dynamic range of the CCD, it is possible to record the spectrum at various optical power levels, with calibrated neutral density (ND) filters of optical densities in the range from 0 to 7. Subsequently, the full dynamic range spectrum was reconstructed by resealing the recorded raw spectra according to the respective attenuation levels. The single-stage VIPA shows an extinction level of about 30 dB over a wide frequency range between 5 and 25 GHz. The extinction can be improved to 55 dB with two stages and to nearly 80 dB with three VIPA etalons in the middle of the frequency range. It might be possible to improve the extinction up to 90 dB by minimization of aberrations in the optical system and improvements of beam shape or profile.
Besides the cross-axis cascading, another approach to improve the extinction ration of a VIPA etalon is to make the intensity profile of the VIPA output close to a Gaussian shape, a technique known as apodization.
In another exemplary variant of the apodization, the VIPA etalon 390 can be made with a gradient coating on the exit surface 394, such that its reflectivity and transmission is varied spatially, producing a rounded, Gaussian-like intensity profile 396, as shown in a right-hand portion of
At the output of the spectrometer light can be collected by a spatially-resolved photo-detector. Given the paucity of the expected Brillouin photons, high sensitivity solutions are required. Several technologies are known in the art for this purpose, CCD cameras, array of photodiodes, array of photomultiplier tubes of small active area, CMOS, or quadrant photodetectors. For some of these solutions, tight control or tracking over the long-term frequency drift of the laser and/or the spectrometer may be necessary.
For exemplary situations in which high spectral resolution and contrast are needed a scanning filter such as a Fabry-Perot interferometer can be used. The Fabry-Perot scanning interferometer can have a free spectral range of about 50 GHz, and finesse of about 1000; it can operate in single-pass configuration or in multipass, fixed or tandem, to enhance contrast. Alternatively, a fixed filter with a bandpass, notch, or edge type may be used, instead of scanning filter, to measure the magnitude of certain frequency component. In such case, the optical frequency of the pump wave can be stabilized or locked with respect to the fixed filter. The exemplary apparatus can employ different detection techniques, including interferometric heterodyne detection, or different illumination sources such as tunable lasers, or double-frequency pump sources.
In another exemplary embodiment of the present disclosure, other spectral modalities can be added in parallel to Brillouin spectroscopy such as Raman and fluorescence spectroscopy. This would allow rapid and complete fingerprint of ocular properties including, mechanical, chemical and functional information. Since fluorescence, Raman and Brillouin spectrum are in different region of the electromagnetic spectrum, after the collection optics, scattered and fluorescent light can undergo a coarse spectral dispersion to separate the fluorescence and Raman information about the sample under test. Preferred solutions for such coarse spectral dispersion are gratings, dichroic mirrors, interferometric bandpass filters. These exemplary solutions facilitate delivering Raman, Brillouin and fluorescent light to the same CCD camera. Different spatial regions of the CCD camera will map to a different spectral scale to allow simultaneous characterization in three different spectral regimes. Alternatively, the different signatures can be directed to different detecting devices operated in synch This multi-modality can be advantageous because the processes involved sample independent and diverse characteristics of a given material thus yielding mechanical as well as optical and chemical information about the analyzed sample.
To facilitate Brillouin imaging, the focused light can be scanned over the sample (eye).
To optimize the imaging speed and sampling resolution, additional beam scanning patterns can be evaluated: XY en face model 560, XZ modes 563, and XYZ volumetric imaging 565, as illustrated in
In this context, a single-stage VIPA spectrometer with apodization reaching ˜60 dB of extinction, can be configured to interrogate multiple spatial points simultaneously in the ocular tissues, thus greatly reducing the acquisition time of a Brillouin image.
Another exemplary method to interrogate multiple spatial points can include a use of a line focus 580. The Brillouin scattered light 582 produced from the line focus 580 can be relayed and coupled to the etalon 390. The spatio-spectral pattern 590 projected on the detection unit can be used to generate the information about the ocular tissues.
According to one exemplary embodiment of the present disclosure using which testing was performed, a light source, imaging optics, a spectrometer, and a computer can be utilized. The light sources for the two prototypes are a frequency-doubled diode-pumped Nd-YAG laser emitting a 532-nm wavelength with a line width of 1 MHz and a grating-stabilized single longitudinal mode laser diode emitting at 780 nm with a line width of about 100 MHz. Light can be focused on a sample through a 35 mm or a 11 mm focal length lens. In the prototypes, it is possible to use a beam scanner and a motorized translation stage to move the sample. It is possible to use the epi-detection configuration so that scattered light is collected through the same lens. A single mode fiber was used as confocal pinhole.
Brillouin-scattered light from the sample can be coupled into the VIPA spectrometer for high spatial separation of the spectral components in the plane of an Electron-Multiplied CCD camera. The spectrometer can employ a 3 nm bandpass filter to block fluorescence light. The optical design which can be used for the prototype can include a combination of a two-stage VIPA spectrometer and a variable attenuation neutral density filter. The spectrometer features a spectral resolution of about 1 GHz, an extinction ratio of about 75 dB, and a total insertion loss of 7 dB with a finesse of about 40.
To test the possibility of measuring the lens elasticity in vivo, Brillouin imaging was performed on laboratory mice 500 (C57BL/6 strain), as illustrated in
Using in vivo Brillouin microscopy, the natural age dependence of the lens modulus has been investigated. The peak Brillouin shift observed at the center of the lens nucleus in a mouse at 18 month old was 16 GHz, whereas the shift in a younger mouse at 1 month old was about 11.5 GHz. The study was extended to 12 mice of different ages to find an evident trend of age-related stiffening. Next, one mouse was imaged every week for two months and obtained consistent age-related data. The exemplary results indicate a quantitative (linear-log) relationship between the hypersonic elastic modulus and the animal age. This exemplary result indicates the first in vivo data evidencing an age-related stiffening of lenses in mice.
The feasibility of measuring the corneal elasticity was also tested. Brillouin imaging was performed on bovine eyes ex vivo.
Certain exemplary biomechanical diagnostic parameters can be extracted from Brillouin images. For example, the slope 730, defined as the rate of change of elastic modulus across depth; the mean modulus 734, defined as the spatial average of the elastic modulus; and the Brillouin stiffness 738 defined as the area under the Brillouin profile curve. These exemplary parameters can be defined for the whole cornea, or limited to specific regions (e.g. anterior, central and posterior).
On bovine corneas, the possible use of the exemplary Brillouin technique has been reviewed as a monitoring tool for corneal procedures. It was tested if Brillouin biomechanical imaging is sensitive to corneal stiffness changes induced by therapeutic procedures known as corneal collagen crosslinking (CXL). CXL is a technique that promotes the formation of covalent bonds between collagen fibers inside corneal stroma through a photosensitizing agent and light irradiation. Enhancing the amount of crosslinks between collagen fibers leads to a stiffer corneal stroma. The CXL procedure was performed on bovine corneal samples. A photosensitizer (Riboflavin) was diffused into the stroma of the cornea after removal of the epithelium and was activated by illuminating blue light.
Using Brillouin parameters, it is possible to quantify the effect of CXL procedure. Using analogous procedure to the one described previously, we obtained corneal axial profiles before and after the CXL. In particular, CXL resulted in a large increase (6-fold) in the downward slope of Brillouin frequency over depth in the stromal region.
The relevance of the extracted biomechanical information for clinical purposes can be important. While it may not be very difficult to find studies that relate mechanical properties to physiological or pathological conditions, the use of the exemplary Brillouin technique and its hypersonic elastic properties and that its sensitivity to measure relevant changes is important, which can be address, by analyzing a human ocular tissue.
As for the crystalline lens, the loss of accommodation power leading to presbyopia likely affects every person beyond the age of 50. Age-related stiffening of the crystalline lens is widely considered to play a primary role in this process. However, data are highly variable in terms of magnitude and rate of elasticity change with age as well as internal regional variation of elastic modulus. The uncertainty in the data has prevented a definitive understanding of the aging lens biomechanics and of the accommodation mechanism thus hindering the progress of lens-based treatments for presbyopia. According to an exemplary embodiment of the present disclosure, the axial profile of 14 fresh human lenses ex vivo between 21 and 67 yrs of age has been tested. The exemplary spatially-resolved data clearly indicate that the age-related changes in the spatial distribution of elasticity are responsible for the age-related stiffening linked to presbyopia onset.
As for the cornea, the relevance of Brillouin measurements for keratoconus and corneal ectasia should be understood. Traditional stress-strain tests have indicated that keratoconic corneas have lower mechanical modulus than normal corneas. The significantly altered collagen networks observed in keratoconic corneas could explain the reduced modulus. Keratoconic corneas also have a substantial decrease in the number of crosslinks between collagen fibers[38]. The exemplary Brillouin modulus can be sensitive to such structural changes in corneal collagen network (see
Exemplary Brillouin axial profiles of exemplary cornea samples have been measured from advanced keratoconus patients and five normal controls. As shown in the exemplary graphs of
Using the exemplary infrared prototype system, it was possible to obtain an early in vivo study of cornea and lens in a human subject. Standard topography and pachymetry did not show any noticeable abnormality before and after Brillouin scans. With low-coherence interferometry, the central corneal thickness was measured to be about 530 μm. For exemplary Brillouin scan of the eye of the volunteer, the beam focus was translated along the optical axis from the cornea to the lens at a speed of 50 μm/s (in air), and the Brillouin spectra were acquired with a CCD frame rate of about 2.5 Hz (i.e. integration time of about 0.4 s). Incident light power at the sample was measured to be about 0.7 mW, which is about 30 times lower than the safety limit.
Therefore, it has been confirmed that the features illustrated in animal studies and in human samples ex vivo are also measurable in vivo, and the instrument is sensitive to detect the elasticity of human cornea and lens.
Finally, the exemplary configuration has been tested beyond the use with the ocular tissue, to demonstrate its applicability and significance with turbid biological tissue, e.g. fresh porcine skin. For example, skin includes two major layers, e.g., epidermis and dermis. The skin tissue sample was cut into thin slices of epidermis and dermis and placed them horizontally on a glass-bottom dish plate with a few saline drops on top to prevent drying. For this measurement, 20 mW of incident light focused onto the tissue through the glass plate was used.
Degradation of skin elasticity is a common problem that can be caused by many factors including age, injuries and diseases. Brillouin confocal microscopy using the apparatus and methods according to the exemplary embodiments of the present disclosure can address skin elasticity below the superficial layer at the microscopic level. With low-power probe beams and rapid spectral readout, the exemplary instrument can serve as a non-invasive non-contact tool to monitor skin elasticity. From a mechanical standpoint, the viable epidermis, below the stratum corneum (˜10 microns), can be composed by living cells and is softer than the dermis, mostly made up by collagen and elastin fibers. In agreement with the micro-structural composition of the skin layers, exemplary Brillouin spectroscopy, as illustrated in the graphs of
A factor determining skin elasticity can be hydration. Dry skin can have higher Young's modulus of elasticity than hydrated skin by two to three orders of magnitude.
Regarding Presbyopia, the exemplary Brillouin imaging according to the exemplary embodiments of the present disclosure can facilitate establishing the age-dependent spatially-varying elastic properties most involved in the decline of accommodation power. This can provide mechanistic insights on the accommodation process. This may facilitate a narrow investigation of the underlying molecular processes connected with the accommodation function through the biomechanical information, to address/improve the effectiveness of current procedures for softening/refilling the lens, and to design, develop and test other approaches to preserve/restore accommodation inspired by the understanding of lens accommodation biomechanics.
With respect to Keratoconus and Crosslinking the exemplary Brillouin imaging according to the exemplary embodiments of the present disclosure can be used for both early identification of keratoconus and customization of treatment. For example, people suffering from keratoconus can present mechanical abnormalities before clinical symptoms are noted. In addition, the exemplary Brillouin imaging according to the exemplary embodiments of the present disclosure can detect localized weakness of the cornea, thus identifying regions that are more susceptible to the degenerative progression of this disease. This will facilitate extending the practice of crosslinking to early keratoconus individuals. This will also facilitate devising crosslinking procedures that are spatially targeted to the weak corneal regions as well as customized to the specific strengthening need of each patient in terms of photosensitizer and light dose to use. For example, the exemplary Brillouin imaging according to the exemplary embodiments of the present disclosure can be used to obtain a mechanical feedback mechanism that can guide the crosslinking procedures in real time. To this end, the quantitative biomechanical parameters described in this invention are poised to become standardized indexes to evaluate the biomechanical health of the cornea and treatment outcome.
Concerning LASIK, the exemplary Brillouin imaging according to the exemplary embodiments of the present disclosure and its extracted parameters can be used for pre-screening test to judge the compatibility of the patient cornea with LASIK surgery. Individuals with low Brillouin elasticity can be probably counseled to avoid the procedure. For example, once the cornea is mapped out biomechanically, patient currently excluded because of thin corneas or because of too high correction/ablation needed, might be considered safe on the basis of their Brillouin corneal strength. In addition, the ablation patterns can be modified so that the biomechanically strong portions of the cornea are left intact during the procedure. This may yield similar vision correction outcomes, with significantly reduced risk to destabilize the mechanical integrity of the cornea.
Turning to Glaucoma, applanation tonometry is generally used to monitor intraocular pressure for glaucoma screening and management. The test is based on Goldmann equations derived from the Imbirt-Fick law. It is well known in the art that this measurement needs to be corrected for corneal thickness, curvature and elasticity. Correction formulae for thickness and curvature has been described. The exemplary procedure according to the exemplary embodiments of the present disclosure can address a correction of the corneal elasticity component of the intraocular pressure results obtained with standard applanation tonometry. To this end, the protocol calls for previous systematic calibration to be performed combining the exemplary Brillouin microscope, an applanation tonometer and a syringe manometer. As first calibration, using ex vivo eyeballs, the IOP can be manipulated by inserting a syringe manometer into the anterior chamber of an intact eye ex vivo. The IOP measured by applanation tonometry can be compared with the “true” IOP controlled by the syringe manometer while monitoring the corneal elasticity from Brillouin microscopy as well as the corneal thickness and curvature. Several measurements should be done in different conditions of corneal elasticity, including: (1) eyes from different species such as rabbit, pig and cow, human; (2) eyes with stiffened corneas by collagen crosslinking; and (3) eyes with softened corneas by loosening the corneal collagen fibers using various chemicals. This can facilitate deriving a quantitative correction formula for applanation tonometry based on the novel corneal elasticity data from the Brillouin microscopy.
In light of these technological advancement and enabling measurements, the exemplary apparatus and method according to the exemplary embodiments of the present disclosure can facilitate ophthalmic research and patient care since they can provide non-invasive, non-contact and microscopic information on ocular biomechanics in situ. The Brillouin ocular analyzer can be proven to be a useful diagnostic tool, facilitating early diagnosis, screening of at-risk patients, monitoring therapeutic responses, developing novel approaches for treatment, and understanding pathogenesis.
The foregoing merely illustrates the principles of the disclosure. Various modifications and alterations to the described embodiments will be apparent to those skilled in the art in view of the teachings herein. Indeed, the arrangements, systems and methods according to the exemplary embodiments of the present disclosure can be used with and/or implement any OCT system, OFDI system, SD-OCT system or other imaging systems, and for example with those described in International Patent Application PCT/US2004/029148, filed Sep. 8, 2004 which published as International Patent Publication No. WO 2005/047813 on May 26, 2005, U.S. patent application Ser. No. 11/266,779, filed Nov. 2, 2005 which published as U.S. Patent Publication No. 2006/0093276 on May 4, 2006, and U.S. patent application Ser. No. 10/501,276, filed Jul. 9, 2004 which published as U.S. Patent Publication No. 2005/0018201 on Jan. 27, 2005, and U.S. Patent Publication No. 2002/0122246, published on May 9, 2002, the disclosures of which are incorporated by reference herein in their entireties. It will thus be appreciated that those skilled in the art will be able to devise numerous systems, arrangements, and procedures which, although not explicitly shown or described herein, embody the principles of the disclosure and can be thus within the spirit and scope of the disclosure. In addition, all publications and references referred to above can be incorporated herein by reference in their entireties. It should be understood that the exemplary procedures described herein can be stored on any computer accessible medium, including a hard drive, RAM, ROM, removable disks, CD-ROM, memory sticks, etc., and executed by a processing arrangement and/or computing arrangement which can be and/or include a hardware processors, microprocessor, mini, macro, mainframe, etc., including a plurality and/or combination thereof. In addition, certain terms used in the present disclosure, including the specification, drawings and claims thereof, can be used synonymously in certain instances, including, but not limited to, e.g., data and information. It should be understood that, while these words, and/or other words that can be synonymous to one another, can be used synonymously herein, that there can be instances when such words can be intended to not be used synonymously. Further, to the extent that the prior art knowledge has not been explicitly incorporated by reference herein above, it can be explicitly being incorporated herein in its entirety. All publications referenced above can be incorporated herein by reference in their entireties.
This application is based upon and claims the benefit of priority from U.S. Patent Application Ser. No. 61/480,885, filed Apr. 29, 2011, the entire disclosure of which is incorporated herein by reference.
Number | Date | Country | |
---|---|---|---|
61480885 | Apr 2011 | US |