[go: up one dir, main page]

Next Article in Journal
A Review of Visual-LiDAR Fusion based Simultaneous Localization and Mapping
Previous Article in Journal
OmniSCV: An Omnidirectional Synthetic Image Generator for Computer Vision
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optical Interferometric Fringe Pattern-Incorporated Spectrum Calibration Technique for Enhanced Sensitivity of Spectral Domain Optical Coherence Tomography

1
School of Electronics Engineering, College of IT Engineering, Kyungpook National University, 80, Daehak-ro, Buk-gu, Daegu 41566, Korea
2
Department of Biomedical Engineering, College of Engineering, Kyungil University, 50, Gamasil-gil, Hayang-eup, Gyeongsan-si, Gyeongsangbuk-do 38428, Korea
3
Department of Autonomous Robot Engineering, College of Smart Engineering, Kyungil University, 50, Gamasil-gil, Hayang-eup, Gyeongsan-si, Gyeongsangbuk-do 38428, Korea
4
Department of Nuclear Energy Convergence, Kyungil University, 50, Gamasil-gil, Hayang-eup, Gyeongsan-si, Gyeongsangbuk-do 38428, Korea
5
School of Architecture and Civil Engineering, Kyungpook National University, 80, Daehak-ro, Buk-gu, Daegu 41566, Korea
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Sensors 2020, 20(7), 2067; https://doi.org/10.3390/s20072067
Submission received: 11 February 2020 / Revised: 30 March 2020 / Accepted: 2 April 2020 / Published: 7 April 2020
(This article belongs to the Section Optical Sensors)
Figure 1
<p>System configuration of the SD-OCT system. Abbreviations: BLS, broadband laser source; C, collimator; DG, diffraction grating; DC, dispersion compensator; E, eye; FC, fiber coupler; GM, Galvano mirror; L, lens; LSC, line scan camera; M, mirror; PC, polarization controller; RSL, retinal scanning lens.</p> ">
Figure 2
<p>Graphical representation of the fundamental steps of the software-based intensity detection and resampling algorithm developed for the spectral calibration.</p> ">
Figure 3
<p>Flow diagram of the fundamental steps performed for dispersion of the system. The red color solid arrows illustrate switching image states (changing of images), while compensating dispersion using decoupled hardware and software methods. The blue color solid arrows illustrate switching image states, while changing the focusing position.</p> ">
Figure 4
<p>Evaluations of the point spread function (PSF)-based spectral interferogram analysis. (<b>a</b>) The depth information table of each PSF position, showing the respective number of points carrying maximum intensity, (<b>b</b>) correlated pixel numbers (pixel positions) at each number of points with maximum intensity prior to interpolation, and (<b>c</b>) the averaged graphical plot of the eight respective graphs in <a href="#sensors-20-02067-f004" class="html-fig">Figure 4</a>b.</p> ">
Figure 5
<p>Quantitative representation of system sensitivity assessments. (<b>a</b>) Depth-dependent system intensity fall-off and (<b>b</b>) depth-dependent sensitivity roll-off (signal-to-noise ratio, SNR) representation. Each PSF peak of <a href="#sensors-20-02067-f004" class="html-fig">Figure 4</a>a corresponds to the correlated pixel numbers of <a href="#sensors-20-02067-f004" class="html-fig">Figure 4</a>b.</p> ">
Figure 6
<p>Full width at half maximum (FWHM) comparison between before (without) and after (with) spectral calibrated depth-dependent PSFs.</p> ">
Figure 7
<p>Spectrum calibrated two dimensional OCT images acquired prior to dispersion compensation. (<b>a</b>) 2D-OCT images of the infrared (IR)-detection card at top and bottom depth levels. (<b>b</b>) 2D-OCT images of five glass cover slips at top and bottom depth levels. Red color square boxes represent a region of interest with non-corrected dispersion. Horizontal scale bar: 1 mm, vertical scale bar: 300 μm.</p> ">
Figure 8
<p>Spectrum calibrated two dimensional OCT images acquired after performing dispersion compensation. (<b>a</b>) 2D-OCT images of IR-detection card at top and bottom depth levels. (<b>b</b>) 2D-OCT images of five glass cover slips at top and bottom depth levels. Green color square boxes represent a region of interest with corrected dispersion. Horizontal scale bar: 1 mm, vertical scale bar: 300 μm.</p> ">
Figure 9
<p>Image comparison with and without spectral calibrated in vivo human retina. (<b>a</b>) In vivo human retinal image before (without) spectral calibration; (<b>b</b>) shows magnified view of red box. (<b>c</b>) In vivo human retinal image after (with) spectral calibration; (<b>d</b>) shows magnified view of green box. (<b>a</b>) Horizontal scale bar: 1 mm, vertical scale bar: 300 μm, (<b>b</b>) horizontal scale bar: 500 μm, vertical scale bar: 150 μm.</p> ">
Figure 10
<p>OCT image verification of <span class="html-italic">in vivo</span> rat posterior chamber along with B-scans corresponding to multiple positions. (<b>a</b>) Before calibration, (<b>b</b>) after calibration, (<b>c</b>) averaged posterior vitreous opacifications. The horizontal and vertical scale bars of P1–P4 (position 1–position 4) are 1 mm and 700 μm, respectively.</p> ">
Figure 11
<p>The OCT image verification of <span class="html-italic">in vivo</span> rat anterior vitreous opacification along with B-scans corresponds to multiple positions. (<b>a</b>) 3D enface visualization of anterior vitreous opacification, (<b>b</b>) before calibration, and (<b>c</b>) after calibration. The horizontal and vertical scale bars of (<b>b</b>) are 500 μm and 300 μm, respectively.</p> ">
Versions Notes

Abstract

:
Depth-visualizing sensitivity can be degraded due to imperfect optical alignment and non-equidistant distribution of optical signals in the pixel array, which requires a measurement of the re-sampling process. To enhance this depth-visualizing sensitivity, reference and sample arm-channeled spectra corresponding to different depths using mirrors were obtained to calibrate the spectrum sampling prior to Fourier transformation. During the process, eight interferogram patterns corresponding to point spread function (PSF) signals at eight optical path length differences were acquired. To calibrate the spectrum, generated intensity points of the original interferogram were re-indexed towards a maximum intensity range, and these interferogram re-indexing points were employed to generate a new lookup table. The entire software-based process consists of eight consecutive steps. Experimental results revealed that the proposed method can achieve images with a high depth-visualizing sensitivity. Furthermore, the results validate the proposed method as a rapidly performable spectral calibration technique, and the real-time images acquired using our technique confirm the simplicity and applicability of the method to existing optical coherence tomography (OCT) systems. The sensitivity roll-off prior to the spectral calibration was measured as 28 dB and it was halved after the calibration process.

1. Introduction

Optical coherence tomography (OCT) enables micrometer-resolution biomedical imaging to be performed in real time [1,2,3]. Due to its capability to perform non-invasive inspection, OCT has been extensively utilized and featured in numerous scientific research and clinical studies [4,5,6,7,8,9,10]. In addition to its medical applications, spectral domain OCT (SD-OCT) based field inspections have also been applied to industrial and agricultural materials [11,12,13,14,15,16]. One issue with SD-OCT, however, is that the sensitivity of its depth-dependent signals tend to be weaker in deeper imaging regions [17,18,19], as a cause of slight technical errors associated with the optical spectrometer.
An optical spectrometer should, therefore, be designed according to a specific center wavelength and a full width at half maximum (FWHM) bandwidth; this will allow it to obtain wavenumber measurements that are uniformly distributed among an array of pixels. If the interference fringe signal of a SD-OCT is subjected to a discrete Fourier transformation, then its performance can be optimized once the Fourier transforming data points are equally distributed [20]. The Fourier transformation operates in the frequency domain, while the camera pixel array produces a linear representation of the spectrum in terms of its wavelength, not frequency. In OCT, the final image is a convolution of the system point spread function (PSF) and the object being imaged. Therefore, an object that is a delta function would provide a direct measurement of the point-spread limited resolution [21]. If the data are not linearized by being inserted into equidistant frequency slots, the depth-dependent PSF can be distorted, and each pixel of the line-scan camera of the optical spectrometer can then integrate with a different spectral width [22,23,24]. Spectral calibration methods have, therefore, always been a primary concern. Numerous approaches were demonstrated by using complex numerical data resampling methods and using bulky optical hardware methods due to unavailability of supremely well-designed spectrometers [25,26,27,28,29].
Resampling schemes involving a graphics processing unit (GPU) have been frequently proposed in relevant studies to convert the wavelength domain to the wavenumber domain enhancing the data acquisition speed [20,30]. In SD-OCT, a technique for linearization of spectral fringe signals has been demonstrated by correcting the ν = c/λ relation with both software and optical hardware elements by extracting accumulative phase evolution of the interference fringes [30,31,32,33]. In the aforementioned techniques, the interference fringe patterns were processed using linear interpolation with fast Fourier transformation (FFT), cubic spline interpolation with FFT, non-uniform discrete Fourier transforms (NDFT), and non-uniform fast Fourier transforms (NFFT). Furthermore, compressed sensing methods that are based on wavenumber domain linearizing techniques have been demonstrated to simplify complex numerical methods [18,26]. Some studies have also focused on optical hardware optimizations by using a wavelength-scanning filter technique, which provides lookup tables instead of interpolating non-linear data [34,35]; other studies have used prism implementations to achieve real-time imaging [17,25,36].
Simultaneous correction of non-linear k-sampling and dispersion mismatch between sample and reference arms was demonstrated by extracting two calibration vectors for the numerical resampling for wavenumber-linearization and phase correction for the dispersion mismatch [37]. Multiple complex numerical steps with linear and non-linear terms were assessed to determine the position of the peak in A-lines and peak broadening distortion, respectively. Additionally, a full-range visualizing technique was integrated along with spectral and dispersion calibration, which enhances the OCT providing a large depth range [38]. To eliminate the necessity of linearization, a dual interferometer-based (master and slave interferometers) approach was developed to operate with a sensitivity similar to SD-OCT providing signal from a single point in depth [39,40] and a similarly complex master slave OCT method was performed to process the interferometer spectrum without Fourier transformation. Semi-automatic spectral mapping is another technique, which was acquired using dual spectrometers, where a set of calibration interference fringes has to be acquired from both spectrometers [41,42]. In these two aforementioned methods [41,42], non-linearity in the unwrapped phase information was sampled by generating linear vectors. For these newly generated linear vectors, a fractional index through interpolation corresponding to its position in the original phase vector can be identified. Thereafter, these fractional indices can be used to resample all subsequent measurements in order to obtain wavenumber linearization samplings. Once wavenumber linearization is completed, A-lines of the interference signal are acquired to find peak axial positions followed by Gaussian filtering, fast Fourier transformation, and further signal processing [37].
In this paper, we present a simply implementable 8-step alternative spectral calibration method to overcome the major drawback of SD-OCT imaging systems. The method only requires raw data of optical PSF resolved interferograms, which is then fed into a customized computational algorithm. The SD-OCT hybrid algorithm encapsulates binarization, intensity detection, interpolation, and indexation, which are performed in real-time. The acquired raw data of multiple adjacent depth positions were well utilized in the developed algorithm to simply functionalize for any broadband laser and line scan camera of SD-OCT. Therefore, the method proposed here can be easily integrated into current SD-OCT systems to overcome spectral broadening limitations and to obtain an enhanced depth-dependent sensitivity fall-off.

2. Materials and Methods

2.1. Description of Optical Coherence Tomography (OCT) System Configuration

Built upon spectral interferometry, SD-OCT uses an optical interference signal consisting of a sample arm beam and a reference arm beam, with an exactly-matched optical path length difference. The experimental setup of our study featured a retinal inspection OCT (shown in Figure 1) that was developed for the direct implementation of the proposed scheme. The system was operated with a broadband light source (EXS210022-02, Exalos, Switzerland) that had a center wavelength of 840 nm and a bandwidth of 50 nm. A 50:50 optical fiber coupler (TW850R5A2, Thorlabs, Newton, NJ, USA) was used to split the optical light beams into a sample arm and a reference arm, and the backscattered optical light beams were recorded by a laboratory-customized spectrometer, containing a diffraction grating with 1800 lines/mm and a line scanning camera (2048 pixels, E2V, UK). During image acquisition, 500 B-scans were acquired in 6 seconds. Each B-scan consisted of 500 A-scans. The detailed optical configurations and system specifications of the retinal OCT instrumentation including retinal scanning lens configuration have been previously reported elsewhere [43].

2.2. Spectral Calibration Algorithm

Figure 2 illustrates the 8 simple fundamental steps of the software-based intensity detection and resampling algorithm, developed for spectral calibration. We obtained PSF incorporated spectrums by using two mirrors as reference and sample arms, starting from 100 pixel depth range, towards 800 pixel depth range, at eight optical path length differences (shown in step one). Step two emphasized the extraction of the interferogram pattern using PSF, which was utilized as baseline data. We then applied an intensity-based peak detection filter to each interferogram to select the maximum intensity points (shown in step 3) and the selected points were later binarized (shown in step 4). Then, the pixel position corresponding to these each binarized maximum intensity point was graphically rearranged and plotted as pixel number (position) against maximum intensity point (shown in step 5). An analogous procedure was repeated for the remaining seven depth levels (200 pixel level to 800 pixel depth range) and the program executed continuously until the condition of step 6 was satisfied. Thus, the algorithm repetitively executed for eight PSF trials. The spectral graphs including spectrally encoded multiple interference fringes from multiple depths were streamlined to enhance the imaging depth and signal-to-noise ratio (SNR) [32]. When PSF-based data acquisition was completed across all of the depth levels, the algorithm executed step 7, which was the cubic spline interpolation [44] step of the maximum intensity points. As the number of maximum intensity points of each depth level was dissimilar, all intensity points had to be interpolated to the total number of pixels of the OCT line scan camera (total number of pixels in current study: 2048). Therefore, during this step 7, the maximum intensity points were resampled using nearest adjacent intensity points coinciding with sampled data at the interpolated nodes. Once this interpolation was completed, the dimensions of each depth level pixel array were similar to each other. In the current experiment, eight PSF trials were performed and therefore, eight curves illustrating correlated pixel numbers (pixel positions) at each number of points with maximum intensity were arranged, summed up, and averaged to generate a single representative curve (shown in step 8). Next, a linear curve was sampled with similar number of data samples and plotted along with the interpolated data, which was linearly fitted later through indexing the non-linear sample points. Finally, a lookup table was generated based on these linearizing index values to rescale real-time raw data enhancing the OCT image. Since the key baseline data required for the developed algorithm operation was the extracted intensity information of PSF based interferogram patterns, the algorithm could be functionalized well despite the FWHM bandwidth of optical broadband laser sources (≥50 nm (FWHM)).
Moreover, cubic-spline interpolation is a complex numerical calibration technique, where the image pixels are stored in a compute unified device architecture (CUDA) array and bound to a designated region in a texture memory. Resampled data points are generated by interpolating the raw fringe data points, which eventually provides evenly distributed pixel information. From the entire array of pixels, raw data were extracted from 0~800 pixel depth range due to weakened fringe information (less intensity) of the extreme corner beyond the 800 pixel range. Therefore, 8 PSF measurements (n = 8) were acquired with a pixel interval of 100 pixels, which provided sufficient raw data of the entire depth range with least missing depth information that eventually led to a precise spectral calibration. Fewer PSF measurements (n < 8) will reduce the preciseness of data with lack of depth information, which has negative impact on distinguishing deep layers of OCT B-scans. Although non-uniform data point distribution was corrected, a slight dispersion was identified on the image due to mismatch between the reference and sample propagating paths, where the light is specularly reflecting from the mirror surface; and secondly, due to depth-dependent dispersion of the sample with several layers of coverslip glass. Thus primarily, the spectral calibration was performed to enhance the depth dependent sensitivity, and secondly the dispersion mismatch was compensated using both hardware (dispersion compensating glass-rods) and numerical operations based software methods. Figure 3 illustrates the flow diagram of the fundamental steps performed after the spectral calibration to further compensate dispersion of the system. To decouple two technical contributions (spectral calibration and dispersion compensation), a single frequency channeled spectrum with no dispersion compensation was acquired in the original state of the image initially. This data stream contains different sampling information between sample and reference arm, which is the main cause for the image dispersion. Next to confirm the exact optical arm with dispersion, BK-7 glass rods (dispersion compensating glass rods) were placed on reference arm. A sample consists of multiple (five) glass cover slips representing multiple depth levels was placed on sample arm and the image was visualized. Raw data of the corresponding image was extracted and arranged to examine dispersion variation. This set of raw data was arranged and inserted to the spectrum calibration algorithm. Instead of using the phase information of the interferogram, fundamental index intensity of the interferogram was varied to compensate dispersion. A numerical multiplication factor with an intensity weight (numerical weight with hundredth decimal point starting from 1) was utilized in the algorithm to match the dispersion of all the corresponding depth levels. The utilization of this numerical multiplication factor results in adjustments of the index weight and sampling of the entire array, which interpolates an index to calibrate the spectrum. The exact numerical factor, which calibrates the spectrum is verified after successful interpolation attempts, and the entire interpolated raw index array is Fourier transformed and log scaled for the image visualization. Simultaneously, the index intensity can be continuously varied to enhance the image quality until the spectrum calibration is satisfactorily accomplished. The performed BK-7 based reference arm dispersion compensation can be successfully confirmed if the acquired image reveals a similar thickness in all depth levels of the sample. When a thickness difference is identified at deeper regions, dispersion of the sample arm has to be corrected using additional BK-7 glass rods until a sufficient dispersion compensation with similar thickness for all depth levels is obtained. Subsequently, a new set of raw data can be arranged and index intensity has to be varied through the numerical factor of the spectrum calibration algorithm to match further dispersion of the entire depth range. The utilization of this numerical factor results in adjustments of the index weight and sampling of the entire array, which simply leads to an index interpolation to compensate dispersion. As a result of this decoupled operation of the spectrum calibration technique and dispersion compensation technique, the dispersion mismatch between reference and sample arms can be compensated by cancelling the frequency-dependent real and the imaginary parts of the spectral fringe pattern, which are used to construct the complex analytical representation of the spectral fringe pattern.

2.3. In Vivo Experimental Procedure

We performed in vivo SD-OCT imaging of the human retina with and without adoption of the developed spectral calibration technique. All human experiments were performed in accordance with the guidelines of the Institutional Animal and Human Care and Use Committee of Kyungpook National University (No. KNU-2018-0100). We captured two image sets, each consists of pair of two dimensional OCT (2D-OCT) images of an exact location of human retina. We further utilized our developed spectral calibration method to confirm the progressive OCT image quality improvement by examining in vivo Sprague Dawley rat posterior chambers. The preparation of the Sprague Dawley rat specimens and the anesthetic procedure were performed according to a previously reported method [45]. The in vivo rat eye specimens were repeatedly scanned to obtain a three-dimensional (3D) enface visualization of their posterior vitreous opacifications. The calibrated volumetric OCT images were averaged for speckle noise reduction and cross-sections were extracted corresponding to four particular positions (position 1–position 4). In addition to the retinal posterior vitreous assessments, in vivo anterior vitreous visualizations were obtained also from in vivo Sprague Dawley rat specimens. B-scans were obtained before and after performing the developed spectral calibration technique. Unlike posterior segment imaging systems, anterior segment scanners require illumination of the eye with a focused optical beam. Therefore, the retinal scanning lens (RSC) of schematic 1 (Figure 1) can be removed temporarily to fulfill this purpose.

3. Results and Discussion

3.1. Analysis of Point Spread Function (PSF) Incorporated Spectral Calibration

The initial experimental procedure of the proposed scheme is the acquisition of spectral interferograms corresponds to multiple PSF trials. The table illustrated in Figure 4a provides depth information for each PSF position starting from the 100th-pixel to 800th-pixel, along with the number of points carrying the maximum intensity in each corresponding interferogram. Since a non-calibrated spectrometer broadens the spectral width along with the depth of the PSF, this broadening effect increases the number of points carrying the maximum intensity. To overcome this dissimilarity, we interpolated each set of maximum intensity points, and converted them to 2048 (the number of pixels in the OCT camera pixel array).
Figure 4b shows the correlated pixel numbers (pixel positions) at each number of points with maximum intensity prior to interpolation acquired during the assessment of each depth-dependent PSF position. The solid black color curve of Figure 4c emphasizes the summated and averaged curve of eight corresponding curves of Figure 4b illustrating pixel position information over each averaged maximum intensity point. From the results, it is seen that the practical point distribution is no longer linear owing to the spectral errors that occur in spectrometer. Therefore, we deliberately generated a calibrated domain, where the non-linear points were re-indexed and linearized. The purple dashed curve shows the deliberately generated linear domain obtained after curve fitting, which is initiated from the lowest averaged maximum intensity point to uppermost point. We used the closest method for a linear pattern by selecting the closest value between the measured averaged intensity point and the linear point, where the index values were later utilized in a real-time image acquisition program.

3.2. Depth-Dependent Sensitivity Assessments after Spectral Calibration

Figure 5a shows the depth dependent intensity fall-off of the system, where we obtained PSFs before and after applying the spectral calibration technique. The black dashed curve shows the intensity roll-off before calibration verifying the broadening effect of the spectral width. The multicolor solid sharp peaks show the intensity roll-off after the spectrum was calibrated. The comparison shown in Figure 5a clearly reveals that the averaged fall-off curve acquired after the spectral calibration exhibits much less decay than the averaged non-calibrated roll-off, across the whole imaging depth range. Similarly, the depth-dependent sensitivity roll-off emphasizing signal-to-noise ratio (SNR) graphs acquired before and after spectral calibration are plotted in Figure 5b. A significant sensitivity roll-off (28 dB) occurs prior to the spectral calibration, while the intensity fall-off was 14 dB after calibration.

3.3. Full Width at Half Maximum (FWHM) Evaluation of Depth-Dependent PSF

We further examined and compared the FWHM between before and after spectral calibrated depth-dependent PSFs as shown in Figure 6. Our primary objective of this evaluation was to gain a better understanding of the depth visualizing sensitivity and pixel-resolution of deep structures, as the FWHM has a direct impact on the pixel-resolution incorporated axial resolution. Thus, the evaluated FWHM of the peaks are presented as pixel units, which have a direct relationship with thickness. According to the graphical representation, Figure 6 suggests that the trend of the non-calibrated PSF width increasingly fluctuates in the deeper regions, while the width of the spectrally-calibrated PSF remains relatively constant across the whole imaging depth range. Henceforth, this comparative evaluation demonstrates that the developed spectral calibration method effectively improves depth-visualizing sensitivity and resolution in deeper imaging regions for SD-OCT systems. Therefore, examined depth-dependent PSF intensity fluctuation and sensitivity roll-off of our method emphasized a higher sensitivity with less noise residual confirming the potential merits over previously reported conventional techniques [17,26].

3.4. Implementation of Dispersion Correction

After performing a spectral calibration process, initial examination of the developed spectrum optimization was assessed using infrared (IR) detection card and a set of five glass cover slips. Figure 7a,b depict 2D-OCT images of the IR-detection card and set of five glass cover slips acquired after spectrum calibration.
Although the individual layer thickness of both samples at the top depth level was similarly distributed at the bottom depth level, sub-surface structural properties were not identified distinguishably due to a dispersion mismatch between the sample and reference arms. BK-7 glass rods were placed on both reference and sample arms, and the number of rods was interchanged until a certain amount of dispersion was compensated for. To overcome the remaining dispersion mismatch between two arms, the fundamental index intensity of the raw index of interference signal was varied using a numerical multiplication factor with hundredth decimal point. The acquired spectral calibration and dispersion correction can be confirmed in Figure 8a,b. The results reveal that the layer thickness in both samples was similarly distributed at both depth levels along with significantly distinguishable sub-surface structural properties.

3.5. In Vivo OCT Image Verification of the Developed Spectral Calibration Method

In Figure 9, two image sets, each consists of pair of 2D-OCT images of an exact similar location of the human retina are presented. Figure 9a illustrates the 2D-OCT image acquired prior to the calibration and Figure 9c shows the image acquired after employing the calibration. Figure 9b,d are the zoomed images of the region of interest indicated on Figure 9a,c. Noticeably, the 2D-OCT image shown in Figure 9a emphasizes a less depth-visualizing sensitivity and less resolution in deep structures (in the case where no spectral calibration was performed) compared to Figure 9c, where the qualitative representation was significantly enhanced after calibration. Thus, in terms of SD-OCT based future clinical imaging applications, the developed technique shows good stability and reproducibility for any type of SD-OCT system.
The developed spectral calibration method was further utilized to confirm the progressive OCT image quality improvement by examining in vivo Sprague Dawley rat posterior vitreous opacifications. Figure 10a,b exhibit structural differences and image enhancements between the images taken before and after the calibration process revealing a higher visibility of blood vessels in 10b. The calibrated 3D-OCT images were further averaged for speckle noise reduction (Figure 10c) and B-scans corresponding to four particular positions (p1~p4) were extracted from both Figure 10a,c as indicated in Figure 10.
Figure 11a depicts the enhanced 3D-OCT enface visualization of the anterior vitreous opacification, along with B-scans that were obtained before and after performing the spectral calibration technique. In particular, the effective depth visualizing sensitivity enhancement of the anterior segment including cornea, iris, and lens region can be clearly observed after the utilization of our spectral calibration technique. Hence, the results of Section 3.4 qualitatively verified that the developed calibration method performed over the entire spectrum remarkably enhanced the visualization of microstructural features of transparent mediums, such as anterior and posterior vitreous owing to the enhanced depth-visualizing sensitivity and image resolution of deep structures.
Similar to previously reported methods [37,38,39,40,41,42], the developed fringe pattern-incorporated spectrum calibration method is applicable to enhance the sensitivity of the SD-OCT system and to compensate the dispersion mismatch as a simultaneous approach. Indeed, it relies on intensity information of the interferometric fringe signals, and the adjustable intensity weight (multiplication factor) confirmed the simple functionality and lower complexity of the proposed method compared to conventional methods. Although the current approach is closely related to the methods reported in [37,38,39] in the introduction, rapid calibration of the spectrum along with dispersion compensation using a single inferometer (without any other complex hardware) is still an important difference. Moreover, the approach described in this manuscript minimized the necessity of complex numerical calculations, which dramatically reduces computational time. Since the developed method is particularly relevant for SD-OCT systems consisting of broadband lasers with any center wavelength at any FWHM bandwidth and capable of real-time rectification of spectrum non-linarites and dispersion mismatch, extended ophthalmological applications can be anticipated.

4. Conclusions

In this study, we introduce an optical PSF resolved interferogram-incorporated spectral calibration technique as an uncomplicated alternative technique to overcome the depth-dependent sensitivity loss commonly experienced in SD-OCT systems. The validation of a comprehensive 8 step-based spectral calibration technique was applied for in vivo human retinal imaging and in vivo retinal and corneal imaging of Sprague Dawley rats. The experimental results, including quantifications and qualitative in vivo assessments, sufficiently confirmed the simplicity and accuracy of our technique along with an axial resolution of 8 μm (measured in air) and maximum imaging depth of 2.6 mm (in air). Moreover, the comparison between non-calibrated and calibrated spectral information revealed its methodological precision, and highlighted that it was able to decrease unfavorable sensitivity fall-off and significantly increase depth-visualizing sensitivity and image resolution in deep cross-sectional regions. Therefore, this calibration method is suitable to use with any standard SD-OCT system, as a simple incorporative spectral calibration technique.

Author Contributions

S.H., Y.H., and R.E.W. conducted the experiment, analyzed and drafted the manuscript. D.J., J.L. (Jaeyul Lee), J.L. (Junsoo Lee), and H.J. designed the experimental setup with software and hardware components. D.-E.L., M.J. and J.K. designed the study and critically revised and finalized the intellectual contents of manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government, MSIT (No. 2018R1A5A1025137).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wojtkowski, M.; Srinivasan, V.J.; Ko, T.H.; Fujimoto, J.G.; Kowalczyk, A.; Duker, J.S. Ultrahigh-resolution, high-speed, fourier domain optical coherence tomography and methods for dispersion compensation. Opt. Express 2004, 12, 2404–2422. [Google Scholar] [CrossRef] [PubMed]
  2. Wojtkowski, M.; Srinivasan, V.; Fujimoto, J.G.; Ko, T.; Schuman, J.S.; Kowalczyk, A.; Duker, J.S. Three-dimensional retinal imaging with high-speed ultrahigh-resolution optical coherence tomography. Ophthalmology 2005, 112, 1734–1746. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Cense, B.; Nassif, N.A.; Chen, T.C.; Pierce, M.C.; Yun, S.-H.; Park, B.H.; Bouma, B.E.; Tearney, G.J.; de Boer, J.F. Ultrahigh-resolution high-speed retinal imaging using spectral-domain optical coherence tomography. Opt. Express 2004, 12, 2435–2447. [Google Scholar] [CrossRef] [PubMed]
  4. Fujimoto, J.G.; Brezinski, M.E.; Tearney, G.J.; Boppart, S.A.; Bouma, B.; Hee, M.R.; Southern, J.F.; Swanson, E.A. Optical biopsy and imaging using optical coherence tomography. Nat. Med. 1995, 1, 970–972. [Google Scholar] [CrossRef] [PubMed]
  5. Zysk, A.M.; Nguyen, F.T.; Oldenburg, A.L.; Marks, D.L.; Boppart, S.A. Optical coherence tomography: A review of clinical development from bench to bedside. J. Biomed. Opthalmol. 2007, 12, 051403–051424. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Spaide, R.F.; Klancnik, J.M.; Cooney, M.J. Retinal vascular layers imaged by fluorescein angiography and optical coherence tomography angiography. JAMA Ophthalmol. 2015, 133, 45–50. [Google Scholar] [CrossRef] [PubMed]
  7. De Rosa, R.; Vasa-Nicotera, M.; Leistner, D.M.; Reis, S.M.; Thome, C.E.; Boeckel, J.-N.; Fichtlscherer, S.; Zeiher, A.M. Coronary atherosclerotic plaque characteristics and cardiovascular risk factors―insights from an optical coherence tomography study. Circ. J. 2017, 81, 1165–1173. [Google Scholar] [CrossRef] [Green Version]
  8. Hariri, L.P.; Adams, D.C.; Wain, J.C.; Lanuti, M.; Muniappan, A.; Sharma, A.; Colby, T.V.; Mino-Kenudson, M.; Mark, E.J.; Kradin, R.L. Endobronchial optical coherence tomography for low-risk microscopic assessment and diagnosis of idiopathic pulmonary fibrosis in vivo. Am. J. Respir. Crit. Care Med. 2018, 197, 949–952. [Google Scholar] [CrossRef]
  9. Wang, K.C.; Astell, C.J.; Wijesinghe, P.; Larcombe, A.N.; Pinniger, G.J.; Zosky, G.R.; Kennedy, B.F.; Berry, L.J.; Sampson, D.D.; James, A.L. Optical coherence tomography-based contact indentation for diaphragm mechanics in a mouse model of transforming growth factor alpha induced lung disease. Sci. Rep. 2017, 7, 1517. [Google Scholar] [CrossRef] [Green Version]
  10. Schuh, S.; Holmes, J.; Ulrich, M.; Themstrup, L.; Jemec, G.B.; De Carvalho, N.; Pellacani, G.; Welzel, J. Imaging blood vessel morphology in skin: Dynamic optical coherence tomography as a novel potential diagnostic tool in dermatology. Dermatol. Ther. 2017, 7, 187–202. [Google Scholar] [CrossRef] [Green Version]
  11. Wijesinghe, R.E.; Park, K.; Jung, Y.; Kim, P.; Jeon, M.; Kim, J. Industrial resin inspection for display production using automated fluid-inspection based on multimodal optical detection techniques. Opt. Las. Eng. 2017, 96, 75–82. [Google Scholar] [CrossRef]
  12. Cho, N.H.; Park, K.; Kim, J.-Y.; Jung, Y.; Kim, J. Quantitative assessment of touch-screen panel by nondestructive inspection with three-dimensional real-time display optical coherence tomography. Opt. Las. Eng. 2015, 68, 50–57. [Google Scholar] [CrossRef]
  13. Shirazi, M.F.; Wijesinghe, R.E.; Ravichandran, N.K.; Kim, P.; Jeon, M.; Kim, J. Quality assessment of the optical thin films using line field spectral domain optical coherence tomography. Opt. Las. Eng. 2018, 110, 47–53. [Google Scholar] [CrossRef]
  14. Wijesinghe, R.E.; Lee, S.-Y.; Ravichandran, N.K.; Han, S.; Jeong, H.; Han, Y.; Jung, H.-Y.; Kim, P.; Jeon, M.; Kim, J. Optical coherence tomography-integrated, wearable (backpack-type), compact diagnostic imaging modality for in situ leaf quality assessment. Appl. Opt. 2017, 56, D108–D114. [Google Scholar] [CrossRef]
  15. Wijesinghe, R.E.; Lee, S.-Y.; Ravichandran, N.K.; Shirazi, M.F.; Kim, P.; Jung, H.-Y.; Jeon, M.; Kim, J. Biophotonic approach for the characterization of initial bitter-rot progression on apple specimens using optical coherence tomography assessments. Sci. Rep. 2018, 8, 1–10. [Google Scholar] [CrossRef] [PubMed]
  16. Wijesinghe, R.E.H.; Lee, S.-Y.; Kim, P.; Jung, H.-Y.; Jeon, M.; Kim, J. Optical sensing method to analyze germination rate of capsicum annum seeds treated with growth-promoting chemical compounds using optical coherence tomography. J. Biomed. Opt. 2017, 22, 091502–091510. [Google Scholar] [CrossRef] [Green Version]
  17. Lan, G.; Li, G. Design of a k-space spectrometer for ultra-broad waveband spectral domain optical coherence tomography. Sci. Rep. 2017, 7, 1–8. [Google Scholar] [CrossRef]
  18. Zhang, N.; Huo, T.; Wang, C.; Chen, T.; Zheng, J.-g.; Xue, P. Compressed sensing with linear-in-wavenumber sampling in spectral-domain optical coherence tomography. Opt. Lett. 2012, 37, 3075–3077. [Google Scholar] [CrossRef]
  19. Yun, S.; Tearney, G.; Bouma, B.; Park, B.; de Boer, J.F. High-speed spectral-domain optical coherence tomography at 1.3 µm wavelength. Opt. Express 2003, 11, 3598–3604. [Google Scholar] [CrossRef]
  20. Van der Jeught, S.; Bradu, A.; Podoleanu, A.G. Real-time resampling in fourier domain optical coherence tomography using a graphics processing unit. J. Biomed. Opt. 2010, 15, 030511–030514. [Google Scholar] [CrossRef] [Green Version]
  21. Agrawal, A.; Pfefer, T.J.; Gilani, N.; Drezek, R. Three-dimensional characterization of optical coherence tomography point spread functions with a nanoparticle-embedded phantom. Opt. Lett. 2010, 35, 2269–2271. [Google Scholar] [CrossRef] [PubMed]
  22. Zhang, K.; Kang, J.U. Real-time 4d signal processing and visualization using graphics processing unit on a regular nonlinear-k fourier-domain oct system. Opt. Express 2010, 18, 11772–11784. [Google Scholar] [CrossRef] [PubMed]
  23. Hu, Z.; Rollins, A.M. Fourier domain optical coherence tomography with a linear-in-wavenumber spectrometer. Opt. Lett. 2007, 32, 3525–3527. [Google Scholar] [CrossRef] [PubMed]
  24. Wang, K.; Ding, Z. Spectral calibration in spectral domain optical coherence tomography. Chin. Opt. Lett. 2008, 6, 902–904. [Google Scholar] [CrossRef] [Green Version]
  25. Gelikonov, V.; Gelikonov, G.; Shilyagin, P. Linear-wavenumber spectrometer for high-speed spectral-domain optical coherence tomography. Opt. Spectrosc. 2009, 106, 459–465. [Google Scholar] [CrossRef]
  26. Xu, D.; Huang, Y.; Kang, J.U. Compressive sensing with dispersion compensation on non-linear wavenumber sampled spectral domain optical coherence tomography. Biomed. Opt. Express 2013, 4, 1519–1532. [Google Scholar] [CrossRef]
  27. Sugita, M.; Brown, R.A.; Popov, I.; Vitkin, A. K-distribution three-dimensional mapping of biological tissues in optical coherence tomography. J. Biomed. 2018, 11, 1–18. [Google Scholar] [CrossRef]
  28. Uribe-Patarroyo, N.; Kassani, S.H.; Villiger, M.; Bouma, B.E. Robust wavenumber and dispersion calibration for fourier-domain optical coherence tomography. Opt. Express 2018, 26, 9081–9094. [Google Scholar] [CrossRef]
  29. Yuan, Q.; Wang, Z.; Song, L.; Ge, C.; Lu, Z.; Yang, T. Ultrafast wavenumber linear-step-swept source based on synchronous lightwave synthesized frequency sweeper. IEEE Photonics J. 2019, 11, 1–9. [Google Scholar] [CrossRef]
  30. Chan, K.K.; Tang, S. High-speed spectral domain optical coherence tomography using non-uniform fast fourier transform. Biomed. Opt. Express 2010, 1, 1309–1319. [Google Scholar] [CrossRef] [Green Version]
  31. Eigenwillig, C.M.; Biedermann, B.R.; Palte, G.; Huber, R. K-space linear fourier domain mode locked laser and applications for optical coherence tomography. Opt. Express 2008, 16, 8916–8937. [Google Scholar] [CrossRef] [PubMed]
  32. Wang, Z.; Yuan, Z.; Wang, H.; Pan, Y. Increasing the imaging depth of spectral-domain oct by using interpixel shift technique. Opt. Express 2006, 14, 7014–7023. [Google Scholar] [CrossRef] [PubMed]
  33. Zhang, X.; Gao, W.; Bian, H.; Chen, C.; Liao, J. Self-spectral calibration for spectral domain optical coherence tomography. Opt. Eng. 2013, 52, 063603–063611. [Google Scholar] [CrossRef] [Green Version]
  34. Jeon, M.; Kim, J.; Jung, U.; Lee, C.; Jung, W.; Boppart, S.A. Full-range k-domain linearization in spectral-domain optical coherence tomography. Appl. Opt. 2011, 50, 1158–1163. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Han, S.; Kwon, O.; Wijesinghe, R.; Kim, P.; Jung, U.; Song, J.; Lee, C.; Jeon, M.; Kim, J. Numerical-sampling-functionalized real-time index regulation for direct k-domain calibration in spectral domain optical coherence tomography. Electronics 2018, 7, 182. [Google Scholar] [CrossRef] [Green Version]
  36. Hagen, N.; Tkaczyk, T.S. Compound prism design principles, iii: Linear-in-wavenumber and optical coherence tomography prisms. Appl. Opt. 2011, 50, 5023–5030. [Google Scholar] [CrossRef] [Green Version]
  37. Attendu, X.; Ruis, R.M.; Boudoux, C.; van Leeuwen, T.G.; Faber, D.J. Simple and robust calibration procedure for k-linearization and dispersion compensation in optical coherence tomography. J. Biomed. Opt. 2019, 24, 056001–056013. [Google Scholar] [CrossRef]
  38. Makita, S.; Fabritius, T.; Miura, M.; Yasuno, Y. 1st Canterbury Workshop on Optical Coherence Tomography and Adaptive Optics. In Full-Range, High-Speed, High-Resolution 1 μm Spectral-Domain Optical Coherence Tomography with bm-Scan Method for the Human Posterior Eye Imaging; International Society for Optics and Photonics: Canterbury, UK, 2008; p. 713914. [Google Scholar]
  39. Podoleanu, A.G.; Bradu, A. Master–slave interferometry for parallel spectral domain interferometry sensing and versatile 3d optical coherence tomography. Opt. Express 2013, 21, 19324–19338. [Google Scholar] [CrossRef]
  40. Bradu, A.; Israelsen, N.M.; Maria, M.; Marques, M.J.; Rivet, S.; Feuchter, T.; Bang, O.; Podoleanu, A. Recovering distance information in spectral domain interferometry. Sci. Rep. 2018, 8, 1–16. [Google Scholar] [CrossRef]
  41. Cui, D.; Liu, X.; Zhang, J.; Yu, X.; Ding, S.; Luo, Y.; Gu, J.; Shum, P.; Liu, L. Dual spectrometer system with spectral compounding for 1-μm optical coherence tomography in vivo. Opt. Lett. 2014, 39, 6727–6730. [Google Scholar] [CrossRef]
  42. Choi, W.; Baumann, B.; Swanson, E.A.; Fujimoto, J.G. Extracting and compensating dispersion mismatch in ultrahigh-resolution fourier domain oct imaging of the retina. Opt. Express 2012, 20, 25357–25368. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Wijesinghe, R.E.; Park, K.; Kim, P.; Oh, J.; Kim, S.-W.; Kim, K.; Kim, B.-M.; Jeon, M.; Kim, J. Optically deviated focusing method based high-speed sd-oct for in vivo retinal clinical applications. Opt. Rev. 2015, 23, 1–9. [Google Scholar] [CrossRef]
  44. Jian, Y.; Wong, K.; Sarunic, M.V. Graphics processing unit accelerated optical coherence tomography processing at megahertz axial scan rate and high resolution video rate volumetric rendering. J. Biomed. Opt. 2013, 18, 026002–026007. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Choi, D.; Lee, J.; Jeon, M.; Kim, J. In vivo fascicle bifurcation imaging of rat sciatic nerve using swept-source optical coherence tomography. IEEE Access 2018, 6, 7713–7718. [Google Scholar] [CrossRef]
Figure 1. System configuration of the SD-OCT system. Abbreviations: BLS, broadband laser source; C, collimator; DG, diffraction grating; DC, dispersion compensator; E, eye; FC, fiber coupler; GM, Galvano mirror; L, lens; LSC, line scan camera; M, mirror; PC, polarization controller; RSL, retinal scanning lens.
Figure 1. System configuration of the SD-OCT system. Abbreviations: BLS, broadband laser source; C, collimator; DG, diffraction grating; DC, dispersion compensator; E, eye; FC, fiber coupler; GM, Galvano mirror; L, lens; LSC, line scan camera; M, mirror; PC, polarization controller; RSL, retinal scanning lens.
Sensors 20 02067 g001
Figure 2. Graphical representation of the fundamental steps of the software-based intensity detection and resampling algorithm developed for the spectral calibration.
Figure 2. Graphical representation of the fundamental steps of the software-based intensity detection and resampling algorithm developed for the spectral calibration.
Sensors 20 02067 g002
Figure 3. Flow diagram of the fundamental steps performed for dispersion of the system. The red color solid arrows illustrate switching image states (changing of images), while compensating dispersion using decoupled hardware and software methods. The blue color solid arrows illustrate switching image states, while changing the focusing position.
Figure 3. Flow diagram of the fundamental steps performed for dispersion of the system. The red color solid arrows illustrate switching image states (changing of images), while compensating dispersion using decoupled hardware and software methods. The blue color solid arrows illustrate switching image states, while changing the focusing position.
Sensors 20 02067 g003
Figure 4. Evaluations of the point spread function (PSF)-based spectral interferogram analysis. (a) The depth information table of each PSF position, showing the respective number of points carrying maximum intensity, (b) correlated pixel numbers (pixel positions) at each number of points with maximum intensity prior to interpolation, and (c) the averaged graphical plot of the eight respective graphs in Figure 4b.
Figure 4. Evaluations of the point spread function (PSF)-based spectral interferogram analysis. (a) The depth information table of each PSF position, showing the respective number of points carrying maximum intensity, (b) correlated pixel numbers (pixel positions) at each number of points with maximum intensity prior to interpolation, and (c) the averaged graphical plot of the eight respective graphs in Figure 4b.
Sensors 20 02067 g004
Figure 5. Quantitative representation of system sensitivity assessments. (a) Depth-dependent system intensity fall-off and (b) depth-dependent sensitivity roll-off (signal-to-noise ratio, SNR) representation. Each PSF peak of Figure 4a corresponds to the correlated pixel numbers of Figure 4b.
Figure 5. Quantitative representation of system sensitivity assessments. (a) Depth-dependent system intensity fall-off and (b) depth-dependent sensitivity roll-off (signal-to-noise ratio, SNR) representation. Each PSF peak of Figure 4a corresponds to the correlated pixel numbers of Figure 4b.
Sensors 20 02067 g005
Figure 6. Full width at half maximum (FWHM) comparison between before (without) and after (with) spectral calibrated depth-dependent PSFs.
Figure 6. Full width at half maximum (FWHM) comparison between before (without) and after (with) spectral calibrated depth-dependent PSFs.
Sensors 20 02067 g006
Figure 7. Spectrum calibrated two dimensional OCT images acquired prior to dispersion compensation. (a) 2D-OCT images of the infrared (IR)-detection card at top and bottom depth levels. (b) 2D-OCT images of five glass cover slips at top and bottom depth levels. Red color square boxes represent a region of interest with non-corrected dispersion. Horizontal scale bar: 1 mm, vertical scale bar: 300 μm.
Figure 7. Spectrum calibrated two dimensional OCT images acquired prior to dispersion compensation. (a) 2D-OCT images of the infrared (IR)-detection card at top and bottom depth levels. (b) 2D-OCT images of five glass cover slips at top and bottom depth levels. Red color square boxes represent a region of interest with non-corrected dispersion. Horizontal scale bar: 1 mm, vertical scale bar: 300 μm.
Sensors 20 02067 g007
Figure 8. Spectrum calibrated two dimensional OCT images acquired after performing dispersion compensation. (a) 2D-OCT images of IR-detection card at top and bottom depth levels. (b) 2D-OCT images of five glass cover slips at top and bottom depth levels. Green color square boxes represent a region of interest with corrected dispersion. Horizontal scale bar: 1 mm, vertical scale bar: 300 μm.
Figure 8. Spectrum calibrated two dimensional OCT images acquired after performing dispersion compensation. (a) 2D-OCT images of IR-detection card at top and bottom depth levels. (b) 2D-OCT images of five glass cover slips at top and bottom depth levels. Green color square boxes represent a region of interest with corrected dispersion. Horizontal scale bar: 1 mm, vertical scale bar: 300 μm.
Sensors 20 02067 g008
Figure 9. Image comparison with and without spectral calibrated in vivo human retina. (a) In vivo human retinal image before (without) spectral calibration; (b) shows magnified view of red box. (c) In vivo human retinal image after (with) spectral calibration; (d) shows magnified view of green box. (a) Horizontal scale bar: 1 mm, vertical scale bar: 300 μm, (b) horizontal scale bar: 500 μm, vertical scale bar: 150 μm.
Figure 9. Image comparison with and without spectral calibrated in vivo human retina. (a) In vivo human retinal image before (without) spectral calibration; (b) shows magnified view of red box. (c) In vivo human retinal image after (with) spectral calibration; (d) shows magnified view of green box. (a) Horizontal scale bar: 1 mm, vertical scale bar: 300 μm, (b) horizontal scale bar: 500 μm, vertical scale bar: 150 μm.
Sensors 20 02067 g009
Figure 10. OCT image verification of in vivo rat posterior chamber along with B-scans corresponding to multiple positions. (a) Before calibration, (b) after calibration, (c) averaged posterior vitreous opacifications. The horizontal and vertical scale bars of P1–P4 (position 1–position 4) are 1 mm and 700 μm, respectively.
Figure 10. OCT image verification of in vivo rat posterior chamber along with B-scans corresponding to multiple positions. (a) Before calibration, (b) after calibration, (c) averaged posterior vitreous opacifications. The horizontal and vertical scale bars of P1–P4 (position 1–position 4) are 1 mm and 700 μm, respectively.
Sensors 20 02067 g010
Figure 11. The OCT image verification of in vivo rat anterior vitreous opacification along with B-scans corresponds to multiple positions. (a) 3D enface visualization of anterior vitreous opacification, (b) before calibration, and (c) after calibration. The horizontal and vertical scale bars of (b) are 500 μm and 300 μm, respectively.
Figure 11. The OCT image verification of in vivo rat anterior vitreous opacification along with B-scans corresponds to multiple positions. (a) 3D enface visualization of anterior vitreous opacification, (b) before calibration, and (c) after calibration. The horizontal and vertical scale bars of (b) are 500 μm and 300 μm, respectively.
Sensors 20 02067 g011

Share and Cite

MDPI and ACS Style

Han, S.; Wijesinghe, R.E.; Jeon, D.; Han, Y.; Lee, J.; Lee, J.; Jo, H.; Lee, D.-E.; Jeon, M.; Kim, J. Optical Interferometric Fringe Pattern-Incorporated Spectrum Calibration Technique for Enhanced Sensitivity of Spectral Domain Optical Coherence Tomography. Sensors 2020, 20, 2067. https://doi.org/10.3390/s20072067

AMA Style

Han S, Wijesinghe RE, Jeon D, Han Y, Lee J, Lee J, Jo H, Lee D-E, Jeon M, Kim J. Optical Interferometric Fringe Pattern-Incorporated Spectrum Calibration Technique for Enhanced Sensitivity of Spectral Domain Optical Coherence Tomography. Sensors. 2020; 20(7):2067. https://doi.org/10.3390/s20072067

Chicago/Turabian Style

Han, Sangyeob, Ruchire Eranga Wijesinghe, Deokmin Jeon, Youngmin Han, Jaeyul Lee, Junsoo Lee, Hosung Jo, Dong-Eun Lee, Mansik Jeon, and Jeehyun Kim. 2020. "Optical Interferometric Fringe Pattern-Incorporated Spectrum Calibration Technique for Enhanced Sensitivity of Spectral Domain Optical Coherence Tomography" Sensors 20, no. 7: 2067. https://doi.org/10.3390/s20072067

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop