基于汉明窗算法的扫频OCT纵向分辨率的优化.docVIP

基于汉明窗算法的扫频OCT纵向分辨率的优化 摘要： 针对扫频光源相干层析成像系统（SSOCT）光谱的非高斯型分布，对采样到的干涉光谱实施了基于汉明窗函数的光谱整形。利用选用的汉明窗整形算法，成功地实现了对样品（镀膜的塑料卡片）上界面层析成像纵向分辨率的提升，使其上界面膜层的分辨率从128.9 μm提高至109.4 μm。由于更好地压缩 纵向扫描（Ascanning）信号的展宽，因此基于汉明窗算法的光谱整形的方法比传统对原始数据直接快速傅里叶变换（FFT），能够实现更高的纵向层析分辨率。基于这个选定的算法，通过光谱整形后，软件呈现了样品的实时结构层析图像。 关键词： 光学相干层析技术； 窗函数； 纵向分辨率 中图分类号： TH 741.4文献标识码： Adoi： 10.3969/j.issn2013.05.012 引言 作为一种新兴的高分辨率层析成像技术，光学相干层析技术（optical coherent tomography，OCT）[1]因具有无损，快速实时成像的特点，结合内窥镜技术，它已被广泛应用于各种疾病（如眼科青光眼、牙釉质疾病、心血管等）的诊疗中。由于传统时域OCT（time domain OCT，TDOCT）系统发展到了无法突破成像速度的技术瓶颈，频域OCT（fourier domain OCT，FDOCT）通过对图像数据并行采集的方式，而具有比TDOCT更高的成像速度，灵敏度和动态范围的优势[2]。其中基于扫频光源的光学相干层析成像系统（swept source optical coherence tomogramphy，SSOCT），除了具备频域FDOCT的快速成像能力的同时，还具备时域TDOCT点探测时高分辨率的特性，因此作为新一代FDOCT的代表，而逐步成为OCT领域研究热点。 1原理 1.1系统结构 图1所示为扫频光源光学相干层析成像（SSOCT）系统结构图，系统装置采用扫频速率50 kHz，扫频范围1 250～1 360 nm，中心波长在1 310 nm，平均输出功率为15 mW的扫频激光光源和光纤系统集成。系统中光源输出的近红外光先后通过光纤环形器和耦合器，其中一路耦合进参考臂，另一路通过样品臂。在样品臂中系统采用基于微机电反射镜（MEMS mirror）的内窥探头[3]。MEMS mirror会将格林透镜出射的聚焦光束从探头的侧壁出射到样品内部，在驱动信号的带动下可完成对样品二维XY方向的横向扫描。之后从样品返回的光从原路再次返回到光纤耦合器。 耦入进参考臂中的光先后经过准直镜、聚焦透镜和平面反射镜，再次返回到光纤耦合器中。从耦合器出射的光一路进入平衡探测器的正输入端，另一路经过光环形器后再进入平衡探测器的负输入端。干涉光强信号在平衡探测器内部转换成电信号，就被NI数据采集卡并行采集。同时由扫频光源产生的一路触发信号对数据采集进行触发，将采集到的数据通过PCI总线存储到计算机的硬盘中。通过软件对存储的数据进行基于1 024个采样点的汉明窗光谱整形，最后经过快速傅里叶变换FFT和图像重建得到样品（一张有镀膜的塑料卡）深度方向的灰度图像。 由系统成像函数公式（7）可知，样品深度坐标z和光束波数k是一对傅里叶变换对，为了确保系统轴向分辨率，需对等间隔分布的波数对应的干涉光强进行FFT。由于在离散数字信号系统中，作为时域波形信号的干涉光谱包络会影响频域信号分辨率和旁瓣[4]，又因为扫频光源输出的光谱形状是非高斯型光谱，这直接导致干涉光谱包络的非高斯性，从而间接影响到了图像的质量和系统分辨率。因此通过对时域干涉光谱直接整形的方法来减小旁瓣，减少图像中噪声，以提高图像的质量。其中，一种可行的整形方法便是采用在光谱域中基于汉明窗函数的光谱整形[5]，它是让原始干涉光谱信号，与一个合适的窗口函数相乘，来确保系统的轴向分辨率的提高和信噪比提升。 2仿真与试验 为了检验光谱整形效果，采用镀膜的塑料卡作为样品。图2为光谱整形前后塑料卡的干涉光谱，而图3为整形前后对干涉信号进行FFT之后得到的轴向扫描信号。可以看出，光谱整形前后干涉光谱包络有很大差异，在光谱整形之前，干涉光谱包络呈现非高斯型（见图2（a）），轴向扫描信号展宽严重（见图3（a）），轴向分辨率较低，且峰值出现错位；光谱整形之后，干涉光谱具有近高斯包络（见图2（b）），轴向扫描信号的半峰全宽（FWHM）大大减小（见图3（b）），轴向分辨率得到改善。 测试其轴向扫描信号，光谱校准前后塑料卡在不同光程差处的轴向扫描信号如图3所示。可以看出，光谱校准之前，原始干涉信号进行FFT之后产生的轴向扫描信号出现展宽现象，信号峰值周边的噪声严重，系统分辨率较低；随着光程差加大，展宽情况更加严重，分辨率下降更快。光谱校准之后的轴向扫描信号，信号峰值周边的噪声被