1. Field of the Invention
The present invention relates to particularly an optical coherence tomography apparatus including an interference optical system which is used in the medical field, an optical coherence tomography method, an ophthalmic apparatus, a method of controlling the ophthalmic apparatus, and a storage medium.
2. Description of the Related Art
Currently, various types of ophthalmic apparatuses using optical devices are used. Such apparatuses include, for example, an anterior ocular segment imaging apparatus, a fundus camera, and a scanning laser ophthalmoscope (SLO). Among them all, an optical coherence tomography (OCT) apparatus (to be referred to as an “OCT apparatus” hereinafter) is an apparatus capable of obtaining a high-resolution tomogram of an object to be examined. This OCT apparatus has been becoming an indispensable apparatus for dedicated retinal outpatient clinics.
For example, the OCT apparatus disclosed in Japanese Patent Laid-Open No. 11-325849 uses low-coherent light as a light source. Light from the light source is split into measurement light and reference light through a splitting optical path such as a beam splitter. Measurement light is light to irradiate an object to be examined such as the eye through a measurement light path. Return light of this light is guided to a detection position through a detection light path. Note that return light is reflected light or scattered light containing information associated with an interface relative to the irradiation direction of light on the object. On the other hand, reference light is light to be guided to the detection position through a reference light path by being reflected by a reference mirror or the like. It is possible to obtain a tomogram of an object to be examined by causing interference between this return light and reference light, collectively acquiring wavelength spectra by using a spectrometer or the like, and performing Fourier transform of the acquired spectra. An OCT apparatus which collectively measures wavelength spectra is generally called a spectral domain OCT apparatus (SD-OCT apparatus).
In an SD-OCT apparatus, a measurement depth Lmax is represented, as an optical distance Lmax, by a pixel count N of the image sensor of a spectrometer and a spectrum width ΔK of the frequency detected by the spectrometer according to equation (1). Note that the spectrum width ΔK is represented by a maximum wavelength λmax and a minimum wavelength λmin. The pixel count N is often an even number, and is generally the factorial of 2, that is 1024 or 2048.
If, for example, a central wavelength of 840 nm, a band of 50 nm, and a pixel count of 1024 are set, λmax=840+50/2=840+25=865 nm, λmin=840−50/2=840−25=815 nm, and N=1024. In this case, optical distance Lmax=3.6 mm. That is, it is possible to perform measurement up to about 3.6 mm on the plus side relative to the coherence gate. The coherence gate is the point at which a reference light path coincides with an optical distance in a measurement light path. When a desired region (a distance in the depth direction) is sufficiently smaller than 3.6 mm (for example, 1 mm or less), the measurement depth can be reduced by decreasing the pixel count of the spectrometer. Decreasing the pixel count is important in order to speed up processing and reduce the data amount. This is because, when measuring a three-dimensional image of the retina, it takes much measurement time and produces a large amount of data. When an object to be examined is a moving object like the eye, in particular, it is required to further shorten the measurement time.
On the other hand, changing the pixel count of a spectrometer is equivalent to changing the resolution of the spectrometer. A problem in this case will be described with reference to
As is obvious from
In consideration of the above problems, the present invention provides a technique of correcting the contrast differences between images which are caused when wavelength resolutions differ (spectrometers differ in resolution in the case of an SD-OCT) in an FD-OCT apparatus such as an SD-OCT apparatus.
According to one aspect of the present invention, there is provided an optical coherence tomography apparatus which acquires a tomogram of a target object based on a light intensity detected for each wavelength by combining return light of measurement light from the target object with reference light corresponding to the measurement light, the apparatus comprising: a normalization unit adapted to normalize the light intensity detected for the each wavelength based on a transfer function corresponding to the wavelength resolution; and an image formation unit adapted to form a tomogram of the target object from the light intensities normalized by the normalization unit.
According to one aspect of the present invention, there is provided an ophthalmic apparatus which acquires a tomogram of an eye to be examined based on a light intensity detected for each wavelength by combining return light of measurement light from the eye with reference light corresponding to the measurement light, from a mirror, the apparatus comprising: a normalization unit adapted to normalize the light intensity detected for the each wavelength based on an attenuation function corresponding to a wavelength resolution; and an image formation unit adapted to form a tomogram of the eye from the light intensities normalized by the normalization unit.
According to one aspect of the present invention, there is provided an optical coherence tomography method of acquiring a tomogram of a target object based on a light intensity detected for each wavelength by combining return light of measurement light from the target object with reference light corresponding to the measurement light the method comprising: normalizing the light intensity detected for the each wavelength based on a transfer function corresponding to the wavelength resolution; and forming a tomogram of the target object from the light intensities normalized in normalizing.
According to one aspect of the present invention, there is provided a method of controlling an ophthalmic apparatus which acquires a tomogram of an eye to be examined based on a light intensity detected for each wavelength by combining return light of measurement light from the eye with reference light corresponding to the measurement light, the method comprising: normalizing the light intensity detected for the each wavelength based on an attenuation function corresponding to the wavelength resolution; and forming a tomogram of the eye from the light intensities normalized in normalizing.
Further features of the present invention will be apparent from the following description of exemplary embodiments with reference to the attached drawings.
An exemplary embodiment(s) of the present invention will now be described in detail with reference to the drawings. It should be noted that the relative arrangement of the components, the numerical expressions and numerical values set forth in these embodiments do not limit the scope of the present invention unless it is specifically stated otherwise.
The arrangement of an optical coherence tomography apparatus using an optical coherence tomography method according to the present invention will be described first with reference to
<Arrangement of Optical System>
The arrangement of an optical coherence tomography apparatus (OCT apparatus) will be described first. The OCT apparatus includes a light source 201, a lens 202, a beam splitter 203, an XY scanner 204, an objective lens 205, a scan lens 206, an imaging lens 207, a spectrometer 208, a reference mirror 209, a mirror driving mechanism 210, a focus driving mechanism 211, and a computer 219. The spectrometer 208 includes a prism 216, a lens 217, and an image sensor 218.
The beam splitter 203 splits the light emitted from the light source 201 into measurement light 212 and reference light 215 through the lens 202. The measurement light 212 reaches an eye 221 (target object) as an object to be examined through the XY scanner 204, the objective lens 205, and the scan lens 206. Return light 213 reflected by a retina 220 of the eye 221 sequentially returns through the scan lens 206, the objective lens 205, the XY scanner 204, and the beam splitter 203. This light reaches the spectrometer 208 through the imaging lens 207. On the other hand, the reference light 215 is reflected by the reference mirror 209. The beam splitter 203 combines the reference light 215 with the return light 213. Note that it is possible to adjust the optical length by causing the mirror driving mechanism 210 to control the reference mirror 209. In addition, the focus driving mechanism 211 allows to control the focus of the measurement light 212.
The light source 201 is an SLD (Super Luminescent Diode) as a typical low-coherent light source. Light from this light source has, for example, a wavelength of 840 nm and a band of 50 nm. Note that the band is an important parameter because it influences the resolution of an obtained tomogram in the optical axis direction. Although an SLD is selected as the light source 201, it is possible to use an ASE (Amplified Spontaneous Emission) source or the like as long as it can emit low-coherent light. Obviously, it is possible to use another kind of light source such as a halogen lamp depending on the type of object to be examined. Note, however, that since a wavelength influences the resolution of an obtained tomogram in the transverse direction, it is good to use a short wavelength, when placing importance on the resolution in the transverse direction.
The spectrometer 208 includes the prism 216, the lens 217, and the image sensor 218 (for example, a line sensor). The spectrometer 208 performs spectroscopy of the measurement light 212 to acquire spectrum data. The image sensor 218 can be a CMOS sensor capable of arbitrarily setting a pixel count N. When decreasing the pixel count N, as the band narrows, the depth resolution of the OCT decreases. For this reason, the band is kept unchanged, and the sampling intervals are increased. Increasing the sampling intervals is equivalent to decreasing the resolution of the spectrometer 208. It is possible to decrease the pixel count N of the image sensor 218 by thinning processing, binning, or the like. It is also possible to change the pixel count N to be imaged by changing the imaging magnification of the lens 217 and changing the spot diameter.
The computer 219 analyzes the wavelength spectrum data captured by the image sensor 218. The computer 219 includes a CPU, a RAM, and a ROM, and has, for example, functions of storing data, displaying images, and issuing a command to perform measurement as well as an analysis function. It is possible to obtain a tomogram of an object to be examined by making the XY scanner 204 raster-scan the measurement light 212 on the object in a direction perpendicular to the optical axis under the control of the computer. Note that an X-Y plane is a plane perpendicular to the optical axis, and Z-axis is an axis parallel to the optical axis. The Z-axis corresponds to an axis in the depth direction of a tomogram.
In this scheme, the mirror driving mechanism 210 modulates the phase of the reference light 215. It is therefore good to allow the mirror driving mechanism 210 to perform control with a precision equal to or higher than 1/100 of a central wavelength. If, for example, the central wavelength is 840 nm, the precision can be about 5 nm. In addition, in order to implement high-speed control, it is possible to combine a coarse motion stage capable of micrometer-order control with a fine motion stage capable of nanometer-order control.
The optical coherence tomography apparatus of the present invention splits light from the light source 201 into the measurement light 212 and the reference light 215 through the splitting optical path. The measurement light 212 irradiates an object to be examined through the measurement light path. In addition, the return light 213 of the measurement light 212 from the object is guided to the detection position through the detection light path. The focus driving mechanism 211 can control the focus of the measurement light 212. On the other hand, the reference light 215 is guided to the detection position through the reference light path. The reference mirror 209 is placed on the reference light path, and it is possible to adjust the coherence gate by controlling the reference mirror 209 using the mirror driving mechanism 210. The light guided to the detection position is decomposed into a wavelength spectrum, which is then analyzed.
A procedure for optical tomography processing according to the present invention will be described with reference to
In step S301, this apparatus starts measurement processing. In step S302, measurement conditions including the resolution and measurement range of the spectrometer 208 are selected. If, for example, there is a plurality of spectrometers having different resolutions, the apparatus may be configured to be able to select one spectrometer (as the spectrometer 208). In addition, the apparatus may be configured to be able to select one image sensor (as the image sensor 218) from a plurality of image sensors with different pixel counts. Furthermore, the apparatus may be configured to be able to move the lens 217 in the optical axis direction so as to change the range in which an image is formed on the image sensor 218, by changing the position of the lens (using a changing unit (not show) such as a stage mounted on the lens). Assume that in this case, a transfer function corresponding to the resolution of each spectrometer (the pixel count of each image sensor or the position of each lens) is set in advance. This transfer function may be an experimentally obtained attenuation function (a function corresponding to the attenuation of light intensity from the coherence gate in the depth direction), as shown in
In step S304, this apparatus calculates a tomogram from the signal acquired in step S303. In step S305, the apparatus normalizes a tomogram by using a transfer function corresponding to the resolution of the spectrometer. That is, the apparatus normalizes the light intensity in a measurement region by using this transfer function. For example, normalization is performed by dividing the light intensity in the measurement region by the transfer function.
In step S306, it is determined whether one or two or more measurement regions are all measured. If it is determined that all the measurement regions are measured (YES in step S306), the process advances to step S307. If it is determined that not all the measurement regions are measured (NO in step S306), the process returns to step S303.
In step S307, the apparatus connects the tomograms normalized in step S305 to form an overall image.
In step S308, after the coherence gate and the focus are returned to the initial positions, the measurement processing is terminated.
<Phase Modulation Interference Method>
In this case, a tomographic measurement method based on a phase modulation interference method of detecting the phase difference between the reference light 215 and the measurement light 212 will be described below. First of all, the light intensity from the position (x, y) of an object to be examined is represented by I(λ, δ) using a wavelength λ and a phase δ determined by the position of the coherence gate. Assume that in the case of interference, the light intensity I(λ, δ) from the object is represented by three variables including an incoherent component I0, a coherent component I1, and a phase φ(λ) of the component according to equation (3).
I(λ,δ)=I0(λ)+I1(λ)cos {φ(λ)+δ} (3)
In this case, δ can be changed by moving the position of the coherence gate. For example, δ=n is equivalent to moving the coherence gate by λm/4 when the central wavelength is λm, in consideration of round-trip optical paths.
The 3 frame method changes δ, a total of three times, to 0, π/2, and π, and measures light intensity at each time. This makes it possible to obtain three variables, that is, the incoherent component I0, the coherent component I1, and the phase φ(λ) of the component by
The 4 frame method changes δ, a total of four times, from 0 to 3π/2, π/2 at a time, and measures light intensity at each time. As in the above case, this makes it possible to obtain three variables, that is, the incoherent component I0, the coherent component I1, and the phase φ(λ) of the component by
The 5 frame method changes δ, a total of five times, from 0 to 2π, π/2 at a time, and measures light intensity at each time. As in the above case, this makes it possible to obtain three variables, that is, the incoherent component I0, the coherent component II, and the phase φ(λ) of the component by
Transforming the light intensities and phases corresponding to the wavelengths, obtained above, into light intensities and phases corresponding to wave numbers k=1/λ can obtain a tomogram F(z) according to equation (7). Assume that in transformation to wave numbers, resampling and the like are performed at equal intervals relative to the wave numbers in consideration of Fourier transform.
F(z)=FT−1{I1(k)exp(jφ(k))} (7)
where FT−1 represents inverse Fourier transform, and j represents an imaginary unit. In general, as the number of frames increases, noise is reduced owing to an averaging effect. On the other hand, as the number of times of measurement increases, the measurement time increases. That is, they have a trade-off relationship. Note that the phase modulation method described above is effective when a coherence gate is placed inside an object to be examined, and a mirror image is generated. That is, when a coherence gate is placed outside the object, since no mirror image is generated, it is possible to obtain the tomogram F(z) by performing Fourier transform of I(λ, δ) as indicated by equation (7). In the second embodiment, a measurement region is divided into a plurality of regions, and a coherence gate is placed at the boundary between the regions. For this reason, a mirror image is generated, and hence the phase modulation method is required. Obviously, the method to be used to remove a mirror image may be a method other than the phase modulation method.
<Image Formation>
Measurement to be performed when the transverse resolution is high will be described next with reference to
In step S401, the apparatus starts measurement processing.
In step S402, the operator selects measurement conditions. The measurement conditions include a depth resolution, a transverse resolution, a measurement range, and the resolution of a spectrometer. In addition, when performing measurement upon dividing the measurement range into one or two or more measurement regions, the selection of measurement conditions may include designation of a specific number by which the measurement region is divided. A depth resolution is determined by the band and focal depth of a light source. A transverse resolution is determined by a beam diameter and the like. Since the higher the transverse resolution, the shallower the focal depth, it is necessary to measure the measurement range upon dividing it into several measurement regions. The measurement time is determined by a measurement range, a resolution, the number of frames in the phase modulation interference method, and the like.
In this case, each region is represented by ZM(i), and the data of a tomogram in each region is represented by SM(i, k). In this case, M is the number of measurement regions, and i is a measurement region number, which ranges from 0 to M−1. If, therefore, the number of measurement regions is 1, M is 1. If the number of measurement regions is 4, M is 4. In addition, k is an element, which is an integer satisfying 0 to N−1. N is the pixel count of an image sensor.
In step S403, the apparatus acquires a signal. In order to measure the desired measurement region ZM(i), the apparatus adjusts a coherence gate and a focus. In this case, the apparatus starts measurement from ZM(0). At the position of the measurement region ZM(0), since the object to be examined is the eye, the coherence gate is set on the cornea side relative to the retina. Note that the apparatus controls the coherence gate to make the position of SM(i−1, nM) of each measurement region ZM(i) coincide with the position of SM(i, 0). In addition, the apparatus performs measurement corresponding to the required number of frames, as described in association with the phase modulation interference method, for each measurement region ZM(i). In this case, nm can be a value satisfying the origin to ½ of the attenuation function represented by equation (2) as well as being smaller than N/2.
In step S404, the apparatus calculates a tomogram from data corresponding to each frame count. Tomograms in the regions ZM(i) range from SM(i, 1) to SM(i, nM). Note that the data of Sm(0, 0) represents the position of the coherence gate, but does not originate from a slice. For this reason, the data of SM(0, 0) may be deleted.
In step S405, the apparatus normalizes the light intensity in each measurement region by using a transfer function predetermined for each measurement condition such as the resolution of the spectrometer. Letting DM(i, k) be the data of a transfer function used for normalization, data HM(i, k) after normalization is expressed as equation (8). This normalization processing, that is, normalization using a transfer function predetermined for each spectrometer resolution, makes it possible to correct the contrast differences between images generated when spectrometers have different resolutions.
HM(i,k)=SM(i,k)/DM(i,k) (8)
A transfer function may be a value theoretically calculated as indicated by equation (2) or a value held in advance upon being experimentally obtained as shown in
In step S406, it is determined whether processing for all the regions is complete. When acquiring a two-dimensional image, for example, the apparatus acquires 512 tomograms by scanning with an X scanner (Y scanner). When acquiring a three-dimensional image, the apparatus acquires 512×512 tomograms by further scanning 512 lines with the Y scanner (X scanner). When acquiring a three-dimensional image while dividing a measuring region into four regions in the depth direction, the apparatus acquires 512×512×4 tomograms. Note that at least a tomogram corresponding to one line is required for alignment, and there is no need to use all tomograms corresponding to 512×512 lines.
In step S407, it is determined whether there is any change in condition. The measurement conditions are changed when, for example, images in the depth direction are acquired in parts upon collective acquisition of them, or when there are both a place to be measured precisely and a place to be measured coarsely.
In step S408, the apparatus adjusts the image. If there is a collectively acquired tomogram as shown in
Image adjustment without using any collectively acquired image in the depth direction will be described with reference to
In step S409, the apparatus connects the normalized data in the respective measurement regions to visualize the data.
In step S410, the apparatus returns the coherence gate and the focus to the initial positions and terminates the processing.
In this case, the apparatus performs calculation upon placing the coherence gate at the boundary between measurement regions. If, however, i of SM(i, k) is a low-order component, an error due to the spectrum of a light source may mix in the data. In such a case, it is possible to measure each measurement region while shifting the coherence gate to the minus side of k.
Although this arrangement is that of an SD-OCT, the arrangement can also be applied to an SS-OCT which sweeps the wavelength of the light source. In this case, the wavelength sampling intervals are changed.
According to the present invention, it is possible to correct the contrast differences between even images obtained by spectrometers having different resolutions. This makes it easy to compare a plurality of images obtained by using spectrometers having different resolutions.
Aspects of the present invention can also be realized by a computer of a system or apparatus (or devices such as a CPU or MPU) that reads out and executes a program recorded on a memory device to perform the functions of the above-described embodiment(s), and by a method, the steps of which are performed by a computer of a system or apparatus by, for example, reading out and executing a program recorded on a memory device to perform the functions of the above-described embodiment(s). For this purpose, the program is provided to the computer for example via a network or from a recording medium of various types serving as the memory device (for example, computer-readable storage medium).
While the present invention has been described with reference to exemplary embodiments, it is to be understood that the invention is not limited to the disclosed exemplary embodiments. The scope of the following claims is to be accorded the broadest interpretation so as to encompass all such modifications and equivalent structures and functions.
This application claims the benefit of Japanese Patent Application No. 2010-129352 filed on Jun. 4, 2010, which is hereby incorporated by reference herein in its entirety.
Number | Date | Country | Kind |
---|---|---|---|
2010-129352 | Jun 2010 | JP | national |
Number | Name | Date | Kind |
---|---|---|---|
6377349 | Fercher | Apr 2002 | B1 |
7466423 | Podoleanu et al. | Dec 2008 | B2 |
7823782 | Yatagai et al. | Nov 2010 | B2 |
20010033364 | Cabib et al. | Oct 2001 | A1 |
20090079993 | Yatagai et al. | Mar 2009 | A1 |
20100226553 | Suehira | Sep 2010 | A1 |
20100226554 | Suehira | Sep 2010 | A1 |
20110058175 | Suehira | Mar 2011 | A1 |
20110096333 | Suehira et al. | Apr 2011 | A1 |
20110098560 | Suehira et al. | Apr 2011 | A1 |
Number | Date | Country |
---|---|---|
11-325849 | Nov 1999 | JP |
Entry |
---|
U.S. Appl. No. 13/204,953, filed Aug. 8, 2011, Shimoyama et al. |
U.S. Appl. No. 13/282,778, filed Oct. 27, 2011, Matsumoto et al. |
U.S. Appl. No. 13/286,481, filed Nov. 1, 2011, Suehira et al. |
U.S. Appl. No. 13/287,380, filed Nov. 2, 2011, Suehira et al. |
Number | Date | Country | |
---|---|---|---|
20110299035 A1 | Dec 2011 | US |