WO2008154349A2 - Method and apparatus for localized polarization sensitive imaging - Google Patents
Method and apparatus for localized polarization sensitive imaging Download PDFInfo
- Publication number
- WO2008154349A2 WO2008154349A2 PCT/US2008/066055 US2008066055W WO2008154349A2 WO 2008154349 A2 WO2008154349 A2 WO 2008154349A2 US 2008066055 W US2008066055 W US 2008066055W WO 2008154349 A2 WO2008154349 A2 WO 2008154349A2
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- light beam
- sample
- polarized
- horizontally
- light
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Ceased
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B9/00—Measuring instruments characterised by the use of optical techniques
- G01B9/02—Interferometers
- G01B9/02041—Interferometers characterised by particular imaging or detection techniques
- G01B9/02044—Imaging in the frequency domain, e.g. by using a spectrometer
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B9/00—Measuring instruments characterised by the use of optical techniques
- G01B9/02—Interferometers
- G01B9/02015—Interferometers characterised by the beam path configuration
- G01B9/02027—Two or more interferometric channels or interferometers
- G01B9/02028—Two or more reference or object arms in one interferometer
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B9/00—Measuring instruments characterised by the use of optical techniques
- G01B9/02—Interferometers
- G01B9/02055—Reduction or prevention of errors; Testing; Calibration
- G01B9/02075—Reduction or prevention of errors; Testing; Calibration of particular errors
- G01B9/02078—Caused by ambiguity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B9/00—Measuring instruments characterised by the use of optical techniques
- G01B9/02—Interferometers
- G01B9/0209—Low-coherence interferometers
- G01B9/02091—Tomographic interferometers, e.g. based on optical coherence
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/47—Scattering, i.e. diffuse reflection
- G01N21/4795—Scattering, i.e. diffuse reflection spatially resolved investigating of object in scattering medium
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B2290/00—Aspects of interferometers not specifically covered by any group under G01B9/02
- G01B2290/70—Using polarization in the interferometer
Definitions
- Embodiments relate generally to the field of biomedical imaging, specifically to methods, apparatuses, and systems associated with optical coherence tomographic imaging.
- Optical coherence tomography is a non-invasive, three- dimensional (3-D) imaging technique, capable of optical ranging within a highly scattering sample, such as biological tissue.
- PS polarization sensitive
- Spectral domain (SD) OCT may provide further advantages because of the significant sensitivity and imaging speed advantages over time domain OCT.
- Current techniques may not be capable of producing images of sufficient depth or clarity.
- conventional methods may require burdensome calibration steps for producing the desired images.
- FIGURE 1 illustrates an SD-PSOCT system in accordance with various embodiments
- FIGURE 2 illustrates an exemplary reference mirror in accordance with various embodiments
- FIGURES 3(A) and 3(B) are images of a chicken breast acquired using an exemplary embodiment of an SD-PSOCT system in accordance with various embodiments;
- FIGURE 4 is a plot of SD-PSOCT sensitivity against imaging depth in accordance with various embodiments
- FIGURES 5(A), 5(B), and 5(C) are images of a chicken breast acquired using an exemplary embodiment of an SD-PSOCT system in accordance with various embodiments;
- FIGURE 6 illustrates another SD-PSOCT system in accordance with various embodiments
- FIGURE 7 illustrates another exemplary reference mirror in accordance with various embodiments
- FIGURES 8(A) and 8(B) are images of a chicken breast acquired using an exemplary embodiment of an SD-PSOCT system in accordance with various embodiments.
- FIGURES 9(A), 9(B), and 9(C) are images of a human nail fold acquired using an exemplary embodiment of an SD-PSOCT system in accordance with various embodiments.
- the description may use perspective-based descriptions such as up/down, back/front, and top/bottom. Such descriptions are merely used to facilitate the discussion and are not intended to restrict the application of the disclosed embodiments.
- a phrase in the form of "A/B" or in the form "A and/or B” means
- a phrase in the form “at least one of A, B and C” means “(A), (B), (C), (A and B), (A and C), (B and C) or (A, B and C).
- a phrase in the form “at least one of A, B, and C” means "(A), (B), (C), (A and B), (A and C), (B and C), or (A, B and C).
- a phrase in the form “(A) B” means "(B) or (A B),” that is, A is optional.
- a computing system may be endowed with one or more components of the disclosed apparatuses and/or systems and may be employed to perform one or more methods as disclosed herein.
- an apparatus comprising a first light path configured to receive a first reference light beam and form, from the first reference light beam, a horizontally-polarized light beam; a second light path configured to receive a second reference light beam and form, from the second reference light beam, a vertically-polarized light beam; and a reference mirror configured to modulate, at a constant modulation frequency, an interference signal formed by the horizontally- and vertically- polarized light beams and the sample beam.
- a spectral domain polarization sensitive optical coherence tomography (SD-PSOCT) system having a single-camera configuration may comprise a reference arm assembly (hereinafter “reference assembly”) for polarization sensitive imaging.
- reference assembly for polarization sensitive imaging.
- FIGURE 1 An exemplary SD-PSOCT system including a reference assembly 102 is illustrated in FIGURE 1.
- the exemplary reference assembly 102 may be configured to operatively couple to an interferometer 104 and a spectrometer 106 including a single charge- coupled device camera (CCD) 108 for spectral interferogram detection to perform polarization sensitive imaging.
- the reference assembly 102 may be configured to separate by a fixed distance OCT images from orthogonal polarization channels at the output plane of Fourier space by introducing a constant frequency modulation in the spatial varying spectral interferograms when a probe beam, provided by a probe 110, is swept over a sample 112.
- the exemplary reference assembly Because the orthogonally polarized beam components travel in the identical beam path in the spectrometer 106, the distortion between spectra may be avoided. Further, because the beam components travel together, calibration of the beams to achieve pixel-by- pixel correspondence may be unnecessary. [0023] Still referring to FIGURE 1 , the exemplary reference assembly
- a light source 114 such as, for example, a broadband superluminescent diode (SLD).
- SLD broadband superluminescent diode
- an exemplary light source may be an SLD with a central wavelength of about 850 nm and full-width half-maximum (FWHM) bandwidth of about 30 nm, yielding a measured axial resolution of about 12 pm in air.
- the light source 114 may be a swept laser source, rather than the SLD.
- the spectrometer 106 may be replaced by an ultra-fast photodiode detector.
- an SD-PSOCT system including a swept laser source as the light source 114 and photodiode for detection may perform substantially similarly to embodiments including an SLD as the light source 114 and the spectrometer 106 for detection.
- Light from the light source 114 may be collimated by one or more lenses 116.
- the diameter of the collimated light may be about 3 mm.
- Light from the light source 114 may be vertically polarized by a linear polarizer 118 and introduced to a polarization sensitive low coherence interferometer 104 including a beam splitter 120.
- the vertically-polarized light leaving the interferometer 104 (hereinafter “reference light”) may then be delivered to the reference assembly 102 and the sample 112.
- the reference light delivered to the reference assembly 102 may be split by a non-polarization beam splitter 122 into two substantially equal beams.
- One of the beams (“first beam”) may pass through a quarter wave plate 124.
- the quarter wave plate 124 may have a fast axis oriented at 45°.
- the first beam may then be reflected by a staggered reference mirror 126.
- the first beam may then pass again through plate 124 to become a horizontally polarized beam.
- the other beam (“second beam”) split by the non-polarization beam splitter 122 may be bent by another mirror 128 and reflected back by the staggered reference mirror 126 to join with the horizontally polarized beam.
- the joined beam may then be routed back into the interferometer 104.
- an SD-PSOCT system described herein may include one or more optical fiber couplers in addition to or instead of a beam splitter. So, in various embodiments, one or both of the beam splitters 120 and 122 may instead be an optical fiber coupler, depending on the particular desired application.
- the reference assembly 102 may include a neutral density filter 130 placed in one or both of the first and second beam paths.
- the neutral density filter 130 may ensure that the light for two orthogonally polarized beam components (horizontal and vertical) impinged onto the CCD 108 has substantially equal strength.
- the reference assembly 102 may provide a fixed optical pathlength difference (OPD) between these two orthogonally polarized lights.
- OPD optical pathlength difference
- the fixed OPD may be adjusted according to the needs for different applications. For example, according to various embodiments, the OPD between the vertically and horizontally polarized lights may be adjusted to about 2.0 mm.
- the OCT images formed by the two orthogonally polarized beam components at the output plane of Fourier space may be separated after processing the mixed spectral interferograms captured by the CCD 108 (as described more fully herein).
- sample light may also be delivered to the sample 112 (hereinafter "sample light").
- sample light may pass through a quarter wave plate 132.
- the quarter wave plate 132 may have a fast axis oriented at 45° to transform the sample light into circularly polarized light.
- the sample light may then be delivered to the sample 112 by an X galvanometer scanner 134 and an achromatic lens 136.
- the achromatic lens 136 may include any focal length suitable for the application. In some embodiments, for example, the focal length may be about 100 mm.
- the sample light may be introduced to the sample by the probe
- the theoretical lateral imaging resolution, depth of focus, power, and other parameters of the probe beam may be variously adjusted depending on the application.
- the probe beam may have a theoretical lateral imaging resolution of about 28 ⁇ m, a depth of focus of about 4.2 mm, and a power of about 0.6 mW.
- the reference light returning from the reference arm 102 and the sample light returning from sample 110 may be recombined through a lens 152 and coupled into a single mode fiber 138 for introduction to the spectrometer 106.
- a polarization controller 140 may be used in the single mode fiber 138 for fine tuning the polarization states of light in the spectrometer system.
- the spectrometer 106 may be variously configured to be suitable for the intended purpose(s). According to the illustrated embodiment, the spectrometer 106 comprises a collimator 142, a polarization-insensitive transmission grating 144, and an achromatic focusing lens 146. According to an exemplary embodiment, the spectrometer 106 may comprise a 30 mm focal length collimator 142, a 1200 lines/mm polarization-insensitive transmission grating 144, and an achromatic focusing lens 146 with 190 mm focal length.
- the CCD 108 may comprise a line scan charge-coupled device camera for capturing both orthogonally polarized beam components (i.e., horizontal and vertical components) in parallel.
- the spectrometer 106 may include any CCD suitable for the purpose, in an embodiment, the CCD 108 may comprise one having 1024 pixels, each about 14x14 microns in size and 12 bits in digital depth, and capable of a 53 kHz line rate.
- the propagating beam paths for vertically and horizontally polarized light components are substantially identical in the spectrometer 106, and their spectra captured by the identical CCD pixel arrays in parallel. This may be particularly advantageous in that careful calibration of the spectral components on the CCD array may not be necessary.
- the spectrometer 106 may be designed to have an imaging depth being the full depth in the Fourier space used for imaging (e.g., approximately 8.2 mm in air), rather than half of the space used in conventional SD-PSOCT systems.
- the CCD 108 may be variously adjusted for imaging.
- the CCD integration time may be set at about 100 microseconds for imaging.
- the spectral data from the CCD 108 (e.g., 1024 pixels, A scan) may then be outputted to a computer 148 for access by a user.
- the image data is downloaded to the computer 148 using CAMERALINKTM and a high-speed frame grabber board (for example, PC1 1428, commercially available from National Instruments, USA).
- the CCD line scan rate may be 10 kHz
- the imaging speed may be limited by the output power of the light source used (about 5 mW)
- the light energy in the reference arm 102 may be adjusted such that the maximum intensity on the linear CCD array detector achieved approximately 40% of the camera's saturation level for each polarization beam component, and the remaining 20% of the dynamic range of the camera may be available to fill up the polarization interference fringes.
- a system comprising a light source; a reference assembly including: a first light path configured to receive a first reference light beam from the light source and form, from the first reference light beam, a horizontally-polarized light beam; a second light path configured to receive a second reference light beam from the light source and form, from the second reference light beam, a vertically-polarized light beam; and a signal modulating device configured to modulate an interference signal formed by the horizontally- and vertically-polarized light beams and a sample light beam to give a constant modulation frequency; a sampling assembly configured to provide light from the light source to a sample to obtain the sample light beam; and a spectrometer configured to detect a spectral interferogram formed by the horizontally-polarized light beam, the vertically-polarized light beam, and the sample light beam.
- the modulation of the interferogram may also be provided by scanner 134 and reference mirror 150 may be maintained in one position.
- scanner 134 may be configured such that the input signal is scanned with an offset reference to the pivot point.
- FIGURE 1 illustrated is a depiction of the relationship between the focused sample beam and the effective positions for the reference mirror 126. Assuming that the coordinate origin along the depth direction, z, is at the focal spot, the optical fields in the reference arm 102 for the vertically and horizontally polarized lights may be expressed as:
- Equation 1 H and V refer to the signals from the horizontal and vertical polarization channels, R is the reflectivity in the reference arm 102, k is the wavenumber, and d is the distance of the reference mirror 126 position relative to the focal spot.
- the optical fields that are decomposed into the vertically and horizontally polarized light may be expressed as:
- Equation 2 S H,V (Z) represents the strengths of backscattering light at optical depth z in the sample in the respective polarization direction.
- Equation 3 the problem of spectral interferogram detected by the CCD 108 at wavenumber k may be bundled and simplified by ignoring the terms that do not contribute to the useful signals, without loss of psychology, as:
- Equation 4 A H , ⁇ ( z ) represents the magnitude of envelopes of interferograms formed by the respective polarization components and ⁇ H,Y( Z ) ⁇ s tne random phase term that is due to the optical properties of the sample 112.
- equation 4 the interferograms formed by the vertical and horizontal polarization light beams may be added together at the CCD 108.
- the sample target may be placed near the zero phase delay line where the sensitivity may be the highest. This may be particularly desirable if the CCD 108 is insufficiently capable of digitizing the spectral components due to decreased sensitivity along the imaging depth.
- OCT images formed by the two orthogonal polarization states may be placed side by side approximately centered at the zero delay line. For some embodiments, the images are placed about 2 mm apart.
- spectral data captured by the CCD 108 may be used for deriving the spectral interference in the polarization sensitive interferometer 104.
- the spectral interferograms captured by the CCD 108 are real-valued components
- Fourier transform thereof may produce a complex conjugate artifact mirroring the desired true sample image about the zero-phase delay in the entire complex space. For some embodiments, this may lead to un- resolved positive and negative distances.
- an in vivo full-range complex Fourier domain OCT technique may be enlisted, such as the one described in U.S. Provisional Patent Application Ser. No. 60/826,998, filed Sept. 26, 2006 (incorporated herein by reference).
- a constant modulation frequency was introduced into the spatial spectral interferograms when the probe beam, provided by the probe 110, is scanned over the sample 112.
- the reference mirror 126 in the reference arm 102 may be mounted onto a linear Piezo-translation stage 150.
- the Piezo- translation stage 150 may be configured to move the reference mirror 126 at a constant velocity across the B scan (i.e., x direction scan). Accordingly, in these embodiments, a constant modulation frequency fo is introduced into equation 4,
- Equation 5 Equation 5
- t is the timing of the x scanner when it scans over the B scan
- Hubert transformation may be first applied to equation 5 against t to obtain the analytic function
- a final Fourier transform against 2z may be used to retrieve the depth information of the sample 112 with the complex conjugate mirror terms removed from the output plane.
- the images formed by the vertically and horizontally polarized beam components may be separated by 2d around the zero delay line that is determined by embodiments of the reference assembly setup as described herein.
- the modulation frequency introduced in the spectral interferograms may be set at 1.25 kHz with the CCD line rate set at 10 kHz.
- the introduced modulation frequency is one eighth of the data acquisition rate, the phenomenon on the spectral interference fringe washout at the CCD 108 may be ignored.
- the OCT images formed by two orthogonally polarized beam components may be computed at once from discrete points measured in the x direction that make up a data matrix.
- 1000 discrete points may be measured to form a data matrix of 1000 by 1024 elements (slice, B scan) from which the OCT images may be computed.
- post-processing may be performed.
- the average of A-lines may be subtracted from the raw data, followed by fixed pattern noise removal. This operation may effectively remove/minimize autocorrelation, self-cross correlation, and camera noise artifacts presented in the final images, which may improve the image quality.
- the subtracted spectral interferograms may be converted into the equal frequency space by use of spline interpolation.
- one of the images may be shifted by the fixed delay (e.g., 2.0 mm in embodiments described herein) to coincide with the other image.
- the A-scan signals of the object structure may be obtained in the form of /A H ,v(z)exp[/ ⁇ H y(z)]. From the magnitude of envelopes Av ⁇ z) and A H ⁇ z) of interferograms, the depth-resolved retardation ⁇ (z) may be expressed as:
- Equation 7 depth-resolved reflectivity R(z)
- Equation 8 Av(z) and A H (z) may be determined by the amplitude of the sample electrical field rather than their phase, measurements of retardation and reflectivity may be substantially accurate.
- the fast axis direction ⁇ (z) may be applicable to various embodiments and may be expressed as:
- Equation 9 The measured values for ⁇ (z) and ⁇ (z) may be unambiguous in the range of 0° to 90° and 0° to 180°, respectively.
- FIGURES 3(A) and 3(B) Illustrated in FIGURES 3(A) and 3(B), for example, are exemplary in vitro results captured from a chicken breast sample.
- FIGURE 3(A) shows the reconstructed in-depth OCT images formed by two orthogonal states of light by use of the standard FDOCT method.
- the disclosed reference assembly separates the images by about 2 mm centered at zero- delay line.
- the true image may be deteriorated by the overlapped complex conjugate image, which may result in the difficult or impossible separation of the two images formed by orthogonally polarized light.
- the images are reconstructed by use of the full range complex imaging technique, i.e., modulation of the spectral interferograms across the B scan at a constant frequency induced by the movement of the reference mirror.
- full range complex imaging may be achieved with the complex conjugate image almost completely removed, where the image formed by the vertically polarized beam component is on the top of that formed by the horizontally polarized light.
- polarization sensitive imaging may be possible by simply shifting one of the images with a known amount of optical delay in the depth direction (in this example, about 2.0 mm).
- the full range complex imaging technique may enable the OCT images formed by two orthogonal polarization beam components to be separated and centered about the zero delay line where the sensitivity of the spectrometer is highest, the sensitivity-fall ing-off characteristic of the system may continue to affect the final determination of the retardation values of the sample.
- Accurate determination of retardation using equation 7 may require the system sensitivity at corresponding depths of two images formed by orthogonal polarization beam components to be equal.
- the measured curve of the system sensitivity may be used to retrospectively correct the OCT signals A(z) along the depth, applied for example to FIGURE 3(B), so that the sensitivity along the depth may be approximately flat.
- FIGURE 4 shows the system sensitivity (hollow circles on FIGURE 4) measured along the depth up to 3.5 mm, where it may be seen that for the imaging depth near to the zero delay line, the sensitivity is about 100 dB and this value drops to about 85 dB at the depth of 3.5 mm.
- a correction curve may be obtained, c(z), by a polynomial curve fitting to a 2 nd order that may then normalized resulting in c(z) having a maximum value of unity (solid curve on FIGURE 4).
- the corrected OCT signals along the depth may be calculated by:
- one of the images may be shifted to coincide the two images formed by orthogonal polarization beam components, respectively, to obtain the depth- resolved retardation ⁇ (z) and reflectivity R(z) (i.e., equations 7 and 8, respectively), without affecting the optic axis calculation of equation 9.
- a method comprising forming a horizontally-polarized light beam from a first reference light beam; forming a vertically-polarized light beam from a second reference light beam; forming an interference signal from the horizontally- and vertically-polarized light beams; and modulating the interference signal at a constant modulation frequency.
- FIGURE 5 shows exemplary final images resulting from various embodiments.
- a simple algorithm was used to find the surface of the sample, and then a signal band strip was isolated which stands for an optical length of about 2 mm from the surface.
- the signals in the other areas were forced to zero as there contains no useful OCT signals for calculation of the polarization sensitive images.
- FIGURE 5(A) gives the structural OCT image of the interrogated sample with the polarization artifact diminished.
- the straight-line shown is an artifact from the zero delay line that is a residue of the autocorrelation signal that was not completely removed by the pre-processing method. The presence of this artifact, however, may have a minimum effect on the interpretation of the sample.
- FIGURE 5(B) is a cumulative round-trip phase retardation image.
- the phase retardation map of the sample is gray-scale coded from 0° to 90°. Because of the cumulative nature, it may appear to have periodic changes along the depth. The non-uniformity of the periodic phase appearance is not an artifact, but rather, may reveal certain organization information of the tissue components that make up the sample.
- the fast axis orientation image in FIGURE 5(C) shows abrupt phase changes that are due to the phase wrapping effect because of the measured range lying between 0° and 180°.
- the separation distance between the two images formed by orthogonally polarized beam components and determined by the reference assembly design may limit the imaging depth that the system may achieve.
- the separation distance For imaging most biological tissues, which are usually opaque (e.g., retina and skin, and also thin samples, such as human cornea), the separation distance of about 2 mm may be sufficient.
- the OCT signals may rapidly fall into the noise floor.
- the separation distance may be increased by adjusting the reference assembly to increase the imaging depth.
- a relatively large depth of focus of the focusing sampling beam volume may be required to keep the system sensitivity as uniform as possible along the depth.
- an objective lens with a focal length of about 100 mm giving a depth of focus of about 4.2 mm may be enlisted. This, however, may reduce the lateral imaging resolution that may be achieved.
- an axicon lens may be used in the sample arm to focus the beam into the sample and to give an improved depth of focus while keeping the lateral resolution high.
- FIGURE 6 Another exemplary SD-PSOCT system including a reference assembly 602 is illustrated in FIGURE 6.
- the exemplary reference assembly 602 may include one or more of various features described herein for various other embodiments.
- the illustrated embodiment similarly to the embodiment illustrated in FIGURE 1 , may be configured to operatively couple to an interferometer 604 and a spectrometer 606 including a single charge-coupled device camera (CCD) 608 for spectral interferogram detection to perform polarization sensitive imaging.
- CCD charge-coupled device camera
- the reference assembly 602 may be configured to separate the OCT images formed by orthogonal polarization channels at the output plane of Fourier space by introducing a constant frequency modulation in the spatial varying spectral interferograms when a probe beam, provided by a probe 610, is swept over a sample 612. Because the orthogonally polarized beam components travel in the identical beam path in the spectrometer 606, the distortion between spectra may be avoided. Further, because the beam components travel together, calibration of the beams to achieve pixel-by-pixel correspondence may be unnecessary. [0051] Still referring to FIGURE 6, the exemplary reference assembly
- the 602 may be configured to receive light from a light source 614.
- the light from the light source 614 may be collimated by one or more lenses 616 and vertically polarized by a linear polarizer 618.
- the light may then be introduced to a polarization sensitive low coherence interferometer 604 including a beam splitter 620.
- the vertically-polarized reference light leaving the interferometer may then be delivered to the reference assembly and to the sample 612.
- the reference light delivered to the reference assembly 602 may be split by a non-polarization beam splitter 622 into two substantially equal beams.
- One of the beams (“first beam”) may pass through a quarter wave plate 624.
- the quarter wave plate 624 may have, in some embodiments, a fast axis oriented at 45°.
- the first beam may then be bent by mirror 626 and reflected by a mobile reference mirror 650.
- the first beam may then pass again through the quarter wave plate 624 to become a horizontally polarized beam.
- the other beam (“second beam”) may be bent by mirrors 628 and 654 and reflected back by the mobile mirror 650 to join with the horizontally polarized beam.
- the joined beam may then be routed back into the interferometer 604.
- the reference assembly 602 may include a neutral density filter 630 placed in one or both of the first and second beam paths for ensuring that the light for the two orthogonally polarized beam components (horizontal and vertical) impinged onto the CCD 608 has substantially equal strength.
- sample light may also be delivered to the sample 612 (hereinafter “sample light”).
- sample light may pass through a quarter wave plate 632, and then may be delivered to the sample 612 by an X galvanometer scanner 634 and an achromatic lens 636.
- the sample light may be introduced to the sample by the probe 610.
- the reference light returning from the reference arm 602 and the sample light returning from the sample 610 may be recombined through a lens 652 and coupled into a single mode fiber 638 for introduction to the spectrometer 606.
- a polarization controller 640 may be used in the single mode fiber 638 for fine tuning the polarization states of light in the spectrometer system.
- the spectrometer 606 may be variously configured to be suitable for the intended purpose(s).
- the spectrometer 606 may comprise one or more of a collimator 642, a polarization-insensitive transmission grating 644, and an achromatic focusing lens 646.
- the spectral data from the CCD 608 may then be outputted to a computer 648 for access by a user.
- FIGURE 7 illustrated is a depiction of the relationship between the focused sample beam and the effective positions for the mobile reference mirror.
- equation 5 may be written as:
- I(k) ⁇ H ⁇ z)cos[2kz + 2 ⁇ fy + ⁇ H ⁇ z)] + A ⁇ z)cos[2kz-2 ⁇ fy + ⁇ ⁇ ⁇ z)]
- Equation 10 Hubert transformation may be first applied to equation 10 against t to obtain the analytic function:
- ⁇ k) ⁇ H (:)exp[i[2k: ⁇ 2 ⁇ f Q t + ⁇ H (i) ⁇
- a final Fourier transform against 2z may be used to retrieve the depth information of the sample with the complex conjugate mirror terms removed from the output plane.
- the images formed by the vertically and horizontally polarized beam components may be separated by and mirrored around the zero delay line that is determined by embodiments of the reference assembly setup as described herein.
- Various embodiments of SD-PSOCT systems including a reference assembly may be capable of providing polarization sensitive imaging of various samples including, for example, various biological samples. Illustrated in FIGURES 8(A) and 8(B), for example, are exemplary in vitro results captured from a fresh chicken breast sample.
- FIGURE 8(A) shows the reconstructed in-depth OCT images formed by two orthogonal states of light by use of the standard FDOCT method.
- the disclosed reference arm assembly separates the images mirrored by the zero-delay line.
- the true image may be deteriorated by the overlapped complex conjugate image, which may result in the difficult or impossible separation of the two images formed by orthogonally polarized light.
- the images are reconstructed by use of the full range complex imaging technique, i.e., modulation of the spectral interferograms across the B scan at a constant frequency induced by the movement of the reference mirror.
- polarization sensitive imaging may be possible by simply flipping one of the images with respect to the zero delay line in the depth direction.
- the analytic function of the spectral interference of both orthogonally polarized lights as described in equation 6 and equation 11 may also be obtained by a Fourier filtering technique, rather than the Hubert transformation described herein.
- the Fourier filtering technique may consist of Fourier transforming the spatial varying intereference signals, applying a windowing function in the Fourier space, and then inverse Fourier transforming back into spatial varying interference signals.
- FIGURE 9 shows various images captured from a proximal nail fold of a human volunteer using an exemplary embodiment of SD-PSOCT.
- the scan covers an area of about 4 mm wide by about 1.3 mm deep (scaled by an average refractive index of about 1.35).
- the imaging rate was set at 10 frames per second with each frame covering 1000 A scans.
- the intensity image shown in FIGURE 9(A) clearly delineates the epidermal (E) and dermal (D) regions near the nail fold, cuticle (C), nail plate (P), nail bed (B), and the lower half (LH) of the nail plate.
- the phase retardation image in FIGURE 9(B) shows the highly birefhngent property of the nail bed that gives banded structures as the amount of phase retardation wraps several times from 0° to 90°.
- the transition from white to black may indicate that the light experiences a change in phase retardation angle with increasing depth due to the presence of birefringent tissue.
- the dermis region near to the nail fold appears homogeneous in the intensity image, whereas it shows structural appearance in the retardation image, indicating that the collagen fiber organization in this region is inhomogeneous.
- the same tissue properties observed in the retardation image are also evident in the fast axis orientation image as shown in FIGURE 9(C).
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Optics & Photonics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
Embodiments include apparatuses and methods for spectral domain polarization sensitive optical coherence tomography including a reference assembly for detection of the polarization sensitive spectral interferograms formed by vertically and horizontally polarized beam components. Interference signals between the reference and sample beams may be modulated at a constant frequency.
Description
METHOD AND APPARATUS FOR LOCALIZED POLARIZATION SENSITIVE IMAGING
CROSS-REFERENCE TO RELATED APPLICATIONS
[0001] The present application claims priority to U.S. Provisional
Patent Application No. 60/942,278, entitled "Method and Apparatus for Localized Polarization Sensitive Imaging," filed June 6, 2007, the entire disclosure of which is hereby incorporated by reference in its entirety for all purposes, except for those sections, if any, that are inconsistent with this specification.
TECHNICAL FIELD
[0002] Embodiments relate generally to the field of biomedical imaging, specifically to methods, apparatuses, and systems associated with optical coherence tomographic imaging.
BACKGROUND
[0003] Optical coherence tomography (OCT) is a non-invasive, three- dimensional (3-D) imaging technique, capable of optical ranging within a highly scattering sample, such as biological tissue. The development of polarization sensitive (PS) OCT has enabled OCT to sense the birefringence properties of biological tissue that are not provided by conventional OCT. Spectral domain (SD) OCT may provide further advantages because of the significant sensitivity and imaging speed advantages over time domain OCT. Current techniques, however, may not be capable of producing images of sufficient depth or clarity. Furthermore, conventional methods may require burdensome calibration steps for producing the desired images.
BRIEF DESCRIPTION OF THE DRAWINGS
[0004] Embodiments will be readily understood by the detailed description in conjunction with the accompanying drawings. Embodiments are illustrated by way of example and not by way of limitation in the figures of the accompanying drawings.
[0005] FIGURE 1 illustrates an SD-PSOCT system in accordance with various embodiments;
[0006] FIGURE 2 illustrates an exemplary reference mirror in accordance with various embodiments;
[0007] FIGURES 3(A) and 3(B) are images of a chicken breast acquired using an exemplary embodiment of an SD-PSOCT system in accordance with various embodiments;
[0008] FIGURE 4 is a plot of SD-PSOCT sensitivity against imaging depth in accordance with various embodiments;
[0009] FIGURES 5(A), 5(B), and 5(C) are images of a chicken breast acquired using an exemplary embodiment of an SD-PSOCT system in accordance with various embodiments;
[0010] FIGURE 6 illustrates another SD-PSOCT system in accordance with various embodiments;
[0011] FIGURE 7 illustrates another exemplary reference mirror in accordance with various embodiments;
[0012] FIGURES 8(A) and 8(B) are images of a chicken breast acquired using an exemplary embodiment of an SD-PSOCT system in accordance with various embodiments; and
[0013] FIGURES 9(A), 9(B), and 9(C) are images of a human nail fold acquired using an exemplary embodiment of an SD-PSOCT system in accordance with various embodiments.
DETAILED DESCRIPTION OF DISCLOSED EMBODIMENTS
[0014] In the following detailed description, reference is made to the accompanying drawings which form a part hereof and in which is shown by way of illustration various embodiments. It is to be understood that other embodiments may be utilized and structural or logical changes may be made without departing from the scope. Therefore, the following detailed description is not to be taken in a limiting sense, and the scope of embodiments is defined by the appended claims and their equivalents.
[0015] Various operations may be described as multiple discrete operations in turn, in a manner that may be helpful in understanding the disclosed embodiments; however, the order of description should not be construed to imply that these operations are order dependent.
[0016] The description may use perspective-based descriptions such as up/down, back/front, and top/bottom. Such descriptions are merely used to facilitate the discussion and are not intended to restrict the application of the disclosed embodiments.
[0017] The description may use the phrases "in an embodiment," or
"in embodiments," which may each refer to one or more of the same or different embodiments. Furthermore, the terms "comprising," "including," "having," and the like, as used with respect to disclosed embodiments, are synonymous.
[0018] A phrase in the form of "A/B" or in the form "A and/or B" means
"(A), (B), or (A and B)." A phrase in the form "at least one of A, B and C" means "(A), (B), (C), (A and B), (A and C), (B and C) or (A, B and C)." A phrase in the form "at least one of A, B, and C" means "(A), (B), (C), (A and B), (A and C), (B and C), or (A, B and C)." A phrase in the form "(A) B" means "(B) or (A B)," that is, A is optional.
[0019] In various embodiments, methods, apparatuses, and systems for imaging, such as biomedical imaging, are provided. In exemplary embodiments, a computing system may be endowed with one or more
components of the disclosed apparatuses and/or systems and may be employed to perform one or more methods as disclosed herein.
[0020] In an embodiment, there is provided an apparatus comprising a first light path configured to receive a first reference light beam and form, from the first reference light beam, a horizontally-polarized light beam; a second light path configured to receive a second reference light beam and form, from the second reference light beam, a vertically-polarized light beam; and a reference mirror configured to modulate, at a constant modulation frequency, an interference signal formed by the horizontally- and vertically- polarized light beams and the sample beam.
[0021] According to various embodiments, a spectral domain polarization sensitive optical coherence tomography (SD-PSOCT) system having a single-camera configuration may comprise a reference arm assembly (hereinafter "reference assembly") for polarization sensitive imaging. Although other applications are envisioned, the disclosed embodiments may be particularly suitable for in vitro and in vivo imaging.
[0022] An exemplary SD-PSOCT system including a reference assembly 102 is illustrated in FIGURE 1. As illustrated, the exemplary reference assembly 102 may be configured to operatively couple to an interferometer 104 and a spectrometer 106 including a single charge- coupled device camera (CCD) 108 for spectral interferogram detection to perform polarization sensitive imaging. According to various embodiments, the reference assembly 102 may be configured to separate by a fixed distance OCT images from orthogonal polarization channels at the output plane of Fourier space by introducing a constant frequency modulation in the spatial varying spectral interferograms when a probe beam, provided by a probe 110, is swept over a sample 112. Because the orthogonally polarized beam components travel in the identical beam path in the spectrometer 106, the distortion between spectra may be avoided. Further, because the beam components travel together, calibration of the beams to achieve pixel-by- pixel correspondence may be unnecessary.
[0023] Still referring to FIGURE 1 , the exemplary reference assembly
102 may be configured to receive light from a light source 114 such as, for example, a broadband superluminescent diode (SLD). Although any light source suitable for the purpose may be enlisted for light source 114, an exemplary light source may be an SLD with a central wavelength of about 850 nm and full-width half-maximum (FWHM) bandwidth of about 30 nm, yielding a measured axial resolution of about 12 pm in air. For various other embodiments, the light source 114 may be a swept laser source, rather than the SLD. In embodiments of an SD-PSOCT system including a swept laser source, the spectrometer 106 may be replaced by an ultra-fast photodiode detector. For various embodiments, an SD-PSOCT system including a swept laser source as the light source 114 and photodiode for detection may perform substantially similarly to embodiments including an SLD as the light source 114 and the spectrometer 106 for detection.
[0024] Light from the light source 114 may be collimated by one or more lenses 116. According to various embodiments, the diameter of the collimated light may be about 3 mm.
[0025] Light from the light source 114 may be vertically polarized by a linear polarizer 118 and introduced to a polarization sensitive low coherence interferometer 104 including a beam splitter 120. The vertically-polarized light leaving the interferometer 104 (hereinafter "reference light") may then be delivered to the reference assembly 102 and the sample 112.
[0026] According to various embodiments, the reference light delivered to the reference assembly 102 may be split by a non-polarization beam splitter 122 into two substantially equal beams. One of the beams ("first beam") may pass through a quarter wave plate 124. The quarter wave plate 124 may have a fast axis oriented at 45°. The first beam may then be reflected by a staggered reference mirror 126. The first beam may then pass again through plate 124 to become a horizontally polarized beam. The other beam ("second beam") split by the non-polarization beam splitter 122 may be bent by another mirror 128 and reflected back by the staggered
reference mirror 126 to join with the horizontally polarized beam. The joined beam may then be routed back into the interferometer 104.
[0027] In various embodiments, an SD-PSOCT system described herein may include one or more optical fiber couplers in addition to or instead of a beam splitter. So, in various embodiments, one or both of the beam splitters 120 and 122 may instead be an optical fiber coupler, depending on the particular desired application.
[0028] The reference assembly 102 may include a neutral density filter 130 placed in one or both of the first and second beam paths. The neutral density filter 130 may ensure that the light for two orthogonally polarized beam components (horizontal and vertical) impinged onto the CCD 108 has substantially equal strength. With such design, the reference assembly 102 may provide a fixed optical pathlength difference (OPD) between these two orthogonally polarized lights. The fixed OPD may be adjusted according to the needs for different applications. For example, according to various embodiments, the OPD between the vertically and horizontally polarized lights may be adjusted to about 2.0 mm. Because the orthogonally polarized beam components travel in the identical beam path in the spectrometer 106, the OCT images formed by the two orthogonally polarized beam components at the output plane of Fourier space may be separated after processing the mixed spectral interferograms captured by the CCD 108 (as described more fully herein).
[0029] As noted above, light from the polarization sensitive low coherence interferometer 104 may also be delivered to the sample 112 (hereinafter "sample light"). The sample light may pass through a quarter wave plate 132. The quarter wave plate 132 may have a fast axis oriented at 45° to transform the sample light into circularly polarized light. The sample light may then be delivered to the sample 112 by an X galvanometer scanner 134 and an achromatic lens 136. According to various embodiments, the achromatic lens 136 may include any focal length suitable
for the application. In some embodiments, for example, the focal length may be about 100 mm.
[0030] The sample light may be introduced to the sample by the probe
110. The theoretical lateral imaging resolution, depth of focus, power, and other parameters of the probe beam may be variously adjusted depending on the application. In an embodiment, for example, the probe beam may have a theoretical lateral imaging resolution of about 28 μm, a depth of focus of about 4.2 mm, and a power of about 0.6 mW.
[0031] The reference light returning from the reference arm 102 and the sample light returning from sample 110 may be recombined through a lens 152 and coupled into a single mode fiber 138 for introduction to the spectrometer 106. According to various embodiments, a polarization controller 140 may be used in the single mode fiber 138 for fine tuning the polarization states of light in the spectrometer system.
[0032] The spectrometer 106 may be variously configured to be suitable for the intended purpose(s). According to the illustrated embodiment, the spectrometer 106 comprises a collimator 142, a polarization-insensitive transmission grating 144, and an achromatic focusing lens 146. According to an exemplary embodiment, the spectrometer 106 may comprise a 30 mm focal length collimator 142, a 1200 lines/mm polarization-insensitive transmission grating 144, and an achromatic focusing lens 146 with 190 mm focal length. The CCD 108 may comprise a line scan charge-coupled device camera for capturing both orthogonally polarized beam components (i.e., horizontal and vertical components) in parallel. Although the spectrometer 106 may include any CCD suitable for the purpose, in an embodiment, the CCD 108 may comprise one having 1024 pixels, each about 14x14 microns in size and 12 bits in digital depth, and capable of a 53 kHz line rate.
[0033] According to various embodiments, the propagating beam paths for vertically and horizontally polarized light components are substantially identical in the spectrometer 106, and their spectra captured by
the identical CCD pixel arrays in parallel. This may be particularly advantageous in that careful calibration of the spectral components on the CCD array may not be necessary. For various embodiments, the spectrometer 106 may be designed to have an imaging depth being the full depth in the Fourier space used for imaging (e.g., approximately 8.2 mm in air), rather than half of the space used in conventional SD-PSOCT systems.
[0034] According to various embodiments, the CCD 108 may be variously adjusted for imaging. For various embodiments, for example, the CCD integration time may be set at about 100 microseconds for imaging. The spectral data from the CCD 108 (e.g., 1024 pixels, A scan) may then be outputted to a computer 148 for access by a user. According to an exemplary embodiment, the image data is downloaded to the computer 148 using CAMERALINK™ and a high-speed frame grabber board (for example, PC1 1428, commercially available from National Instruments, USA). For an exemplary embodiment, the CCD line scan rate may be 10 kHz , the imaging speed may be limited by the output power of the light source used (about 5 mW), and the signal sensitivity of 100 dB may be measured at z=+0.5 mm and dropped to about 85 dB at z=+3.5 mm when the camera integration time is set at about 100 microseconds. In this example, the light energy in the reference arm 102 may be adjusted such that the maximum intensity on the linear CCD array detector achieved approximately 40% of the camera's saturation level for each polarization beam component, and the remaining 20% of the dynamic range of the camera may be available to fill up the polarization interference fringes.
[0035] In an embodiment, there is provided a system comprising a light source; a reference assembly including: a first light path configured to receive a first reference light beam from the light source and form, from the first reference light beam, a horizontally-polarized light beam; a second light path configured to receive a second reference light beam from the light source and form, from the second reference light beam, a vertically-polarized light beam; and a signal modulating device configured to modulate an interference signal formed by the horizontally- and vertically-polarized light
beams and a sample light beam to give a constant modulation frequency; a sampling assembly configured to provide light from the light source to a sample to obtain the sample light beam; and a spectrometer configured to detect a spectral interferogram formed by the horizontally-polarized light beam, the vertically-polarized light beam, and the sample light beam.
[0036] In various embodiments, the modulation of the interferogram may also be provided by scanner 134 and reference mirror 150 may be maintained in one position. In an exemplary embodiment, scanner 134 may be configured such that the input signal is scanned with an offset reference to the pivot point.
[0037] Turning now to FIGURE 2, with continuing reference to
FIGURE 1 , illustrated is a depiction of the relationship between the focused sample beam and the effective positions for the reference mirror 126. Assuming that the coordinate origin along the depth direction, z, is at the focal spot, the optical fields in the reference arm 102 for the vertically and horizontally polarized lights may be expressed as:
(hereinafter "equation 1 ") where H and V refer to the signals from the horizontal and vertical polarization channels, R is the reflectivity in the reference arm 102, k is the wavenumber, and d is the distance of the reference mirror 126 position relative to the focal spot. In the sample arm (i.e., X galvanometer scanner 134, an achromatic lens 136, and probe 110), the optical fields that are decomposed into the vertically and horizontally polarized light may be expressed as:
(hereinafter "equation 2") where SH,V(Z) represents the strengths of backscattering light at optical depth z in the sample in the respective polarization direction. Thus, considering that the vertically and horizontally
polarized beam components are independent, the spectral interferogram detected at the CCD 108 may be described as:
(hereinafter "equation 3"). Because of the linear relationship between the contributions from different optical depth z from within a sample to equation 2, the problem of spectral interferogram detected by the CCD 108 at wavenumber k may be bundled and simplified by ignoring the terms that do not contribute to the useful signals, without loss of generosity, as:
!{k ) - An i% ) cos[2*C t + d) 4 φJf (2)}+ 4.(z)co$[2lt- ™ «0 * ?1 C41
(hereinafter "equation 4") wherein AH,γ(z) represents the magnitude of envelopes of interferograms formed by the respective polarization components and ΨH,Y(Z) ιs tne random phase term that is due to the optical properties of the sample 112. Thus, the interferograms formed by the vertical and horizontal polarization light beams may be added together at the CCD 108.
[0038] According to various embodiments, the sample target may be placed near the zero phase delay line where the sensitivity may be the highest. This may be particularly desirable if the CCD 108 is insufficiently capable of digitizing the spectral components due to decreased sensitivity along the imaging depth. According to embodiments, OCT images formed by the two orthogonal polarization states may be placed side by side approximately centered at the zero delay line. For some embodiments, the images are placed about 2 mm apart.
[0039] According to various embodiments, spectral data captured by the CCD 108 may be used for deriving the spectral interference in the polarization sensitive interferometer 104. In embodiments wherein the spectral interferograms captured by the CCD 108 are real-valued components, Fourier transform thereof may produce a complex conjugate artifact mirroring the desired true sample image about the zero-phase delay in the entire complex space. For some embodiments, this may lead to un-
resolved positive and negative distances. Accordingly, for various embodiments, an in vivo full-range complex Fourier domain OCT technique may be enlisted, such as the one described in U.S. Provisional Patent Application Ser. No. 60/826,998, filed Sept. 26, 2006 (incorporated herein by reference). In that disclosure, a constant modulation frequency was introduced into the spatial spectral interferograms when the probe beam, provided by the probe 110, is scanned over the sample 112. According to various embodiments, the reference mirror 126 in the reference arm 102 may be mounted onto a linear Piezo-translation stage 150. The Piezo- translation stage 150 may be configured to move the reference mirror 126 at a constant velocity across the B scan (i.e., x direction scan). Accordingly, in these embodiments, a constant modulation frequency fo is introduced into equation 4,
(hereinafter "equation 5") where t is the timing of the x scanner when it scans over the B scan, and t=0 is the start of the B scan. To eliminate the complex conjugate mirror images from the output plane of Fourier space, Hubert transformation may be first applied to equation 5 against t to obtain the analytic function,
hk) *> AH(z)txpH2k{x ut) + 2π fo + φu(z)]}
(hereinafter "equation 6").
[0040] A final Fourier transform against 2z may be used to retrieve the depth information of the sample 112 with the complex conjugate mirror terms removed from the output plane. The images formed by the vertically and horizontally polarized beam components may be separated by 2d around the zero delay line that is determined by embodiments of the reference assembly setup as described herein.
[0041] According to various embodiments, the modulation frequency introduced in the spectral interferograms may be set at 1.25 kHz with the
CCD line rate set at 10 kHz. For these embodiments, because the introduced modulation frequency is one eighth of the data acquisition rate, the phenomenon on the spectral interference fringe washout at the CCD 108 may be ignored. For various embodiments, the OCT images formed by two orthogonally polarized beam components may be computed at once from discrete points measured in the x direction that make up a data matrix. In an exemplary embodiment, 1000 discrete points may be measured to form a data matrix of 1000 by 1024 elements (slice, B scan) from which the OCT images may be computed.
[0042] After the data has been acquired, post-processing may be performed. According to various embodiments, the average of A-lines may be subtracted from the raw data, followed by fixed pattern noise removal. This operation may effectively remove/minimize autocorrelation, self-cross correlation, and camera noise artifacts presented in the final images, which may improve the image quality. For various embodiments, the subtracted spectral interferograms may be converted into the equal frequency space by use of spline interpolation. In some embodiments, full range complex imaging processing may be applied to the processed data to remove the complex conjugate images presented in the Fourier output plane, leaving two OCT images side by side separated by 2d = 2mm as described herein. The two images may represent the images formed in the horizontal and vertical polarization channels, respectively. In still further embodiments, one of the images may be shifted by the fixed delay (e.g., 2.0 mm in embodiments described herein) to coincide with the other image. In the images, the A-scan signals of the object structure may be obtained in the form of /AH,v(z)exp[/ ΦHy(z)]. From the magnitude of envelopes Av{z) and AH{z) of interferograms, the depth-resolved retardation δ(z) may be expressed as:
(hereinafter "equation 8"). Since Av(z) and AH(z) may be determined by the amplitude of the sample electrical field rather than their phase, measurements of retardation and reflectivity may be substantially accurate. The fast axis direction θ(z) may be applicable to various embodiments and may be expressed as:
(hereinafter "equation 9"). The measured values for δ(z) and θ(z) may be unambiguous in the range of 0° to 90° and 0° to 180°, respectively.
[0043] Various embodiments of SD-PSOCT systems including a reference assembly may be capable of providing polarization sensitive imaging of various samples including, for example, various biological samples. Illustrated in FIGURES 3(A) and 3(B), for example, are exemplary in vitro results captured from a chicken breast sample. FIGURE 3(A) shows the reconstructed in-depth OCT images formed by two orthogonal states of light by use of the standard FDOCT method. As expected, the disclosed reference assembly separates the images by about 2 mm centered at zero- delay line. However, due to the real-valued function of spectral interferograms captured by the CCD camera, the true image may be deteriorated by the overlapped complex conjugate image, which may result in the difficult or impossible separation of the two images formed by orthogonally polarized light. Shown in FIGURE 3(B), however, the images are reconstructed by use of the full range complex imaging technique, i.e., modulation of the spectral interferograms across the B scan at a constant frequency induced by the movement of the reference mirror. As may be seen, full range complex imaging may be achieved with the complex conjugate image almost completely removed, where the image formed by the vertically polarized beam component is on the top of that formed by the
horizontally polarized light. Accordingly, for various embodiments, polarization sensitive imaging may be possible by simply shifting one of the images with a known amount of optical delay in the depth direction (in this example, about 2.0 mm).
[0044] Although the full range complex imaging technique may enable the OCT images formed by two orthogonal polarization beam components to be separated and centered about the zero delay line where the sensitivity of the spectrometer is highest, the sensitivity-fall ing-off characteristic of the system may continue to affect the final determination of the retardation values of the sample. Accurate determination of retardation using equation 7 may require the system sensitivity at corresponding depths of two images formed by orthogonal polarization beam components to be equal. For various embodiments, the measured curve of the system sensitivity may be used to retrospectively correct the OCT signals A(z) along the depth, applied for example to FIGURE 3(B), so that the sensitivity along the depth may be approximately flat. FIGURE 4 shows the system sensitivity (hollow circles on FIGURE 4) measured along the depth up to 3.5 mm, where it may be seen that for the imaging depth near to the zero delay line, the sensitivity is about 100 dB and this value drops to about 85 dB at the depth of 3.5 mm. From the measured data, a correction curve may be obtained, c(z), by a polynomial curve fitting to a 2nd order that may then normalized resulting in c(z) having a maximum value of unity (solid curve on FIGURE 4). Thereafter, in various embodiments, the corrected OCT signals along the depth may be calculated by:
[0045] According to various embodiments and as described herein, one of the images may be shifted to coincide the two images formed by orthogonal polarization beam components, respectively, to obtain the depth- resolved retardation δ(z) and reflectivity R(z) (i.e., equations 7 and 8, respectively), without affecting the optic axis calculation of equation 9.
[0046] In an embodiment, there is provided a method comprising forming a horizontally-polarized light beam from a first reference light beam; forming a vertically-polarized light beam from a second reference light beam; forming an interference signal from the horizontally- and vertically-polarized light beams; and modulating the interference signal at a constant modulation frequency.
[0047] FIGURE 5 shows exemplary final images resulting from various embodiments. For obtaining FIGURE 5, a simple algorithm was used to find the surface of the sample, and then a signal band strip was isolated which stands for an optical length of about 2 mm from the surface. The signals in the other areas were forced to zero as there contains no useful OCT signals for calculation of the polarization sensitive images. FIGURE 5(A) gives the structural OCT image of the interrogated sample with the polarization artifact diminished. The straight-line shown is an artifact from the zero delay line that is a residue of the autocorrelation signal that was not completely removed by the pre-processing method. The presence of this artifact, however, may have a minimum effect on the interpretation of the sample. FIGURE 5(B) is a cumulative round-trip phase retardation image. The phase retardation map of the sample is gray-scale coded from 0° to 90°. Because of the cumulative nature, it may appear to have periodic changes along the depth. The non-uniformity of the periodic phase appearance is not an artifact, but rather, may reveal certain organization information of the tissue components that make up the sample. The fast axis orientation image in FIGURE 5(C) shows abrupt phase changes that are due to the phase wrapping effect because of the measured range lying between 0° and 180°.
[0048] For some embodiments, the separation distance between the two images formed by orthogonally polarized beam components and determined by the reference assembly design may limit the imaging depth that the system may achieve. For imaging most biological tissues, which are usually opaque (e.g., retina and skin, and also thin samples, such as human cornea), the separation distance of about 2 mm may be sufficient. For
optical depths beyond about 2 mm, the OCT signals may rapidly fall into the noise floor. For relatively transparent biological samples, however, the separation distance may be increased by adjusting the reference assembly to increase the imaging depth.
[0049] For some embodiments wherein the two virtual reference mirrors are spaced about 2 mm apart, a relatively large depth of focus of the focusing sampling beam volume may be required to keep the system sensitivity as uniform as possible along the depth. In an embodiment, an objective lens with a focal length of about 100 mm giving a depth of focus of about 4.2 mm may be enlisted. This, however, may reduce the lateral imaging resolution that may be achieved. For various embodiments, an axicon lens may be used in the sample arm to focus the beam into the sample and to give an improved depth of focus while keeping the lateral resolution high.
[0050] Another exemplary SD-PSOCT system including a reference assembly 602 is illustrated in FIGURE 6. As illustrated, the exemplary reference assembly 602 may include one or more of various features described herein for various other embodiments. For example, the illustrated embodiment, similarly to the embodiment illustrated in FIGURE 1 , may be configured to operatively couple to an interferometer 604 and a spectrometer 606 including a single charge-coupled device camera (CCD) 608 for spectral interferogram detection to perform polarization sensitive imaging. According to various embodiments, the reference assembly 602 may be configured to separate the OCT images formed by orthogonal polarization channels at the output plane of Fourier space by introducing a constant frequency modulation in the spatial varying spectral interferograms when a probe beam, provided by a probe 610, is swept over a sample 612. Because the orthogonally polarized beam components travel in the identical beam path in the spectrometer 606, the distortion between spectra may be avoided. Further, because the beam components travel together, calibration of the beams to achieve pixel-by-pixel correspondence may be unnecessary.
[0051] Still referring to FIGURE 6, the exemplary reference assembly
602 may be configured to receive light from a light source 614. The light from the light source 614 may be collimated by one or more lenses 616 and vertically polarized by a linear polarizer 618. The light may then be introduced to a polarization sensitive low coherence interferometer 604 including a beam splitter 620. The vertically-polarized reference light leaving the interferometer may then be delivered to the reference assembly and to the sample 612.
[0052] According to various embodiments, the reference light delivered to the reference assembly 602 may be split by a non-polarization beam splitter 622 into two substantially equal beams. One of the beams ("first beam") may pass through a quarter wave plate 624. The quarter wave plate 624 may have, in some embodiments, a fast axis oriented at 45°. The first beam may then be bent by mirror 626 and reflected by a mobile reference mirror 650. The first beam may then pass again through the quarter wave plate 624 to become a horizontally polarized beam. The other beam ("second beam") may be bent by mirrors 628 and 654 and reflected back by the mobile mirror 650 to join with the horizontally polarized beam. The joined beam may then be routed back into the interferometer 604.
[0053] The reference assembly 602 may include a neutral density filter 630 placed in one or both of the first and second beam paths for ensuring that the light for the two orthogonally polarized beam components (horizontal and vertical) impinged onto the CCD 608 has substantially equal strength.
[0054] As noted above, light from the interferometer 604 may also be delivered to the sample 612 (hereinafter "sample light"). The sample light may pass through a quarter wave plate 632, and then may be delivered to the sample 612 by an X galvanometer scanner 634 and an achromatic lens 636. The sample light may be introduced to the sample by the probe 610.
[0055] The reference light returning from the reference arm 602 and the sample light returning from the sample 610 may be recombined through
a lens 652 and coupled into a single mode fiber 638 for introduction to the spectrometer 606. According to various embodiments, a polarization controller 640 may be used in the single mode fiber 638 for fine tuning the polarization states of light in the spectrometer system.
[0056] The spectrometer 606 may be variously configured to be suitable for the intended purpose(s). The spectrometer 606 may comprise one or more of a collimator 642, a polarization-insensitive transmission grating 644, and an achromatic focusing lens 646. The spectral data from the CCD 608 may then be outputted to a computer 648 for access by a user.
[0057] Turning now to FIGURE 7, illustrated is a depiction of the relationship between the focused sample beam and the effective positions for the mobile reference mirror. The same analysis described herein with reference to equations 1 to 4, in part or in whole, may be applied here, but with d=0. Consequently, equation 5 may be written as:
I(k) = ΛH{z)cos[2kz + 2πfy + φH{z)] + AΛz)cos[2kz-2πfy + φιΛz)]
(hereinafter "equation 10"). To eliminate the complex conjugate mirror images from the output plane of Fourier space, Hubert transformation may be first applied to equation 10 against t to obtain the analytic function:
ϊ{k) = ΛH (:)exp[i[2k: ÷ 2π fQt + φH(i)}}
+ A1. (z ) exp {i[-2/b -2,τfj + <pr { z )] )
(hereinafter "equation 11 ").
[0058] A final Fourier transform against 2z may be used to retrieve the depth information of the sample with the complex conjugate mirror terms removed from the output plane. The images formed by the vertically and horizontally polarized beam components may be separated by and mirrored around the zero delay line that is determined by embodiments of the reference assembly setup as described herein.
[0059] Various embodiments of SD-PSOCT systems including a reference assembly may be capable of providing polarization sensitive imaging of various samples including, for example, various biological samples. Illustrated in FIGURES 8(A) and 8(B), for example, are exemplary in vitro results captured from a fresh chicken breast sample. FIGURE 8(A) shows the reconstructed in-depth OCT images formed by two orthogonal states of light by use of the standard FDOCT method. As expected, the disclosed reference arm assembly separates the images mirrored by the zero-delay line. However, due to the real-valued function of spectral interferograms captured by the CCD camera, the true image may be deteriorated by the overlapped complex conjugate image, which may result in the difficult or impossible separation of the two images formed by orthogonally polarized light. Shown in FIGURE 8(B), however, the images are reconstructed by use of the full range complex imaging technique, i.e., modulation of the spectral interferograms across the B scan at a constant frequency induced by the movement of the reference mirror. As may be seen, full range complex imaging may be achieved with the complex conjugate image almost completely removed, where the image formed by the vertically polarized beam component is in the negative space of the Fourier output plane, and that formed by the horizontally polarized light is in the positive space. Accordingly, for various embodiments, polarization sensitive imaging may be possible by simply flipping one of the images with respect to the zero delay line in the depth direction.
[0060] In various embodiments, the analytic function of the spectral interference of both orthogonally polarized lights as described in equation 6 and equation 11 may also be obtained by a Fourier filtering technique, rather than the Hubert transformation described herein. The Fourier filtering technique may consist of Fourier transforming the spatial varying intereference signals, applying a windowing function in the Fourier space, and then inverse Fourier transforming back into spatial varying interference signals.
[0061] Various embodiments are suitable for use in in vivo applications. FIGURE 9 shows various images captured from a proximal nail fold of a human volunteer using an exemplary embodiment of SD-PSOCT. For the embodiment, the scan covers an area of about 4 mm wide by about 1.3 mm deep (scaled by an average refractive index of about 1.35). The imaging rate was set at 10 frames per second with each frame covering 1000 A scans. The intensity image shown in FIGURE 9(A) clearly delineates the epidermal (E) and dermal (D) regions near the nail fold, cuticle (C), nail plate (P), nail bed (B), and the lower half (LH) of the nail plate. The phase retardation image in FIGURE 9(B) shows the highly birefhngent property of the nail bed that gives banded structures as the amount of phase retardation wraps several times from 0° to 90°. At the epidermal-dermal boundary (EDB) of the nail fold, the transition from white to black may indicate that the light experiences a change in phase retardation angle with increasing depth due to the presence of birefringent tissue. The dermis region near to the nail fold appears homogeneous in the intensity image, whereas it shows structural appearance in the retardation image, indicating that the collagen fiber organization in this region is inhomogeneous. The same tissue properties observed in the retardation image are also evident in the fast axis orientation image as shown in FIGURE 9(C).
[0062] Although certain embodiments have been illustrated and described herein for purposes of description of the preferred embodiment, it will be appreciated by those of ordinary skill in the art that a wide variety of alternate and/or equivalent embodiments or implementations calculated to achieve the same purposes may be substituted for the embodiments shown and described without departing from the scope. Those with skill in the art will readily appreciate that embodiments may be implemented in a very wide variety of ways. This application is intended to cover any adaptations or variations of the embodiments discussed herein. Therefore, it is manifestly intended that embodiments be limited only by the claims and the equivalents thereof.
Claims
1. An apparatus comprising: a first light path configured to receive a first reference light beam and form, from the first reference light beam, a horizontally-polarized light beam; a second light path configured to receive a second reference light beam and form, from the second reference light beam, a vertically-polarized light beam; and a reference mirror configured to modulate, at a constant modulation frequency, an interference signal formed by the horizontally- and vertically- polarized light beams and a sample light beam.
2. The apparatus of claim 1 , further comprising a beam splitter configured to receive a source light beam and split the source light beam into the first reference light beam and the second reference light beam.
3. The apparatus of claim 2, wherein the beam splitter is configured to receive a vertically-polarized and/or a horizontally-polarized source light beam.
4. The apparatus of claim 1 , wherein the first light path and the second light path are configured to form the horizontally-polarized light beam and the vertically-polarized light beam with a fixed optical pathlength difference.
5. The apparatus of claim 4, further comprising a scanner configured for sweeping the sample light beam across a sample to provide the constant modulation frequency.
6. The apparatus of claim 1 , wherein the first light path includes a quarter wave plate.
7. The apparatus of claim 6, wherein the quarter wave plate includes a fast axis oriented at about 45°.
8. The apparatus of claim 1 , wherein the reference mirror comprises a mobile reference mirror.
9. The apparatus of claim 8, wherein the reference mirror is mounted on a linear moving stage configured to move the reference mirror at a constant velocity across a B scan.
10. The apparatus of claim 1 , wherein the first light path and/or the second light path includes a neutral density filter.
11. The apparatus of claim 1 , wherein the first light path and the second light path, and the reference mirror form a reference assembly of a spectral domain optical coherence tomography system.
12. A system comprising: a light source; a reference assembly including: a first light path configured to receive a first reference light beam from the light source and form, from the first reference light beam, a horizontally-polarized light beam; a second light path configured to receive a second reference light beam from the light source and form, from the second reference light beam, a vertically-polarized light beam; and a signal modulating device configured to modulate, at a constant modulation frequency, an interference signal formed by the horizontally- and vertically-polarized light beams and a sample light beam; a sampling assembly configured to provide light from the light source to a sample to obtain the sample light beam; and a spectrometer configured to detect a spectral interferogram formed by the horizontally-polarized light beam, the vertically-polarized light beam, and the sample light beam.
13. The system of claim 12, wherein the first light path and the second light path are configured to form the horizontally-polarized light beam and the vertically-polarized light beam with a fixed optical pathlength difference.
14. The system of claim 12, wherein the spectrometer includes a line- scan charge-coupled device camera.
15. The system of claim 12, wherein the signal modulating device comprises a reference mirror.
16. The system of claim 12, wherein the signal modulating device comprises a scanner configured for sweeping the sample light beam across a sample to provide the constant modulation frequency.
17. A method comprising: forming a horizontally-polarized light beam from a first reference light beam; forming a vertically-polarized light beam from a second reference light beam; forming an interference signal from the horizontally- and vertically- polarized light beams; and modulating the interference signal at a constant modulation frequency.
18. The method of claim 17, further comprising providing a source light beam to a sample to obtain a sample light beam, and detecting a spectral interferogram formed by the horizontally-polarized light beam, the vertically- polarized light beam, and the sample light beam.
19. The method of claim 18, wherein modulating the interference signal at a constant modulation frequency comprises offsetting an input signal at a scanner in the sample beam.
20. The method of claim 17, wherein modulating the interference signal at a constant modulation comprises moving a reference mirror at a constant frequency.
21. The method of claim 17, further comprising splitting a source light beam to form the first reference light beam and the second reference light beam.
22. The method of claim 17, wherein the horizontally-polarized light beam and the vertically-polarized light beam are formed with a fixed optical pathlength difference.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US12/602,822 US9086264B2 (en) | 2007-06-06 | 2008-06-06 | Polarization sensitive spectral domain OCT using an interference signal modulated at a constant frequency and a two-path reference arm with a single reference mirror |
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US94227807P | 2007-06-06 | 2007-06-06 | |
| US60/942,278 | 2007-06-06 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| WO2008154349A2 true WO2008154349A2 (en) | 2008-12-18 |
| WO2008154349A3 WO2008154349A3 (en) | 2009-03-19 |
Family
ID=40130443
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| PCT/US2008/066055 Ceased WO2008154349A2 (en) | 2007-06-06 | 2008-06-06 | Method and apparatus for localized polarization sensitive imaging |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US9086264B2 (en) |
| WO (1) | WO2008154349A2 (en) |
Families Citing this family (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110237999A1 (en) * | 2010-03-19 | 2011-09-29 | Avedro Inc. | Systems and methods for applying and monitoring eye therapy |
| JP2015534056A (en) | 2012-09-24 | 2015-11-26 | トルネード スペクトラル システムズ,インコーポレイテッド | On-chip multi-function spectrometer with a single detector array |
| US9279660B2 (en) * | 2013-05-01 | 2016-03-08 | Canon Kabushiki Kaisha | Method and apparatus for processing polarization data of polarization sensitive optical coherence tomography |
| US9885557B2 (en) | 2013-10-11 | 2018-02-06 | Case Western Reserve University | Polarization sensitive optical coherence tomography using multiple polarization sensitive semiconductor optical amplifiers |
| US9291500B2 (en) * | 2014-01-29 | 2016-03-22 | Raytheon Company | Configurable combination spectrometer and polarizer |
| CN104297885B (en) * | 2014-11-10 | 2016-12-07 | 长春理工大学 | A kind of reference mirror is bonding and self-calibrating device |
| AU2015367283B2 (en) | 2014-12-14 | 2021-05-20 | Alcon Inc. | Multichannel optical receivers |
| WO2016187675A1 (en) * | 2015-05-28 | 2016-12-01 | Cylite Pty Ltd | High resolution 3-d spectral domain optical imaging apparatus and method |
| CA3012490A1 (en) | 2016-02-05 | 2017-08-10 | Synaptive Medical (Barbados) Inc. | System and method for providing surgical guidance based on polarization-sensitive optical coherence tomography |
| CN105559756B (en) * | 2016-02-05 | 2019-11-15 | 浙江大学 | Microangiography method and system based on full spatial modulation spectrum segmentation and angle compounding |
| CN108956533A (en) * | 2018-07-13 | 2018-12-07 | 福州大学 | A kind of OCT conjugation mirror image removal device and method for eliminating polychrome error |
| US11473897B2 (en) * | 2018-10-12 | 2022-10-18 | The General Hospital Corporation | Method and apparatus for measuring depth-resolved tissue birefringence using single input state polarization sensitive optical coherence tomography |
| US11175126B2 (en) | 2019-04-08 | 2021-11-16 | Canon U.S.A., Inc. | Automated polarization control |
| EP3730897A1 (en) * | 2019-04-25 | 2020-10-28 | Nokia Technologies Oy | Apparatus, systems and methods for detecting light |
| US12277731B2 (en) * | 2022-07-21 | 2025-04-15 | Canon U.S.A., Inc. | Methods and systems for system self-diagnosis |
Family Cites Families (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB8903725D0 (en) * | 1989-02-18 | 1989-04-05 | Cambridge Consultants | Coherent tracking sensor |
| US5861952A (en) | 1992-11-16 | 1999-01-19 | Canon Kabushiki Kaisha | Optical inspection method and apparatus including intensity modulation of a light beam and detection of light scattered at an inspection position |
| WO1998043068A1 (en) * | 1997-03-26 | 1998-10-01 | Kowa Company, Ltd. | Optical measuring instrument |
| US20050140981A1 (en) * | 2002-04-18 | 2005-06-30 | Rudolf Waelti | Measurement of optical properties |
| JP2003337007A (en) | 2002-05-21 | 2003-11-28 | Canon Inc | Interference device |
| JP4505807B2 (en) * | 2004-08-09 | 2010-07-21 | 国立大学法人 筑波大学 | Multiplexed spectral interferometric optical coherence tomography |
| JP4563130B2 (en) | 2004-10-04 | 2010-10-13 | 株式会社トプコン | Optical image measuring device |
| EP1785690A1 (en) * | 2005-11-10 | 2007-05-16 | Haag-Streit Ag | Method and device for determining an object s geometric characteristics |
| JP5149196B2 (en) * | 2005-12-06 | 2013-02-20 | カール ツァイス メディテック アクチエンゲゼルシャフト | Sample measurement by interferometry |
| WO2008091961A2 (en) * | 2007-01-23 | 2008-07-31 | Volcano Corporation | Optical coherence tomography implementation |
| JP5339828B2 (en) * | 2007-10-04 | 2013-11-13 | キヤノン株式会社 | Optical coherence tomography apparatus and optical coherence tomography method |
-
2008
- 2008-06-06 WO PCT/US2008/066055 patent/WO2008154349A2/en not_active Ceased
- 2008-06-06 US US12/602,822 patent/US9086264B2/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| US20100182609A1 (en) | 2010-07-22 |
| WO2008154349A3 (en) | 2009-03-19 |
| US9086264B2 (en) | 2015-07-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9086264B2 (en) | Polarization sensitive spectral domain OCT using an interference signal modulated at a constant frequency and a two-path reference arm with a single reference mirror | |
| EP2341823B1 (en) | System and method for providing full jones matrix-based analysis to determine non-depolarizing polarization parameters using optical frequency domain imaging | |
| US7697145B2 (en) | System for fourier domain optical coherence tomography | |
| CN101214145B (en) | Frequency domain optical coherence tomography method and system with large detection depth | |
| US8081316B2 (en) | Process, system and software arrangement for determining at least one location in a sample using an optical coherence tomography | |
| EP1971848B1 (en) | Systems and methods for generating data based on one or more spectrally-encoded endoscopy techniques | |
| US8514404B2 (en) | Multiple path interferometer and method | |
| US8605287B2 (en) | Extended range imaging | |
| US20070038040A1 (en) | Arrangements, systems and methods capable of providing spectral-domain polarization-sensitive optical coherence tomography | |
| CN201139554Y (en) | Frequency domain optical coherence tomography device with large detection depth | |
| US7812961B2 (en) | Method, apparatus, and program for processing tomographic images | |
| GB2407155A (en) | Spectral interferometry method and apparatus | |
| JP2019512086A (en) | Apparatus and method for high speed and long depth range imaging using optical coherence tomography | |
| Fan et al. | Spectral domain polarization sensitive optical coherence tomography achieved by single camera detection | |
| CN103792192A (en) | Polarization frequency domain optical coherence tomography system based on single detector | |
| KR101919416B1 (en) | Device for polarization sensitive optical coherence tomography on one-unit detector | |
| JP2014206433A (en) | Optical tomographic image-capturing device | |
| Sharma et al. | Human ex-vivo oral tissue imaging using spectral domain polarization sensitive optical coherence tomography | |
| Baradaran et al. | Visible light polarization-sensitive optical coherence tomography with balanced detection | |
| JP2008151734A (en) | Optical tomographic imaging method, apparatus and program, and optical tomographic imaging system | |
| Gelikonov et al. | Suppression of image artifacts in the spectral-domain optical coherence tomography | |
| Gelikonov et al. | New approaches in broadband fiber-optical interferometry for optical coherent tomography | |
| Arbab et al. | Terahertz imaging and spectroscopic measurements in the presence of scattering for biophotonics applications | |
| Harris et al. | Polarization-sensitive fast THz-TDS scanner in a portable form-factor | |
| Jiao et al. | Characterization of the polarization properties of biological tissues with fiber-based Mueller-matrix optical coherence tomography |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| 121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 08770288 Country of ref document: EP Kind code of ref document: A2 |
|
| WWE | Wipo information: entry into national phase |
Ref document number: 12602822 Country of ref document: US |
|
| NENP | Non-entry into the national phase |
Ref country code: DE |
|
| 122 | Ep: pct application non-entry in european phase |
Ref document number: 08770288 Country of ref document: EP Kind code of ref document: A2 |
