The invention concerns an imaging apparatus for in vivo full-field interference imaging of a scattering sample.
The invention relates to the field of in vivo imaging, and more particularly to in vivo imaging of ophthalmic tissues.
Optical coherence tomography (OCT) is an interferometric technique which has become a powerful imaging modality. OCT has applications in a broad spectrum of areas, and in particular in applications of biomedical fields such as ophthalmology, dermatology, cardiovascular field, gastroenterology, and so on. Optical coherence tomography includes, for instance, full-field OCT.
For instance, document WO 2020/065068 A1 discloses an imaging apparatus for full-field OCT designed for in vivo imaging of ophthalmic tissues.
The imaging apparatus includes a full-field OCT (FF-OCT) system. Furthermore, the imaging technique implemented with this apparatus is based on the use of light backscattered by a sample (i.e., in the case of a biological sample, by the microscopic cell and tissue structures of said sample) when it is illuminated by a probe beam emitted by a light source of the FF-OCT system, said light source having low spatial and temporal coherence length. The sample is, typically, the eye of a subject.
A microscope objective is also provided in front of the sample to focus light originating from the light source onto said sample, specifically onto the cornea. Alternatively, in the case of retinal imaging, the microscope objective is removed, and the cornea and lens of the eye ensure focusing of the probe beam on the retina.
The FF-OCT system further comprises an interferometer, including a reference arm and a sample arm separated by a beam splitter. In this case, the sample is arranged in the sample arm of the interferometer. The presence of the interferometer makes it possible to generate, by an interference phenomenon, an interference optical signal representative of the light originating selectively from a given slice of the sample volume, and to eliminate the light originating from the rest of the sample. This results from the low temporal coherence of the light source, which allows to isolate the light backscattered by a virtual slice depth-wise in the sample.
The FF-OCT system further includes a camera that captures, over time, 2D interference images representative of the interference optical signal.
Then, processing of interference images acquired based on the optical interference signal allows to retrieve at least one FF-OCT image, said FF-OCT image being an en-face image of the sample, that is, an image of the sample in a plane orthogonal to an optical axis of the sample arm.
Furthermore, a spectral-domain OCT (SD-OCT) system is provided to track the position of the eye. This SD-OCT system includes a secondary light source, a secondary beam of which is coupled through the illumination arm of the FF-OCT apparatus so as to propagate along the probe beam and impinge onto the sample. Analysis of the photons backscattered by the sample in response to illumination by the secondary beam gives access, in real time, to a location of the sample. Knowledge of the position of the eye allows to compensate for eye movement (especially in the propagation direction of the probe beam), in real time.
Moreover, by scanning the secondary beam, using a scanner, cross-sections of the sample are also available.
However, such imaging apparatus is not fully satisfactory.
Indeed, in such imaging apparatus, the FF-OCT system and the SD-OCT system are strongly interrelated, as several optical elements are shared, which renders assembly and alignment of such apparatus very challenging.
Furthermore, in such known architecture, the scanner of the SD-OCT system is located far (typically 30 cm) from the microscope objective. As a result, only light propagating under small angles with respect to the propagation direction of the probe beam will enter the entrance pupil of the microscope objective. Consequently, the cross-sectional image provided by the SD-OCT system only has a few hundreds of pixels in the lateral dimension (i.e., a direction orthogonal to the propagation direction of the probe beam in the sample arm of the interferometer). This results in the SD-OCT system providing SD-OCT images of unsatisfying quality. Moreover, such SD-OCT system has a small field of view (around 1 mm).
Moreover, in such configuration, the scanner is not optically conjugated to the back focal plane of the objective. As a result, the secondary beam is not incident perpendicular to the sample surface, which leads to low-signal in the SD-OCT images, since perpendicular direction produces the highest back-reflectivity. Consequently, the signal-to-noise ratio in the cross-sections is too low to evaluate, in real-time, the health of the ocular layers.
Furthermore, in such known apparatus, achieving conjugation of the scanner to the image back focal plane of the objective is very impractical. Indeed, even though the addition of a 4F set-up, on the SD-OCT system side, could allow such conjugation, this would result, due to the extended path of the SD-OCT beam in the FF-OCT system (about 20 cm), in an unacceptably long 4F set-up (around 80 cm) which complicates any commercial exploitation.
Furthermore, since the SD-OCT system is coupled through the illumination arm of the FF-OCT apparatus, the beam splitter arranged in the interferometer of said FF-OCT apparatus causes a 75% signal loss on the SD-OCT system side, compared to a situation where the secondary beam output by the SD-OCT system directly impinges onto the sample. This is prejudicial to the computation of high-quality cross-sectional images of the sample, as well as to precise location of the moving sample.
A purpose of the invention is to provide an imaging apparatus that overcomes these problems.
To this end, the present invention is an imaging apparatus of the aforementioned type, comprising:
Indeed, such architecture provides a full-field OCT imaging system and a sample locating system that share a minimal number of components (namely, the sample beam splitter), rendering these two systems almost optically independent. This enables easier manufacturing, assembling and maintenance.
Furthermore, in such design, due to the very short distance between the sample locating system and the sample, the optical conjugation between the sample locating system and focal plane in the sample is improved leading to close to perpendicular direction of scanning angles and improvement of the signal.
Furthermore, this design allows to build a conjugation device with a short footprint between the beam scanner and the sample for the best optical performance and signal. As a result, scanning of the sample can be performed in a very precise way, thereby providing better location of the sample. Moreover, a larger field of view is achieved (about 4-5 mm).
Moreover, since the sample locating system is not coupled through the illumination arm of the full-field OCT imaging system, the aforementioned reduction of the signal on the sample locating system side is not experienced.
According to other advantageous aspects of the invention, the imaging apparatus includes one or more of the following features, taken alone or in any technically possible combination:
The invention also relates to an imaging method for in vivo full-field interference imaging of a scattering sample, the method comprising:
According to other advantageous aspects of the invention, the imaging method includes one or more of the following features, taken alone or in any technically possible combination:
The invention will be better understood with the attached figures, in which:
An imaging apparatus 2 according to the invention is shown on
The imaging apparatus 2 is particularly suitable for in vivo imaging of a scattering sample 4, especially a moving sample, typically a subject's eye. More precisely, the imaging apparatus 2 is configured for full-field interference imaging of the sample 4.
In particular, the imaging apparatus 2 is configured to determine at least one en-face image (and, preferably, at least one cross-sectional image) of at least part of the sample 4, such as the cornea 6 of the in vivo eye 4. Preferably, the imaging apparatus 2 is also configured to provide at least one 3D image of at least part of the sample 4, for instance based on the determined en-face image and cross-sectional image.
The imaging apparatus 2 comprises a full-field OCT imaging system 8 (also referred to as “FF-OCT imaging system”), an optical sample locating system 10, and a processing unit 12.
The FF-OCT imaging system 8 is configured to provide at least one interference image representative of an en-face image of the sample 4, said en-face image being computed by the processing unit 12.
According to the invention, an en-face image of the sample 4 is an image of an in-depth section of the sample 4. Such in-depth section corresponds to a section of the sample 4 according to a plane that is orthogonal to an optical axis Z-Z of the FF-OCT imaging system 8.
Furthermore, the optical sample locating system 10 is configured to output a location signal representative of a location of the sample 4 along the optical axis Z-Z. The location signal is, for instance, an optical spectrum, or a temporal signal having an evolution that is representative of the movements of the sample.
Preferably, the processing unit 12 is also configured to compute at least one cross-sectional image of at least part of the sample 4 based on the location signal.
Preferably, the processing unit 12 is also configured to compute the at least one 3D image of at least part of the sample 4 based on a time sequence of en-face images and associated cross-sectional images.
Preferably, the imaging apparatus 2 further comprises a main movable platform 13 for mounting the FF-OCT imaging system 8 and the optical sample locating system 10 (for instance a motorized platform, or a manual platform that is hand-controlled by an operator), the main movable platform being movable along three independent (for instance perpendicular) directions so as to jointly translate the FF-OCT imaging system 8 and the optical sample locating system 10.
A description of the FF-OCT imaging system 8 will now be provided.
The FF-OCT imaging system 8 comprises a primary light source 14, a primary interferometer 16, a camera 18 and a sample beam splitter 19.
The primary light source 14 is configured to output a probe beam 20 (dashed line on the figure), and is arranged at an input of the primary interferometer 16.
Preferably, the spectrum of the probe beam 20 has a central peak at a wavelength belonging to the [0.4 μm; 1 μm] interval, even more preferably the [0.75 μm; 1 μm] interval for eye examination.
The primary light source 14 is a source having low spatial coherence (around half the wavelength, for instance lower than 1.5 μm when using a second lens system 36 having a 0.3 numerical aperture) and a low temporal coherence length (for instance less than 20 μm). For instance, the primary light source 14 is a thermal light source (such as a halogen lamp), a light-emitting diode, or a laser followed by mode-mixing multimode fiber.
Preferably, an optical element 22 is arranged at the output of the primary light source 14 to achieve Köhler illumination.
The primary light source 14 is arranged so that the probe beam 20 emitted by said primary light source 14 reaches a primary beam-splitter 24 (described below) of the primary interferometer 16, thereby having a first part thereof propagating in a sample arm 26 of the primary interferometer 16, and a second part thereof (distinct from the first part of the probe beam 20) propagating in a reference arm 28 of the primary interferometer 16.
The primary interferometer 16 includes the aforementioned primary beam splitter 24, sample arm 26 and reference arm 28. The sample arm 26 is intended to receive the sample 4.
The primary interferometer 16 is configured to generate a primary optical interference signal resulting, when the sample 4 is arranged in the sample arm 26, from interferences produced, at each point of an imaging field of the FF-OCT imaging system 8, between a sample wave obtained from the sample arm 26, and a reference wave received from the reference arm 28. The sample wave and the reference wave will be further described below.
Due to the low-temporal coherence of the primary light source 14 the interference phenomena only happens when the sample and reference arms 26 and 28 of the primary interferometer 16 match with micron precision (typically 1-10 microns). In the present embodiment, when the sample 4 is placed inside the sample arm 26 of the primary interferometer 16, the interference occurs only for the light reflected from the few-micron-thick sample layer. Light reflected from the other sample layers does not interfere with the reference light coming from the referent arm 28.
Preferably, as shown on
Each element of the primary interferometer 16 will now be described in more detail.
The primary beam splitter 24, which separates the sample arm 26 and the reference arm 28, is, for instance, a non-polarizing splitter cube.
As previously stated, the sample arm 26 is intended to receive the sample 4. The sample arm has a respective optical axis, which defines the optical axis Z-Z of the FF-OCT imaging system 8.
The sample arm 26 is arranged so that the sample 4 receives, from the primary beam splitter 24, the aforementioned first part of the probe beam 20.
Furthermore, the aforementioned sample wave is obtained by scattering of incident light waves of said first part of probe beam 20 by the sample 4. More precisely, at each point of the imaging field, the corresponding sample wave is obtained by scattering (e.g., backscattering) of incident light waves of said first part of the probe beam 20 by an elementary volume of a slice of the sample 4, said elementary volume being conjugate to said point of the imaging field, all elementary volumes having the same given depth, i.e., the same coordinate along the optical axis Z-Z.
The reference arm 28 has an optical axis X-X, and comprises a reference reflection surface 30.
The reference arm 28 is arranged so that reference reflection surface 30 receives, from the primary beam splitter 24, the aforementioned second part of the probe beam 20, distinct from the first part.
In this case, the aforementioned reference wave is obtained by reflection of incident light waves of said second part of probe beam 20 on the reference reflection surface 30. More precisely, at each point of the imaging field, the corresponding reference wave is obtained by reflection of incident light waves of said second part of the probe beam 20 on an elementary surface of the reference reflection surface 30 conjugate to said point of the imaging field.
Preferably, the reference reflection surface 30 is flat. For instance, the reference reflection surface 30 is a metallic mirror, a neutral density filter glass, or a glass plate.
Alternatively, the reference reflection surface 30 is curved, preferably with a curvature that matches the convex curvature of the cornea. This feature is advantageous, as it results in an en-face image that does not lie in a plane orthogonal to the Z-Z axis, but rather represents a whole cell layer of the cornea at once.
Advantageously, the reference reflection surface 30 is mounted on a piezoelectric stage 32. Such piezoelectric stage 32 may be used for phase modulation, as will later be described in relation to the processing unit 12. Preferably, the piezoelectric stage 32 is configured to operate in the [100 Hz; 10 KHz] range.
The camera 18 is arranged at an output of the primary interferometer 16.
The camera 18 is configured to acquire at least one interference image representative of the aforementioned primary optical interference signal resulting, when the sample 4 is arranged in the sample arm 26, from interferences in the imaging field between the reference waves and the sample waves.
Each interference image is a 2D image which, once processed by the processing unit 12, allows to retrieve an en-face image of the sample 4, that is, an image of the sample in a plane orthogonal to the optical axis Z-Z of the sample arm 28 and conjugate to the imaging field.
Preferably, the camera 18 is configured to acquire the interference images at a high rate, for example with a frequency comprised between 100 Hz and 1 kHz, or even higher, depending on the dynamics of the sample 4, and more specifically the dynamics of the movements within the sample 4.
For instance, the camera 18 is a charge-coupled device (CCD) camera or a complementarity metal oxide semiconductor (CMOS) camera.
The sample beam splitter 19 is arranged in the sample arm 26 and crosses the optical axis Z-Z of the sample arm 26. More precisely, the sample beam splitter 19 is arranged downstream the primary beam splitter 24 in a direction from the primary light source 14 to the primary beam splitter 24. In other words, as shown on
The sample beam splitter 19 is, preferably, configured to transmit light having a wavelength within a predetermined range around the central wavelength of the primary light source 14, and to reflect light having a wavelength within a predetermined range around a central wavelength of a secondary light source 46 (described below) of the optical sample locating system 10. For instance, the sample beam splitter 19 is a dichroic mirror.
Advantageously, the primary interferometer 16 also comprises a reference beam splitter 33 arranged in the reference arm 28, between the primary beam splitter 24 and the reference reflection surface 30. This is advantageous, as such reference beam splitter 33 introduces, in the reference arm 28, dispersion similar to that cause by the sample beam splitter 19 in the sample arm 26, thereby improving the ability to extract en-face images from the interference images acquired thanks to the FF-OCT imaging system 8.
The reference beam splitter 33 has the same optical properties as the sample beam splitter 19. For instance, the reference beam splitter 33 is identical to the sample beam splitter 19.
Preferably, the reference beam splitter 33 and the sample beam splitter 19 are oriented so that they define a same angle with a propagation direction of light waves of the probe beam 20 propagating through each of them. In other words, the reference beam splitter 33 and the sample beam splitter 19 are oriented so that an angle between the plane into which the sample beam splitter 19 extends and a propagation direction of the first part of the probe beam 20 is equal to an angle between the plane into which the reference beam splitter 33 extends and a propagation direction of the second part of the probe beam 20.
More precisely, for each point of an illumination field of the probe beam 20:
This feature is advantageous, since the FF-OCT imaging system 8 may be sensitive to small mismatches between the sample arm 26 and the reference arm 28 (contrary to conventional SD-OCT systems). Indeed, by having the angles between the identical beam splitters 19, 33 in the sample arm 26 and the reference arm 28 match, the dispersion of the arms 26, 28 of the primary interferometer 16 are balanced. As a result, with appropriate tuning, the optical path length difference between the sample arm 26 and the reference 28 becomes smaller than the coherence length of the primary light source 14 (for instance, 8 μm), thereby allowing the FF-OCT imaging system 8 to acquire reliable interference images.
This feature also allows to provide lateral (i.e., perpendicular to the propagation of light) optical matching in the sample arm 26 and the reference arm 28, thereby further increasing the ability of the FF-OCT imaging system 8 to acquire reliable interference images.
Preferably, the primary interferometer 16 further comprises a first lens system 34 and a second lens system 36, arranged respectively on each of the sample arm 26 and the reference arm 28.
More precisely, the first lens system 34 is arranged so that the focal plane of said first lens system 34 matches the section of the sample 4 that is imaged, when the sample 4 is positioned in the sample arm 26. Furthermore, the first lens system 34 is arranged so that a corresponding optical axis coincides with the optical axis Z-Z of the sample arm 26.
In practice, the depth of focus of the first lens system 34 is always more extended, along the optical axis Z-Z, than the slice width which actually causes interference, due to the low temporal coherence of the primary light source 14. In other words, the camera 18 captures light from a relatively thick region of the sample 4 but optical interference signals resulting in the en-face image are located in a thinner sub-region of the depth of focus, the thickness of which is determined by the primary light source 14. By acquiring several camera frames with different phases, one can extract only the thin interference fringe section.
Preferably, the first lens system 34 is arranged downstream the sample beam splitter 19 in a direction from the primary beam splitter 24 to the sample beam splitter 19.
Furthermore, the second lens system 36 is arranged so that the reference reflection surface 30 is located at the focal plane of the second lens system 36.
Moreover, the second lens system 36 is arranged so that a corresponding optical axis coincides with the optical axis X-X of the reference arm 28.
Advantageously, each lens system 34, 36 has a high numerical aperture (typically around 0.3) to ensure cellular resolution while providing a large field-of-view (typically around 1 mm). For instance, each of the first lens system 34 and the second lens system 36 is a microscope objective.
Advantageously, the reference reflection surface 30 and the second lens system 36 are mounted on a motorized platform 40, movable along the optical axis X-X of the reference arm 28.
This is advantageous, as displacement of the motorized platform 40 along the optical axis X-X allows to position the imaged slice based on the position of the sample. More precisely, when imaging retina of the eye 4, displacement of the motorized platform 40 along the optical axis X-X allows to select the depth of the imaged slice (i.e., the plane orthogonal to the optical axis Z-Z that is imaged by the FF-OCT imaging system 8). Furthermore, when imaging cornea, displacement of the motorized platform 40 allows to compensate for defocus caused, in the eye, by Snell's law, and which depends on the position of the eye relative to the FF-OCT imaging system 8, so that the imaged slice falls in the current focus plane.
Preferably, the primary interferometer 16 further comprises a primary spectral filter 42 arranged between an output of said primary interferometer 16 and the camera 18.
The primary spectral filter 42 is configured to transmit light having a wavelength within a predetermined range around the central wavelength of the primary light source 14, and to block light having a wavelength outside said predetermined range, such as light emitted by the optical sample locating system 10. This prevents the camera 18 from detecting photons originating from the locating system 10.
For instance, the primary spectral filter 42 is an optical bandpass filter.
Preferably, the primary interferometer 16 also includes a third lens system 44, such as an achromatic doublet, arranged between an output of said primary interferometer 16 and the camera 18.
The third lens system 44 has a focal length selected to conjugate a detecting surface of the camera 18 with the planes situated at the foci of the first and second optical systems 34, 36. For instance, the focal length of the third lens system 44 is a few hundreds of millimeters, typically 300 mm.
Such third lens system 44 (typically a tube lens) also allows to modify the magnification and field-of-view that the camera sees 18, which is advantageous. Indeed, the first and second lens systems 34, 36 are generally “infinity corrected”. For instance, the system of Nikon microscope objective lens 34 (having the Tube Lens Focal Length equal to 200 mm) together with the 200 mm tube lens 44 would produce 10× magnified view, while the system with the 300 mm tube lens 44 would give a 15× magnification.
Preferably, the primary interferometer 16 further includes a tunable (i.e., adaptive) lens 45 arranged in the sample arm 26.
The tunable lens 45 is operable by the processing unit 12, advantageously in combination with the motorized platform 40, to compensate for defocus on the FF-OCT imaging system side, due to the fact that the eye does not have the same refractive index as air. More precisely, the tunable lens 45 is operable to control the position of the focus point of the first part of the probe beam 20, which impinges on the sample 4.
Optionally, the tunable lens 45 may be configured to provide similar functionalities as those of the third lens system 44. In this case, the third lens system 44 may be removed.
A description of the locating system 10 will now be provided.
As stated previously, the optical sample locating system 10 (also referred to as “locating system”) is configured to output a location signal representative of a location of the sample 4 along the optical axis Z-Z.
The locating system 10 is arranged outside the sample arm 26, and has a respective optical axis Δ intersecting the sample beam splitter 19. Due to the presence of the sample beam splitter 19, a location beam 51 emitted by the locating system 10 along its optical axis Δ is reflected by the sample beam splitter 19 so as to propagate substantially parallel to the optical axis Z-Z of the FF-OCT imaging system 8. Therefore, after reflection on the sample beam splitter 19, the optical axis Δ and the optical axis Z-Z are co-located.
It is to be noted that the direction of propagation of the location beam 51 (i.e., the optical axis Δ) is truly parallel to the optical axis Z-Z only for a single orientation of a deflection element 74 (described below) of the locating system 10. For all the other positions of the deflection element 74 (all of them being used during acquisition of the location signals leading to a 2D cross-sectional image of the sample 4), there exists an angle between the optical axis Δ and the optical axis Z-Z.
Conversely, light from the sample 4 (with appropriate wavelength) is able to reach the locating system 10 after being reflected by the sample beam splitter 19.
Preferably, the locating system 10 includes an OCT imaging system (such as a spectral-domain OCT system (also referred to as “SD-OCT system”), a swept-source OCT system or a time-domain OCT system), as will be described below.
In this case, the locating system 10 includes the aforementioned secondary light source 46, a secondary interferometer 48 and a detector 50.
Preferably, the locating system 10 further includes a conjugation device 53, shown on
The secondary light source 46 is configured to output the aforementioned location beam 51 (dotted line on the figure), and is arranged at an input of the secondary interferometer 48.
Advantageously, the location beam has a spectrum that at least partially does not overlap with a spectrum of the light emitted by the primary light source 14.
For instance, the spectrum of the location beam 51 has a central peak at a wavelength belonging to the [870 nm; 1000 nm] interval while the primary light source 14 emits light in the [750 nm; 860 nm] interval.
Preferably, the secondary light source 46 is a spatially coherent light source. For instance, the secondary light source 46 is a super-luminescent diode or a swept laser source (when using a swept-source OCT system).
Preferably, in the case where the locating system 10 includes an SD-OCT system, the secondary light source 46 has low temporal coherence (i.e., a broad spectrum compared to its central wavelength). Furthermore, in the case of a swept laser source, at each moment, the location beam has quite large temporal coherence, but is swept over a broad spectrum.
The secondary light source 46 is coupled to an input of the secondary interferometer 48, for instance through a first optical fiber 52.
The secondary light source 46 is arranged so that the location beam 51 emitted by secondary light source 46 reaches a secondary beam-splitter 54 (described below) of the secondary interferometer 48, thereby having a first part thereof propagating in a secondary sample arm 56 of the secondary interferometer 48, and a second part thereof (distinct from the first part of the location beam 51) propagating in a secondary reference arm 58 of the secondary interferometer 48.
The secondary interferometer 48 includes the aforementioned secondary beam splitter 54, secondary sample arm 56 and secondary reference arm 58.
The secondary beam splitter 54 separates the secondary sample arm 56 and the secondary reference arm 58. Moreover, the secondary sample arm 56 is intended to receive the sample 4.
The secondary interferometer 48 is configured to produce a secondary optical interference signal resulting, when the sample 4 is arranged in the secondary sample arm 56, from interferences produced between a secondary sample wave obtained from the secondary sample arm 56, and a secondary reference wave received from the secondary reference arm 58. The secondary sample wave and the secondary reference wave will be further described below.
Due to the architecture of the imaging apparatus 2, arranging the sample 4 in the sample arm 26 of the primary interferometer 16 is equivalent to arranging the sample 4 in the secondary sample arm 56 of the secondary interferometer 48.
Preferably, the secondary reference arm 58 includes an optical fiber 60, an optional lens 62, an optional dispersion compensation plate 64 and a secondary reflector 66 (such as a metalized mirror) arranged in series.
In this case, the aforementioned secondary reference wave is obtained by reflection of incident light waves of the second part of the location beam on the secondary reflector 66.
Advantageously, the length of the optical fiber 60 is chosen so that an optical path between the secondary beam splitter 54 and the secondary reflector 66 is equal to an optical path between the secondary beam splitter 54 and the cell layer in the sharp focus of first lens system 34, that is the cell layer that is currently imaged. Preferably, the length of the optical fiber 60 is chosen so that the optical path between the secondary beam splitter 54 and the secondary reflector 66 is equal to an optical path between the secondary beam splitter 54 and a predetermined position where the sample 4 is supposed to be located during operation of the imaging apparatus 2.
Preferably, the secondary sample arm 56 includes an optical fiber 68, a beam scanner 70 and a secondary spectral filter 72 arranged in series.
In this case, the aforementioned secondary sample wave is obtained by scattering of incident light waves of the first part of the location beam by the sample 4. More precisely, the sample wave is obtained by scattering of incident light waves of the first part of the location beam by elementary volumes of the sample 4 located along the direction of propagation of the first part of the location beam after deflection by the beam scanner 70.
The beam scanner 70 is arranged between the sample beam splitter 19 and the secondary beam splitter 54. The beam scanner 70 is configured to scan the location beam in a plane orthogonal to the optical axis X-X. More precisely, the beam scanner 70 includes a deflection element 74 (typically, a mirror) for deflecting the location beam 51. Consequently, due to the reflection of the location beam on the sample beam splitter 19, operation of the beam scanner 70 (i.e., movement of the deflection element 74) results, at the sample 4, in a scanning of the location beam in a plane orthogonal to the optical axis Z-Z, more precisely in a scanning of the location beam along multiple lines in a plane orthogonal to the optical axis Z-Z, each line corresponding to the direction of propagation of the first part of the location beam 51.
For instance, the beam scanner 70 is configured to scan a line every 10 us approximately. In other words, the beam scanner 70 is configured to scan the eye 4 during a time interval that is shorter than a fluctuation time constant of the eye position, that is to say shorter than the typical time during which the eye position can fluctuate with a large range (greater than 10 micrometers), typically around 10 ms.
The secondary spectral filter 72 is configured to transmit light having a wavelength within a predetermined range around the central wavelength of the secondary light source 46, and to block light having a wavelength outside said predetermined range, such as light emitted by the FF-OCT imaging system 8. For instance, the secondary spectral filter 72 is an optical bandpass filter. This prevents the detector 50 from detecting photons originating from the FF-OCT imaging system 8.
The detector 50 is coupled to an output of the secondary interferometer 48.
The detector 50 is configured to detect the secondary optical interference signal and to output the corresponding location signal over time.
For instance, in the case where the locating system 10 is a time-domain OCT system of a swept-source OCT system, the detector 50 is a photodetector such as a photodiode. Alternatively, in the case where the locating system 10 is a spectral-domain OCT system, the detector 50 is a spectrometer.
The conjugation device 53 includes two lenses arranged in a 4F configuration, as shown on
More precisely, the conjugation device 53 is arranged so that a back focal plane of the first lens system 34 coincides with a first focal plane F1 of the conjugation device 53, and the deflection element 74 of the beam scanner 70 is arranged in a second focal plane F2 of the conjugation device 53.
Thanks to the presence of the conjugation device 56, the location beam 51 is scanned, by the beam scanner 70, in a plane orthogonal to the optical axis X-X and located at about the same depth as the imaging field determined by the first lens system 34. This is shown by the crosshatched beam and the dotted beam on
As mentioned before, due to the low-temporal coherence of the primary light source 14, the interference phenomena only happens when the sample and reference arms 26 and 28 of the primary interferometer 16 match with micron precision (typically 1-10 microns).
From the first impression it may seem that due to the low-temporal coherence of the secondary light source 46 of the SD-OCT system, the interference phenomena will happen only, when the secondary sample and secondary reference arms 56 and 58 of the secondary interferometer 48 match with micron precision (like in the FF-OCT system). However, the inventors discovered that it is not the case in the SD-OCT version of the optical sample locating system 10 where the detector is a spectrometer that is composed of a diffraction grating and a detector. In this case, grating separates the broadband low temporally coherent spectrum of the reflected secondary light source 46 (coming respectively from the secondary sample and secondary reference arms 56 and 58) into the thin highly temporally coherent bands. Each of these bands are then detected by a separate detector. Therefore, for each detector the interference appears to happen even if the secondary sample and secondary reference arms 56 and 58 of interferometer are largely mismatched (by 1-10 millimeters).
In order to still reconstruct a tomographic image, SD-OCT relies on the fact that the interference frequency is different for different mismatch distances. The interference frequency is small, if the secondary sample and secondary reference arms 56 and 58 are closely matching, and gets larger the bigger is the mismatch. By taking a Fourier transform on the measured spectrum all the frequencies are extracted and then can be converted to distances. As a result, SD-OCT can reconstruct a 1D in-depth profile (A-line) of the scatterers inside the sample.
2D cross-sectional image can be reconstructed by placing a beam scanner in the secondary sample arm 56, like the deflection element 74. Notice that there is no scanning in the secondary reference arm 58 and different secondary sample beams (coming from the secondary sample arm 56) interfere with a single secondary reference beam (coming from the secondary reference arm 58).
Thus, SD-OCT is a point-by-point imaging system that is not affected by the wavefront mismatches in the secondary sample and secondary reference arms 56 and 58, contrary to FF OCT system 8, which captures simultaneously all the pixels of en-face tomographic image.
Furthermore, the inventors discovered that SD-OCT system and FF-OCT system 8 are not similarly affected by the insertion of a dichroic mirror in the sample and secondary sample arms 26 and 56, like the sample beam splitter 19. For instance, in the case of a dichroic mirror which would be a 1 mm glass plate, the glass plate results in a 1*1.5−1*1=0.5 mm (thickness*refractive index) optical path difference between the sample (respectively secondary sample) and reference (respectively secondary reference) arms. It appears that the SD-OCT system still exhibit interference and the final image is not affected except by a depth shift and a bit of dispersion leading to reduced axial resolution. Thus, wavefront of the SD-OCT system does not require correction because it's a point imaging system.
On another hand in the FF-OCT system 8, the inventors discovered that the 0.5 mm mismatch appears to destroy all the interference signals because the optical arms require to be aligned with few microns precision. No image could be retried. Additionally, as FF-OCT system 8 is a full-field and not point-by-point imaging system the signal, it suffers from the wavefront mismatch.
Furthermore, the inventors discovered that a glass plate placed in the FF-OCT system 8 leads to curvature of the wavefront. One way to compensate optical mismatch, dispersion and wavefront mismatch in FF-OCT is to insert exactly the same dichroic mirror (as described later in the second embodiment) in the reference arm 28 as in the sample arm 26 and rotate it to exactly the same angle.
Achieving exactly the same angle with micron optical path precision is not easy, therefore it requires an additional aligning step. First, the FF-OCT system 8 interferometer is aligned in a usual way, but with two dichroic mirrors 19 and 33 that are placed in the sample and reference arms 26 and 28 perpendicular to the optical paths of the light. After obtaining an interference image on the screen, one dichroic mirror 19 or respectively 33 is tilted. The resulting mismatch leads to interference signal disappearing. Then the second dichroic mirror 33 or respectively 19 is tilted to recover the signal. Then, this process step-by-step is continued step by step: small tilt of one dichroic mirror followed by the small tilt of another dichroic mirror to recover the signal. The final alignment is achieved, when both dichroic mirrors 19 and 33 are tilted at 45 degrees.
It is to be noted that, in regard to the FF-OCT system of document WO 2020/065068 A1, the problem of signal loss is not applicable, because the dichroic mirror was placed in the illumination arm (i.e. the probe beam 20), thus sample and reference arms stayed aligned. Due to the sensitivity of the FF OCT system compare to the SD-OCT system, another configuration that solves the technical problems cited in the background of the invention seemed to be impossible before the present invention.
According to the present invention, the expression “processing unit” should not be construed to be restricted to hardware capable of executing software, and refers in a general way to a processing device, which can for example include a microprocessor, an integrated circuit, or a programmable logic device (PLD). The processing unit may also encompass one or more Graphics Processing Units (GPU) or Tensor Processing Units (TSU), whether exploited for computer graphics and image processing or other functions. Additionally, the instructions and/or data enabling to perform associated and/or resulting functionalities may be stored on any processor-readable medium such as, e.g., an integrated circuit, a hard disk, a CD (Compact Disc), an optical disc such as a DVD (Digital Versatile Disc), a RAM (Random-Access Memory) or a ROM (Read-Only Memory). Instructions may be notably stored in hardware, software, firmware or in any combination thereof.
The processing unit 12 is configured to compute each en-face image based on each interference image output by the FF-OCT imaging system 8.
To do so, the processing unit is configured to control the piezoelectric stage 32 so as to displace the reference reflection surface 30 along the optical axis X-X, thereby modulating the reference wave. The amplitude of such modulation is set according to a predetermined phase demodulation scheme, for instance 2-phase demodulation scheme or 4-phase demodulation scheme.
For instance, in a 2-phase scheme, the processing unit 12 is configured to determine the en-face image from the two interference images acquired by the camera 18: the first interference image is captured with the reference reflection surface 30 is in a first position, and the second interference image is captured with the reference reflection surface 30 shifted by an optical phase of π, i.e., half the central wavelength of the probe beam 20.
As previously said, the interference occurs only for the light reflected from the few-micron-thick sample layer. Light reflected from the other sample layers does not interfere with the reference light coming from the referent arm 28. This modulation of the reference reflection surface 30 in time enables to acquire several interference patterns (optical phases). In the simplest demodulation scheme the two opposite (pi phase shifted) interference patterns are subtracted, which rejects the non-interfering light from the other sample layers and doubles the signal from the layer of interest (tomography).
As another example, in 4-phase scheme, the processing unit 12 is configured to determine the en-face image from the four interference images acquired by the camera 18, the images being captured sequentially. If one assumes that the phase of the initial frame is 0, then the processing unit 12 is configured to control the piezoelectric stage 32 and the camera 18 so that the phase shifts for the next frames, relatively to the first frame, are π/2, π and 3π/2.
The processing unit 12 is further configured to compute a reflectivity of the sample structures within the coherence volume (i.e., the plane of interest) based on a result of subtracting said phase shifted images by pairs. The resulting en-face image, an example of which is shown on
Moreover,
Imaging trabecular meshwork is key for diagnosing glaucoma, a condition that is detrimental to vision and that is extremely difficult to diagnose at an early stage. The ability to image the trabecular meshwork at high resolution (which is not provided by known imaging systems) allows to detect glaucoma at an early stage and avoid the risks of losing vision.
Alternatively, or in addition, in the case of a sample moving along the Z-Z axis (typically, the eye exhibiting the natural physiological eye movements), the processing unit 12 is configured to compute each en-face image based on each interference image output by the FF-OCT imaging system 8, without modulation of the position of the reference reflection surface 30. In this case, small movements of the sample along the Z-Z axis result in interference images that are phase-shifted with respect to each other, thereby allowing computation of the en-face images.
Preferably, the processing unit 12 is also configured to control the position of the motorized platform 40, based on the location signal, to compensate for sample movement, so that a same sample slice is imaged over time.
More precisely, for a given position of the deflection element 74, the detector 50 acquires a location signal representative of a line corresponding to the direction of propagation of the first part of the location beam 51. The processing unit 12 is classically configured to reconstruct the portion of the sample that is intersected by said line.
Preferably, the processing unit 12 further controls the beam scanner 70 so that the location signal for various positions of the deflection element 74 (that is, various neighboring lines) is acquired. The processing unit 12 is further configured to compute, based on said plurality of location signals, at least one cross-sectional image of at least part of the sample 4. An example of such cross-sectional image of the cornea is shown on
Furthermore, the processing unit 12 is configured to segment the cross-sectional image so as to identify a position of at least one given structure (such as a cellular structure) of the sample 4. As a result, by comparing each segmented cross-sectional images to a predetermined target central position, the processing unit 12 is able to determine, over time, the position, along the optical axis Z-Z, the position of each segmented structure.
Preferably, the processing unit 12 is further configured to perform image processing on each computed en-face image to correct noise arising from the operation of the FF-OCT imaging system 8.
Indeed, as shown on
The Gaussian noise originates from the detection, that is from the camera 18. This type of noise, which is not specific to FF-OCT, is present in any images taken by a camera. Gaussian noise is uncorrelated, which means that noise in one pixel is independent from noise in any other pixel.
Interference fringe noise, which can be seen, on
Interference fringes are correlated: for example, if one pixel is very bright, there is high chance that the nearby pixels are bright as well.
In order to perform noise correction, the processing unit 12 is configured to implement an artificial intelligence model, such as Noise2Noise deep learning model.
More precisely, the model has been previously trained using a set of training images, for instance in vivo human corneal FF-OCT images (such as stroma and endothelium).
More precisely, using the full-field OCT imaging system 10, at least one training set of en-face images of a training sample is acquired. Each training set includes a first en-face image and a second en-face image of a same section of the training sample 4 at distinct acquisition times. Due to the high speed of the camera 18 (around 500 frames/second), said en-face images can be regarded as if the eye movements were frozen (i.e., in the same location) for a few frames. In this case, for each training set, the first en-face image and the second en-face image exhibit different noises: this is due, on the one hand, to Gaussian noise (which is different from a camera frame to the other), and, on the other hand, to interference noise (which is caused by imperfect demodulation scheme and/or sub-micron movements of the eye from frame to frame leading to phase fluctuations from an acquired interference image to the other).
Then, to train the Noise2Noise model, for each training set, the first en-face image is provided as an input, and the second en-face image is provided as a target output.
This is advantageous. Indeed, during such training, the Noise2Noise model tries to learn how to create one specific noisy pattern from another noisy pattern. Since noise is random, the model cannot efficiently learn this transformation, and it ends up learning to transform the noisy pattern into the average of all noisy patterns, which is a clean image. Thank to such training, both spatially uncorrelated noise (Gaussian) and correlated noise (interference fringes) are removed, as shown on
Operation of the imaging apparatus 2 will now be described.
During a configuration step, alignment of the FF-OCT imaging system 8 and the optical sample locating system 10 is performed.
Moreover, during the configuration step, a modulation/demodulation scheme is preferably chosen by an operator and stored into the processing unit 12.
Preferably, in the case where the processing unit 12 is intended to perform noise correction, configuration of the Noise2Noise neural network is performed during the configuration step. More precisely, the neural network is trained, as described previously.
Then, the sample 4 is arranged in the sample arm 26 of the FF-OCT imaging system 8 and the secondary sample arm 56 of the optical sample locating system 10.
Then, the processing unit 12 controls the piezoelectric stage 32 to displace the reference reflection surface 30 along the optical axis X-X, thereby modulating the reference wave. The processing unit 12 also controls the camera 18 so that it acquires interference images in accordance with the predetermined modulation scheme, that is when the reference reflection surface 30 is in positions leading to a predetermined phase shift of the reference wave.
The processing unit 12 also controls the beam scanner 70 so that the location beam 51 scans the sample 4. The processing unit 12 simultaneously controls the detector 50 so that the location signal representative of each line being currently scanned by the location beam 51 is output.
Then, the processing unit 12 computes each en-face image based on each interference image output by the FF-OCT imaging system 8.
The processing unit 12 also computes, based on said plurality of location signals, the cross-sectional image of the sample 4. Furthermore, the processing unit 12 segments the computed cross-sectional image so as to identify a position of at least one given structure of the sample 4.
The processing unit 12 further compares the position of the segmented cross-sectional image to a predetermined target central position to determine whether the position of each segmented structure changes over time.
If the position changes, the processing unit 12 controls the position of the motorized platform 40 to compensate for sample movement, so as to compensate for defocus.
Preferably, the processing unit 12 further implements the trained Noise2Noise neural network to perform noise correction on each computed en-face image. More precisely, each computed en-face image is provided at the input of the neural network, which outputs a denoised image.
Alternatively, as previously mentioned, modulation of the reference reflection surface 30 is not performed. In this case, the processing unit 12 relies on phase shifts introduced, from one interference image to the other, by natural movement of the sample along the optical axis Z-Z to compute the en-face images.
A second embodiment of the imaging apparatus according to the invention is shown on
The imaging apparatus 2 of
To do so, the imaging apparatus 2 of
Preferably, the imaging apparatus further includes an adaptive lens 76 (for example a liquid lens) and/or a rotating glass plate 78 inserted in the sample arm 26 to compensate for aberrations and dispersion mismatch introduced by the eye 4 with respect to the dispersion of the second lens system 36 in the reference arm 28.
If the optional conjugation device 53 is provided, then the conjugation device 53 is arranged so that the first focal plane F1 of the conjugation device 53 coincides with the cornea-lens of the eye 4.
Furthermore, in such imaging apparatus, the processing unit 12 is advantageously configured so that, in operation, if the processing unit 12 determines that the sample 4 has moved (for instance by determining that a segmented structure of a cross-section of the sample has moved with respect to a predetermined target central position), the processing unit 12 controls the motorized platform 40 so that the plane of the retina 7 from which the interference image (of the FF-OCT imaging system 8) arises remains at the same position within the sample 4.
A third embodiment of the imaging apparatus according to the invention is shown on
The imaging apparatus 2 of
In this case, the reference arm 28 does not include a reference beam splitter, since the probe beam 20, in the sample arm 26, does not undergo propagation through the sample beam splitter 19 that would cause dispersion requiring compensation.
Such configuration can also be applied using the FF-OCT imaging system 8 of the imaging apparatus 2 of
A fourth embodiment of the imaging apparatus according to the invention is shown on
The imaging apparatus 2 of
For instance, as shown on
Alternatively, as shown on
A fifth embodiment of the imaging apparatus according to the invention is shown on
In this embodiment, the locating system 10 includes an active position sensing device configured to emit an interrogation wave 90 along the axis Δ (or “sensing axis”). The interrogation wave is reflected by the sample beam splitter 19 so that it impinges onto the sample 4.
The active position sensing is further configured to detect a response wave 92 backscattered by the sample 4 and reflected by the sample beam splitter 19 to propagate along the sensing axis Δ towards the locating system 10. Consequently, based on the detected response wave 92, the locating system 10 outputs the location signal representative of the position of the sample 4 along the Z-Z axis.
The interrogation wave may be light (e.g., laser light) or an acoustic wave.
For instance, the locating system 10 is configured to determine the location signal based on time of flight.
Number | Date | Country | Kind |
---|---|---|---|
22305066.7 | Jan 2022 | EP | regional |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/EP2023/051400 | 1/20/2023 | WO |