SPECKLE REDUCTION IN OPTICAL COHERENCE TOMOGRAPHY BY PATH LENGTH ENCODED ANGULAR COMPOUNDING

Abstract

Accordingly, exemplary embodiments of an apparatus obtaining information associated with a structure can be provided. For example, with a first arrangement, it is possible to separate at least one first electro-magnetic radiation into a first radiation and a second radiation forwarded to the structure, with the first and second radiations having different path lengths. Using a second arrangement, it is possible to receive third and fourth radiations from the structure associated with the first and second radiations and a fifth radiation received from a reference. Further, with a third arrangement, it is possible to generate data corresponding to an amount of a ranging depth within the structure associated with the second arrangement. For example, a difference between the path lengths of the first and second radiations is equal or greater than the ranging depth.

Description

FIELD OF THE INVENTION

The present invention relates to optical imaging methods and apparatus useful in medical diagnosis and imaging. In particular, the present invention relates to a high speed method for implementing angular compounding, e.g., angular compounding by path length encoding (“ACPE”) for reducing speckle in optical coherence tomography images.

BACKGROUND INFORMATION

Optical coherence tomography (“OCT”) is a technique for obtaining high resolution cross-sectional images of biological tissues. Clinical OCT studies conducted in the gastrointestinal tract and cardiovascular system have shown that OCT is capable of providing images of the architectural (>20 μm) microanatomy of a variety of epithelial tissues, including the layered structure of squamous epithelium and arterial vessels. However, for certain medical applications, such as the early detection of high-grade dysplasia in Barrett's esophagus and the identification of inflammation within atherosclerotic plaques, visualization of structures that are on a size scale of <20 μm may be preferable. OCT systems, with typical axial resolutions ranging from 8-12 μm, have the potential to resolve many of these structures, including nuclei, individual glands, and macrophages. Unfortunately, speckle, which occurs on the same size scale as these features, may prohibit unambiguous including nuclei, individual glands, and macrophages. Unfortunately, speckle, which occurs on the same size scale as these features, may prohibit unambiguous identification of the cellular and subcellular tissue components required for widespread clinical utilization of such technology.

Catheter or endoscope access and high-speed imaging is used in order to perform OCT in the internal organs of patients. In order to minimize diameter, most catheter-based OCT probes employ a single optical fiber to illuminate the sample and detect the signal from the tissue. High frame rates (typically 4-10 frames per second) are preferred for performing OCT imaging while minimizing artifacts caused by patient motion. A way to reduce speckle in OCT images that does not significantly increase the complexity of single optical fiber probe designs while maintaining high frame rates may be beneficial for applying OCT to accurately detect and quantify key microscopic tissue structures in patients.

The reduction of speckle in the OCT images speckle has been previously described. A publication by J. M. Schmitt, “Array Detection for Speckle Reduction in Optical Coherence Microscopy,” Phys. Med. Biol. 42, 1427-1429, 1997, the entire disclosure of which is incorporated herein by reference, describes a procedure for a speckle reduction by averaging multiple images acquired at different angles, known as angular compounding. In this publication, multiple (N) detectors receive images that have been acquired from different angles. The images are averaged incoherently, providing an improvement (√{square root over (N)}) in the signal to noise ratio (“SNR”). While this technique has the advantage that the measurements may be performed in real-time, the experimental apparatus as therein would not be compatible with a single fiber optic catheter.

Another publication, M. Bashkansky and J. Reintjes, “Statistics and reduction of Speckle in Optical Coherence Tomography,” Opt. Lett. 25, 545-547, 2000, the entire disclosure of which is incorporated herein by reference, describes an alternative technique for angular compounding to reduce speckle. In this method, a retroreflector apparatus is translated in front of the objective lens to change the angle of the incident beam on the tissue. N successive images are acquired and added incoherently to reduce speckle, again improving the SNR by a factor of √{square root over (N)}. While this method may be less complex than the use of multiple detectors, the time needed thereby to acquire the images is increased by N. In addition, the implementation of this method within the confines of a small diameter catheter or endoscope could be difficult.

SUMMARY OF THE INVENTION

According to one exemplary embodiment of the present invention, a method for performing angular compounding to reduce speckle within OCT images is provided, which may use angular compounding by path length encoding (“ACPE”). With ACPE, a high-speed acquisition can be maintained, and modifications to standard OCT catheter optics are likely to be minimal. In another exemplary embodiment according to the present invention, an apparatus for imaging is provided that uses ACPE.

Accordingly, apparatus probe catheter, and method according to exemplary embodiments of the present invention are provided for irradiating a sample. In particular, an interferometer may forward forwarding an electromagnetic radiation. In addition, a sample arm may receive the electromagnetic radiation, and can include an arrangement which facilitates a production of at least two radiations from the electromagnetic radiation so as to irradiate the sample. Such arrangement can be configured to delay a first radiation of the at least two radiations with respect to a second radiation of the at least two radiations.

According to another exemplary embodiment of the present invention, a reference arm provides a further electromagnetic radiation. In particular, the interferometer receives the first, second and further radiations, and forms a resultant signal based on the first, second and further radiations. A processing arrangement may also be provided for generating a first image based on the first radiation and a second image based on the second radiation, such that the first and second images are different from one another. The further image may be generated based on the first and second images. In addition, the further image may have a noise that is smaller than a noise of the first image and a noise of the second image. Also, the further image may have a signal to noise ratio that is improved according to the equation:

${\mathrm{SNR}}_{\mathrm{ACPE}}=\frac{〈S_{\mathrm{OCT}}〉}{\sqrt{\mathrm{var}\left[S_{\mathrm{OCT}}\right]}}\propto \frac{\sum _{i=1}^Nu_i}{\sqrt{\sum _{i=1}^Nu_i^2}},$

where SNRACPE is the signal to noise ratio, SOCT is an amplitude of a high-pass filtered OCT signal, m is a thickness of the arrangement, ui is an amplitude of a demodulated OCT signal at a spatial location, and N=2m−1. Further, m=2 and N=3 images associated with the at least two radiations may be obtained. The further image may also be generated based on a mathematical combination of the first and second images.

According to yet another exemplary embodiment of the present invention, the sample can be irradiated by the first irradiation at a first angle, and by the second radiation at a second angle, such that the first and second angles different from one another. For example, the first and second angle may be different from one another based on the delay and at least one of a phase and a incident angle of each of the first and second radiations. It is also possible to utilize a detector which detects the first electromagnetic energy, and forwards the detected energy to the processing arrangement.

In still another exemplary embodiment of the present invention, the arrangement can include two sections, each being configured to delay a respective one of the first and second radiations. In particular, a delay of the first radiation is preferably greater than a delay of the second radiation. In addition, the delay of a path of the first radiation compared to a path of the second radiation is at least 500 μm in air. Further, the delay of a path of the first radiation compared to a path of the second radiation is at least 1 mm in air. The arrangement may have a refractive index of at least 1.5 or at least 3.0. The arrangement can include silicon and/or an anti-reflective coating on at least one surface thereof. Such arrangement can be an anti-reflection-coated BK 7 glass. The glass may have a thickness of from about 1.6 mm to about 7.7 mm, and a refractive index of from about 1.51 to about 3.5.

According to yet another exemplary embodiment of the present invention, an apparatus for imaging is provided. The apparatus includes a sample arm that receive an electromagnetic radiation. The sample arm includes an arrangement which facilitates a production of at least two radiations from the electromagnetic radiation so as to irradiate a sample. The arrangement is configured to delay a first radiation of the at least two radiations with respect to a second radiation of the at least two radiations. The apparatus also includes a device for receiving the first and second radiations from the sample arm and at least one third radiation from a reference arm, such that the first and second radiations interfere with the third radiation. Further, the apparatus includes at least one of spectral separating unit which separates spectrum of at least one of the first, second and third radiations into frequency components, and at least one detection arrangement including a plurality of detectors. Each detector is capable of detecting at least a portion of at least one of the frequency components.

In a further exemplary embodiment of the present invention, an apparatus is provided. Such apparatus includes at least one first arrangement that provides at least one first electro-magnetic radiation to a sample arm and at least one second electro-magnetic radiation to a non-reflective reference arm. A frequency of radiation provided by the first arrangement varies over time. The sample arm receives the first electromagnetic radiation, and includes an arrangement which facilitates a production of at least two radiations from the electromagnetic radiation so as to irradiate a sample. The arrangement is configured to delay a first radiation of the at least two radiations with respect to a second radiation of the at least two radiations. The apparatus also includes at least one second arrangement detecting an interference between the first and second radiations generated at the sample arm and the second electro-magnetic radiations generated at the reference.

In still another exemplary embodiment of the present invention, an apparatus is provided. The apparatus includes a at least one first arrangement providing at least one first electro-magnetic radiation to a sample arm and at least one second electro-magnetic radiation to a reference arm. The first and/or second electro-magnetic radiation have a spectrum which changes over time. The spectrum contains multiple frequencies at a particular time. The sample arm receives the first electromagnetic radiation, and includes an arrangement which facilitates a production of at least two radiations from the electromagnetic radiation so as to irradiate a sample. The arrangement is configured to delay a first radiation of the at least two radiations with respect to a second radiation of the at least two radiations. The apparatus also includes at least one second arrangement detecting an interference between the first and second radiations generated at the sample arm and the second electro-magnetic radiations generated at the reference.

Other features and advantages of the present invention will become apparent upon reading the following detailed description of embodiments of the invention, when taken in conjunction with the appended claims.

BRIEF DESCRIPTION OF THE DRAWINGS

The invention is illustrated in the drawings in which like reference characters designate the same or similar parts throughout the figures of which:

FIG. 1A is a schematic diagram of an exemplary embodiment of an ACPE-OCT apparatus according to the present invention.

FIG. 2A is a schematic diagram of an exemplary embodiment of an Intralipid-Agar phantom according to the present invention.

FIG. 2B is an exemplary OCT image of the phantom of FIG. 2A without the insertion of a BK7 ACPE element.

FIG. 2C is an illustration of exemplary ACPE OCT images obtained by splitting a sample beam in two parts using an exemplary 3.1 mm BK7 glass element.

FIG. 3A is an exemplary OCT image of a ventral forearm obtained in vivo prior to the insertion of the BK7 ACPE element.

FIG. 3B is an exemplary compounded ACPE OCT image of the skin acquired at the same location as that in FIG. 3A.

DETAILED DESCRIPTION

I. Materials and Methods

FIG. 1A shows an exemplary embodiment of an apparatus according to the present invention which uses ACPE. An optical glass 20 is placed in an imaging path 10 of a conventional OCT imaging apparatus 5, splitting the incident field into two beamlets 1 and 2 which are provided on beam paths 30 and 40, respectively. For the purpose of the present discussion, a beamlet can be defined as a portion of a beam. An optical element (e.g., the optical glass 20) causes a portion of the incident beam (beamlet 2) to experience a greater path length delay than beamlet 1. In addition, beamlet 2 illuminates the sample at a different angle than beamlet 1. A lens may be provided to converge the beamlets 1 and 2 onto the sample (e.g., a tissue). As a result, multiple OCT images of a specimen in or on the sample, each acquired at a different angle, can appear simultaneously on the OCT display (see FIG. 1B). In particular, FIG. 1B shows that the top image (i.e., labeled as 1+1) corresponds to the image fainted from path 1 (i.e., incident and reflected), the middle image (i.e., labeled as 1+2, 2+1) corresponds to the image formed from path 1 incident, path 2 reflected and path 2 incident, path 1 reflected. Further, the bottom image (i.e., 2+2) of FIG. 1B corresponds to the image formed from path 2 (incident and reflected).

These beamlets 1 and 2, thus encoded by optical path length, also illuminate the sample at different angles. As a result, multiple OCT images, each acquired at different angles, are preferably present in a single OCT frame (as shown in FIG. 1B). For example, when the optical element contains m distinct thicknesses of glass, with each optical thickness a multiple of the others, 2m−1 OCT images can be obtained in one OCT frame. Each image may be separated by a group delay of D(n−1)/2, where D is the thickness and n is the refractive index of the optical material. The distinct OCT images are then averaged to produce a composite OCT image with significantly reduced speckle. Since all of the images are acquired in one OCT frame, single frame acquisition time can be maintained. In addition, the modifications to the OCT probe may involve only the insertion of at least one small optical element in the beam path of the distal optics. These advantages of ACPE allow speckle averaging to be performed within the confines of a small diameter catheter or endoscope, and without compromising the acquisition speed.

Previous attempts of utilizing OCT speckle reduction by image compounding has shown that the addition of N images of the same intensity provides an SNR increase by a factor of √{square root over (N)}. For example, the Bashkansky publication described above describes that the speckle distribution in OCT takes the form of the probability density function:

$\begin{array}{cc}p\left(S_{\mathrm{OCT}}\right)=\frac{1}{\sqrt{2\pi }\kappa }\exp \left(\frac{-S_{\mathrm{OCT}}^2}{2{\kappa }^2}\right)& \left(1\right)\end{array}$

where κ=2A_Rσ, with A_R is the amplitude of the reference field, S_OCT is the amplitude of the high-pass filtered OCT signal, and σ its standard deviation. For this probability density function, it can be demonstrated that the SNR improvement obtained by averaging N images of the same amplitude is also a factor of √{square root over (N)}, a result that has been experimentally described in the Schmitt and Bashkansky publications. For ACPE, the distinct OCT subimages generally may not have equal amplitudes. However, with the assumption of isotropic backscattering, these subimages are related to the original OCT image, S0_OCT, by β/m², where β is the number of path length combinations that contribute to a distinct subimage. As a result, the SNR for ACPE may be defined as

$\begin{array}{cc}{\mathrm{SNR}}_{\mathrm{ACPE}}=\frac{〈S_{\mathrm{OCT}}〉}{\sqrt{\mathrm{var}\left[S_{\mathrm{OCT}}\right]}}\propto \frac{\sum _{i=1}^Nu_i}{\sqrt{\sum _{i=1}^Nu_i^2}},& \left(2\right)\end{array}$

where u_i is the amplitude of the demodulated OCT signal at a spatial location, and N=2m−1. In the case of m=2, N=3 images are obtained and the relationships between the amplitudes of the ACPE OCT subimages are S1_OCT=S3_OCT=¼S0_OCT and S2_OCT=2S1_OCT=½S0_OCT. The potential SNR improvement of the compounded m=2 ACPE image can then become SNR_ACPE/SNR₀=1.63, where SNR₀ is the signal to noise ratio of S0_OCT.

EXAMPLES

Example 1

An exemplary polarization-diverse OCT system according to the present invention can be utilized for the exemplary apparatus (e.g., as shown in FIG. 1A) that implements ACPE. Using such exemplary system according to the present invention, OCT images may be acquired at 2 frames per second (e.g., 500 axial pixels×500 transverse pixels), which may be displayed with an inverse gray-scale lookup table, and digitally archived. The optical source which can be used in such exemplary OCT system may hive a center wavelength of 1310 nm and a bandwidth of 70 nm, and thus providing an axial resolution of approximately 8 μm in the tissue.

A modified hand-held galvanometer probe (e.g., the apparatus of FIG. 1A) can be inserted in the sample arm of the OCT system 5. The objective lens 50 may have a focal length of 25 mm and a numerical aperture (NA) of 0.11, providing a measured 1/e² focal spot diameter of 23 μm. A square, antireflection-coated D=3.1 nun BK7 glass (n=1.51) (e.g., the optical glass 20, may be inserted between the optical fiber collimator and the objective lens 50 (see FIG. 1A) so that overlaps with half of the illuminating beam. With this exemplary configuration, an OCT image separation of approximately 800 μm can be achieved. When the glass plate is inserted, the spot diameter perpendicular to the glass edge may increase by a factor of two (46 μm). In the plane of the OCT image, however, the transverse resolution is likely preserved.

Example 2

A solid phantom consisting of 1% Intralipid solution and Agar can be used to measure a reduction in speckle provided by ACPE. For example, four (4) hairs can be embedded in the Intralipid-Agar gel at different transverse positions and depths. A schematic of the phantom is depicted in FIG. 2A. The corresponding OCT images with and without the BK7 glass plate are shown in FIGS. 2B and 2C, respectively. ROI's labeled 1-5 represent exemplary locations where the SNR improvement by ACPE is preferably measured. The insertion of the BK7 glass plate 20 (with a 3.1 mm BK7 glass element) in the sample arm can produce three copies of the original OCT image (see FIG. 2C), with each image being acquired at a different illumination angles and separated by group delay increments of, e.g., 800 μm. The amplitudes of the signals in the top and bottom images of FIG. 2C can be approximately half of the center image (S1_OCT/S2_OCT=S3_OCT/S2_OCT=1:2). The compounded ACPE image shown in FIG. 2D can best obtained by incoherently averaging the three images that are shown in FIG. 2C. A substantial reduction of speckle in the compounded image can be visualized in FIG. 2D. Compared to the original OCT image, the average ACPE SNR improvement for the five regions depicted in FIG. 2B is preferably 1.54±0.12 (mean±standard deviation).

Example 3

In order to demonstrate SNR improvement in vivo, ACPE OCT imaging can be performed on a ventral forearm of a subject. FIGS. 3A and 3B shows one representative set of images. FIG. 3A shows an exemplary OCT image of the ventral forearm obtained in vivo prior to the insertion of the BK7 ACPE element. Visual assessment of these pictures demonstrates a qualitative improvement in the compounded ACPE image (FIG. 3B). The boundary between the epidermis (E) and dermis (D) is more clearly demarcated with ACPE as shown in FIG. 3B. In addition, horizontal structures consistent with dermal vasculature are more readily identified in the ACPE image of FIG. 3B. The SNR can be measured for the ACPE and original OCT images in FIG. 3B, providing an exemplary SNR improvement of 1.56.

The apparatus and method according to the exemplary embodiments of the present invention are provide to reduce speckle in OCT images that does not decrease the OCT frame rate and requires only minor modifications to the OCT probe. The implementation of ACPE, implements potential compromises between speckle reduction and three other OCT system parameters: a) sample arm transverse resolution, b) total reference arm path length, and c) OCT image sensitivity. In particular, for any given objective lens, ACPE may compromise the transverse resolution in one dimension by underfilling the lens aperture for each individual beamlet. In most cases, increasing the numerical aperture (NA) of the objective can compensate for this resolution loss.

Speckle averaged ACPE-OCT images can be obtained at the same rate as conventional OCT images by scanning an increased reference arm path length delay at substantially the same frequency. In order to acquire each of the individual OCT subimages, the new scan length of the ACPE-OCT system can preferably be L(2m−1), where L is the original scan length of the OCT system. Using phase control RSOD lines, scan ranges up to 10 mm are possible, thus enabling m=3, L=2 mm, and a maximum predicted SNR improvement of ˜2.1.

Increasing the reference arm path length scan range while maintaining the scan rate, may increase the electronic bandwidth and decrease the sensitivity of the OCT system. Also, because ACPE splits the sample arm power into 2m−1 subimages, each subimage contains a fraction of the original sample arm power. When imaging a human tissue, these losses may primarily affect the penetration depth of the OCT image. Since many features of clinical relevance, such as-nuclei in patients with Barrett's esophagus or macrophages in atherosclerotic plaques, may bee present at tissue surfaces, for modest m, the improvements in image quality provided by ACPE likely outweigh sensitivity losses. Moreover, ongoing technical developments towards more efficient interferometer designs and higher power, clinically viable OCT sources may render ACPE sensitivity losses a non-issue.

The thickness of the path length encoding optical element described herein in the above-referenced Examples (3.1 mm BK7) may not be sufficient for OCT imaging in some tissues since the provide 800 μm separation between individual subimages. Increasing the thickness of the BK7 glass to 7.7 nun may allow a path length separation of 2 mm. This thickness may be adequate for a free-space, hand-held OCT probes, but can be problematic in small diameter, flexible catheters, where minimizing the rigid length can be important. In order to increase the optical thickness of the path length encoding element, a higher refractive index material such as silicon (n=3.5) may be used. To create a 2 mm delay with silicon, preferably approximately 1.6 mm of the material can be used. When using high refractive index glass, dispersion imbalances between the reference and sample arms should be considered. For high-resolution OCT imaging (Δλ/λ>10%), appropriate selection of the optical material used for path length encoding will depend on the center wavelength and bandwidth of the source.

Advantages of the exemplary embodiment of the method and apparatus according to the present invention may include the fact that the OCT frame rate is not increased, and the addition of only a single passive element in the OCT probe should be implemented. These features of ACPE make the system and method of the present invention compatible with OCT imaging in internal organ systems in patients. While implementation of ACPE may facilitate tradeoffs between speckle reduction and system sensitivity, the problems caused by speckle noise are likely more significant for clinical diagnosis than the penetration depth of modem OCT systems, especially at 1300 nm. Since difficulties in interpreting features on the size scale of 20 μm or less is in part a result of speckle noise in OCT images, ACPE would likely significantly improve the capabilities of OCT for the diagnosis of important diseases such as, but not limited to, dysplasia and inflammation in atherosclerosis.

Although only particular exemplary embodiments of the present invention have been described in detail herein above, those skilled in the art will readily appreciate and understand that many modifications are possible in the exemplary embodiments without materially departing from the novel teachings and advantages of this invention. Accordingly, all such modifications are intended to be included within the scope of the present invention as defined in the following claims.

Claims

1-26. (canceled)

27. An apparatus for obtaining information associated with a structure, comprising; a first arrangement configured to provide a first radiation and a second radiation that are forwarded to the structure, wherein the first and second radiations have different path lengths;a second arrangement configured to receive third and fourth radiations from the structure associated with the first and second radiations and a fifth radiation received from a reference; anda third arrangement configured to generate the information corresponding to an amount of a ranging depth within the structure associated with the second arrangement.

28. The apparatus according to claim 27, wherein the first arrangement includes an optical section with a refractive index of at least 1.5, the optical section being structured to propagate the first and second radiations.

29. The apparatus according to claim 27, wherein the first arrangement includes a section which has silicon that Is structured to propagate the first and second radiations.

30. The apparatus according to claim 27, wherein the first arrangement includes a section which has an anti-reflective coating on at least one surface thereof, the section being structured to propagate the first and second radiations.

31. The apparatus according to claim 27, wherein the third arrangement is further configured to generate additional data corresponding to an amount of the ranging depth within the structure associated with the first arrangement, and wherein the ranging depth is a penetration depth in the structure which includes an anatomical structure.

32. The apparatus according to claim 27, wherein the third arrangement is further configured to generate: second data associated with the first radiation which includes a first speckle pattern and second data associated with the second radiation which includes a second speckle pattern, andthird data associated with the structure based on the first and second data that has a speckle pattern the amount of which is smaller than the amount of at least one of the first speckle pattern or the second speckle pattern, andwherein the data, the second data and the third data are images, the third image data having a noise that is smaller than a noise of the first data and a noise of the second data.

33. The apparatus according claim 27, wherein the third arrangement is configured to generate the third image data based on a mathematical combination of the first and second image data.

34. The apparatus according to claim 27, further comprising a fourth arrangement which is configured to irradiate the structure by the first irradiation at a first angle, and by the second radiation at a second angle, the first and second angles being different from one another.

35. The apparatus according to claim 34, wherein the first and second angle are different from one another as a function of a delay between the first and second radiations and at least one of a phase or an incident angle of each of the first and second radiations.

36. The apparatus according to claim 27, further comprising a detection arrangement which is configured to detect the first radiation, and forward the detected radiation to the second arrangement.

37. The apparatus according to claim 27, further comprising a further arrangement which includes two sections, each section being configured to delay a respective one of the first and second radiations, and wherein a delay of the first radiation is greater than a delay of the second radiation.

38. The apparatus according to claim 27, further comprising a further arrangement configured to receive a plurality of radiation, irradiate the structure. and delay at least one radiation of the first and second radiations with respect to one another.

39. The apparatus according to claim 28, wherein the refractive index of the optical section is at least 3.0.

40. The apparatus according to claim 39, wherein at least one of the first section or the second section comprises an anti-reflection-coated BK 7 glass.

41. The apparatus according to claim 40, wherein the glass has a thickness of from about 1.6 mm to about 7.7 mm.

42. The apparatus according to claim 28, wherein at least one of the first section or the second section comprises an anti-reflection-coated glass that has a refractive index of from about 1.51 to about 3.5.

43. The apparatus according to claim 27, wherein a difference between the path lengths of the first and second radiations is equal or greater than the ranging depth.

44. A method for obtaining information associated with a structure, comprising; providing a first radiation and a second radiation to the structure, wherein the first and second radiations have different path lengths;receiving third and fourth radiations from the structure associated with the first and second radiations and a fifth radiation received from a reference; andgenerating the information corresponding to an amount of a ranging depth within the structure associated with the second arrangement.

45. The method according to claim 44, wherein a difference between the path lengths of the first and second radiations is equal or greater than the ranging depth.