基于界面信号的扫频光学相干层析成像系统相位矫正方法

由于扫频光源的采集触发信号和采样时钟信号存在时间上的随机延时,导致扫频光学相干层析成像(SS-OCT)系统干涉信号光谱的整体错移,进而引发OCT空间域信号的相位跳变,阻碍了基于相位信息的功能成像.为了获得稳定的相位,便于开展功能OCT的研究,提出了一种基于界面信号的数字相位矫正方法.对界面附近相邻A-line间同一深度的相位信号进行差分运算,计算得到相位跳变的A-line位置与光谱错移量(以像素为单位),然后在原始干涉信号上对齐光谱,重新傅里叶逆变换,得到矫正后的复信号.该数字矫正算法不会引入额外的相位噪声,可以实现OCT信噪比受限的相位探测.通过对反射镜、荧光板和小鼠脑皮层血流的多普勒成像验证了该方法的可行性.

Phase correction metho d based on interfacial signal in swept source optical coherence tomography

There are intrinsic phase errors in swept source optical coherence tomography (SS-OCT), which severely influences the functional imaging. To overcome this difficulty, a numerical correction method is presented in this paper to correct the phase artifacts due to wavenumber shift among the spectral interferograms, resulting from the random delay variance between the sampling trigger and the clock of the swept source laser. This correction method is based on the linear relationship of phase difference to the depth of the sample and the wavenumber shift. The detailed procedure to eliminate the phase artifacts is as follows. Firstly, we figure out the complex OCT signals through inverse Fourier transform of the initial interferograms. Then we fit the upper surface of the sample with the intensity information of the B-scan by setting a floating threshold. After that the wavenumber shifts of each A-line are determined by two steps with the phase information of the sample surface: the relative wavenumber shifts between adjacent A-lines are first obtained according to the phase difference between the adjacent A-lines, the signal depth, and the linear relationship mentioned above; then we figure out the absolute wavenumber shifts between each A-line and the first A-line of the B-scan by an iteration algorithm. With the information about the wavenumber shift, we align the initial interferograms, and obtain the corrected complex signal through re-inverse Fourier transform of the aligned interferograms. This method introduces no extra noise, realizing phase measurement limited by the signal-to-noise ratio. It is noted that we take the average phase information of several axial pixels near the sample surface to diminish the noise influence when calculating the wavenumber shifts. Besides, this corrected algorithm acquires oversampling along the scanning direction to ensure the signal correlation between adjacent A-lines.The SS-OCT system in the paper is set up with a vertical cavity surface emitting laser with a center wavelength of 1297 nm. The system measurement range is 12 mm in lateral direction, the axial resolution is 17 μm, and the lateral resolution is 24 μm. And the feasibility of this method is verified by Doppler imaging of a mirror, an infra-red detection card and the cerebral cortex of a mouse.