Optimal spectral reshaping for resolution improvement in optical coherence tomography

We analyze the resolution limit that can be achieved by means of spectral reshaping in optical coherence tomography images and demonstrate that the resolution can be improved by means of modelessly reshaping the source spectrum in postprocessing. We show that the optimal spectrum has a priory surprising “crater-like” shape, providing 0.74 micron axial resolution in free-space. This represents ~50% improvement compared to resolution using the original spectrum of a white light lamp. ©2006 Optical Society of America OCIS codes: (100.2980) Image enhancement; (170.4500) Optical coherence tomography; (350.5730) Resolution References and links 1. D. Huang, E. A. Swanson, C. P. Lin, J. S. Schuman, W. G. Stinson, W. Chang, M. R. Hee, T. Flotte, K. Gregory, C. A. Puliafito, and J. G. Fujimoto, “Optical coherence tomography,” Science, 254, 1178 (1991). 2. Brett E. Bouma, Guillermo J. Tearney, Handbook of optical coherence tomograph, (Marcel Dekker, New York, 2002) 3. B. Bouma, G. Tearney, S. Boppart, M. Hee, M. Brezinski, and J. Fujimoto, “High-resolution optical coherence tomographic imaging using a mode-locked Ti:Al2O3 laser source,” Opt. Lett. 2


Introduction
Optical coherence tomography is a promising non-invasive in vivo imaging technique, which has undergone rapid development since its invention in 1991 [1,2].In an OCT system, the axial resolution is determined by the temporal coherence length l c of a light source.If the spectrum of the light source has a Gaussian shape, l c is proportional to λ 0 2 /Δλ with λ 0 the central wavelength and Δλ the bandwidth of the emission spectrum, respectively [2].
Depending on the choice of a light source, the axial resolution of a typical OCT system varies from ~1 to 30 µm.Several techniques have been proposed to improve the resolution of OCT.One possible strategy is to develop new light sources with shorter coherence lengths.For example, a Kerr-lens mode-locked Ti:sapphire laser [3], Ti:sapphire pumped super-continuum [4,5], and thermal light sources [6,7] were used to obtain 3.7 [3], 0.5 [4], ~1 [5], 0.7 [6], and ~1 µm [7] axial resolution in biological tissue, respectively.Alternatively, resolution can be improved by combining multiple light sources, e.g. after optimizing the power ratios of three LEDs, resolution in OCT images was improved from 12 µm to 7 µm [8].Moreover, by digitally reshaping the source spectra to known modes, such as Gaussian, white, and Hamming windowed shapes, OCT resolution was enhanced from over 10 to a few microns [9,10,11].A modeless spectral reshaping, in which a spectral profile can be changed to any arbitrary shape, was used to effectively reduce sidelobes of in OCT images [12].
In this Letter we analyze the resolution limit that can be achieved by means of the modeless spectral reshaping approach, determine the shape of the optimal spectrum, and use this to improve the resolution of OCT.

Theory and simulation
The modeless spectral reshaping we propose here is based on a spectral domain OCT modality in which the interference signal I(k,r) is given by [2] )} where k=2π/λ is the wave number; n is the refractive index of the sample; z 0 is the path length from the beam splitter to the top surface of the sample, z is the depth of the sample, r is the path difference between the reference mirror and the top surface of the sample; a(z) is the depth-resolved backscattered amplitude of the sample; P(k), Q(k), R(k) and S(k) are the source spectrum, spectral response of the detector, reflection spectrum of the reference mirror, and the scattering spectrum in the sample (dispersion within the sample was neglected), respectively; describes the mutual interference of the light scattered from the sample; (2) The bandwidth of an OCT system is limited by the spectrum of a light source as well as the spectral range of the detector and other components.Assuming that the minimum and maximum detectable wave numbers are k 1 and k 2 , respectively, the sum of the interference signal weighted with a distribution V(k) within the range [k 1 , k 2 ] can be expressed as follows: where is the point spread function (PSF) that determines axial resolution and V(k) can be reshaped in postprocessing.In order to improve the resolution, the width of PSF has to be minimized.
Equation (4) can be rewritten as where ⊗ denotes convolution operation, rect(k 1 , k 2 ) is the rectangle function which equals to 1 within the range [k 1 , k 2 ] and zero elsewhere, and is a high frequency oscillation term whose envelope is a sinc function.If V(k) is a constant (i.e.white light source), the first term in Eq. ( 5) would be a delta function, and PSF would equal to U(Z).The peak of this PSF is narrow; however the sidelobes are large.If V(k) is Gaussian or Hamming windowed, the sidelobes can be significantly compressed, however the peak of the PSF is broadened.In this paper we optimized V(k) so that U(Z) functions with different shifts compensate each other to create a narrower PSF peak than that of white light, at the same time the sidelobes are suppressed to a lower level than that of white light also.
It is hard to get an analytical solution for the optimal spectrum V(k).We discretize the optical path difference (OPD) and wave number as Z i (i=1, 2,…M) and k j (j=1, 2 ,…W), respectively, and rewrite Eq. ( 5) in matrix notations: where f i is the PSF for Z i , v j =V(k j ), and . In order to find the optimal spectrum V, we solved the nested optimization problem: , where FWHM represents the full-width half-maximum.The optimization procedure is shown as Fig. 1.Since the shape of the PSF of a broad-band source is similar to the sinc function, we selected the sinc function with a variable width T, sinc(Z/T), as the target PSF.The initial width T was set to 1/k 1 .In the m-th iteration, the target function F m was chosen as sinc(Z/T m ), the half-optimal spectrum V m was found so that is minimal, and the FWHM of its corresponding PSF, D m , was calculated by Hilbert transform.In order to optimize F, we applied a 0.1% perturbation to T m to obtained another FWHM D m ', and in the following iteration, T m was set as to give a PSF with narrower peak.γ is a proportional coefficient.Normally a smaller value of γ leads to more stable convergence but also increases the complexity of computations.In our computations, the convergence criterion was set as γ(D-D')<10 -5 T m , and γ was chosen as 10 to ensure the convergence after 100-1000 iterations.The reflective Newton iterative method [13], which allows setting the range of v j , was chosen to solve the matrix-form least-squares problem, and set the lower limit of v j to zero to satisfy the physics of real light sources.The system is spectrally sensitive, and the signal-to-noise ratio (SNR) values differ significantly at different wave numbers.In order to limit the use of noisy spectral data, we applied an additional condition α and its FWHM D m ' Increase m by 1 value, the SNR has larger effect while FWHM has lower effect on the optimization procedure.
After a few attempts we selected α as W, the number of the wave numbers.
The obtained optimal spectrum was shown in Fig. 2. Different from the conventional symmetric bell-like reshaping spectrum, the optimal spectrum has a priory surprising asymmetric crater-like shape with two peaks close to the lower and upper bounds of the spectrum, and the peak at shorter wavelength is higher.It can be understood with this way: the resolution is primarily limited by the shortest and longest detectable wavelengths, especially the shortest detectable wavelength, so the side-spectrum enhancement, especially the enhancement at the shorter wavelength side, helps improve the resolution.In Fig. 2 we also show five other spectra: the modified optimal spectrum which is more symmetric with the amplitude of the two peaks being the same, Gaussian profile with FWHM = 300 nm, white light, white light weighted with Hamming window, and the xenon lamp spectrum.To be consistent with our experimental system, a wavelength band from 421 nm to 895 nm was used.The envelopes of theoretical PSFs of the six spectra are shown in Fig. 4(a), and their FWHM values are given in the notation text.Evidently, the optimal spectrum provides the highest FWHM resolution, 0.44 μm in free space, with sidelobes lower than those of the PSF generated by the white light spectrum.The PSF peak of the modified optimal spectrum is slightly broader than that of the optimal spectrum, which proves that an asymmetric spectral profile may lead to a higher resolution than a symmetric spectral profile.To quantitatively compare the sidelobes effect of the PSFs, we calculated the root-mean square width (RMSW) which represents the deviation of the PSF from the central position z=0 [12].The theoretical FWHM and RMSW values are shown in Table 1. the optimal spectrum provides the highest FWHM as well as adequate RMSW.

Experimental results
The approach reported in section 2 was also validated experimentally.Our spectroscopic OCT system setup is shown in Figure 4.A 250W xenon lamp (Oriel) was used as the light source.The output of the xenon lamp was first collimated by two 250-mm-focal-length lenses L1 and L2 together with a 1 mm pinhole, spatially filtered by an adjustable aperture, and then separated into the reference and sample beams by a polarization independent beam splitter BS.We used the largest facet of an uncoated right-angle prism as the reference mirror M1.Two identical 4x achromatic objective lenses were used to focus the light beams to M1 and the sample, and collected the light back reflected/scattered by M1 and the sample.The collected light was analyzed by means of a spectrograph (Acton SP2150i) coupled into a CCD camera (Princeton instruments PIMAX 1KHQ), and magnification of the sample to the CCD was 20.The width of the entrance slit of the spectrograph was adjusted to be 10μm, that is, the vicinity of 10/20=0.5μm of the reference mirror or the sample was regarded as a single dot.light from the .The detection system was sensitive to wavelengths from 421nm to 895 nm, with a spectral resolution of 2 nm.When biological tissue was used as a sample, the CCD was set to a higher gain and longer integration time, and in order to avoid the CCD saturation, a neutral density filter with 10 dB single-trip loss was inserted in the reference arm.A glass was placed in the sample arm to compensate the additional optical path introduced by the neutral density filter.Reference mirror M1 and lens L3 were moved as one unit to change the OPD between the reference and sample beams with a step size of 0.1 μm.Interference spectrum I(λ, r) was recorded at each OPD step.For calibration, I(λ, r) reflected by the reference mirror and the sample were also recorded at both the beginning and the end of each imaging experiment.We then followed Eq. ( 2) to obtain I d (k,r) for different k and r values, and applied different spectra V(k) to get the corresponding I d (r) values according to Eq. (3).A piece of microscope cover glass was used as the sample first.The light was focused on the top surface of the glass and a(z) can be regarded as a delta function.Thus the resulted I d (r) gave the PSFs for different source spectra V(k).Fig. 3(b) shows the envelope of the experimental PSFs.There are multi-peaks in the theoretical PSF, with about 100nm between the adjacent peaks.In the experiments, the un-flatness of components surfaces and the dispersion induced optical path mismatch, which are in the order of 100nm, smooth the peaks.So there is no multiple peaks in the experimental PSFs.The corresponding FWHM and RMSW values are given in Table 1.As shown in Fig. 3(b) and Table 1, the modeless spectral reshaping indeed provides resolution superior to those obtained with other conventionally used spectra including the original xenon-lamp spectrum, white light, Gaussian spectrum, and Hamming windowed profile, with the FWHM reduced by 51%, 20%, 28%, and 30%, respectively.Although the spectrum was not purposely reshaped to optimize the sidelobes, we also notice that the optimal spectrum suppress sidelobes compared to the original xenon-lamp spectrum with RMSW reduced by 54%.The sidelobes are larger, however, than those of corresponding to the Gaussian and Hamming profiles.
As a further illustration of the improved resolution, some onion root tip tissue was used the sample.The experiment was repeated by 100 times, and the results were averaged to give a higher SNR.Figs.5(a, b) show the B-scan images before and after spectral reshaping, respectively, and Figs.5(c, d) show the corresponding A-scan traces taken from the images at 75 μm lateral position.Evidently, the image/trace obtained using the optimal spectrum (Figs.5(b, d)) has an improved higher resolution compared to those obtained using the original xenon-lamp spectrum (Figs.5(a, c)).

Conclusions
In summary, we showed that the optimal spectrum has a priory surprising crater-like spectral shape, which is different from the conventionally adopted bell-shape spectra.Although the resolution provided by the optimal spectrum is primarily determined by the spectral range, we further demonstrated, for the first time to our knowledge, the optimal spectral reshaping results in an extra improvement in resolution.In particular, in our experiments, the resolution of 0.74 μm in free space was obtained, which is 51% higher than the one provided by the original xenon-lamp spectrum.Therefore, the modeless spectral reshaping technique provides a flexible way to improve the resolution of OCT and may help achieve sub-micron resolution for subcellular OCT imaging.
signal-to-noise ratio at wave number w v which was obtained by dividing the signal level at at wave number w v by the amplitude of background white noise.α is a threshold controlling the trade-off of the SNR and the FWHM.For a less α

FFig. 3 .
Fig. 3.The optimal PSF provides improved OCT resolution for the source spectra shown in Fig. 3. (a) Envelopes of theoretical PSFs (b) Envelopes of experimentally measured PSFs.PSFs in the insets are normalized to conserve the L2 norm.

Table 1 .
Optimal PSF has improved FWHM and RMSW compared to other PSF given by white light, Gaussian, and Hamming windowed spectra