基于全空间调制谱分割角度复合的微血管造影方法与系统的制作方法
【技术领域】
[0001 ]本发明设及光学相干层析成像技术(Optical Coherence Tomogra地y，0CT)W及 在此基础之上的无标记S维光学微血管造影技术(OCT Angiogra地y,Angio-OCT)，尤其设 及基于全空间调制谱分割角度复合的微血管造影方法与系统。 技术背景
[0002] OCT技术是自从上世纪90年代W后逐步发展起来的在生物医学领域具有重大应用 价值的一种新型成像技术。凭借着无标记物、非侵入性、非接触性、高成像分辨率和高探测 灵敏度等优点，它吸引着越来越多的国内外科研工作者更加深入地对其进行研究，目前在 医学临床中得到了广泛应用。OCT技术主要依据所探测到的由于生物样品的光学不均匀性 所导致的样品背向散射光光强变化，获取样品内部的反射率信息，进而重构出样品的断层 结构图像。然而通常在疾病的早期阶段，由于正常的和病变的生物组织间光散射特性的差 异性较小，运用传统的结构型OCT较难加 W检测并区分，故而在临床应用上其存在诸多的局 限性。为了能够获得除生物样品组织形态结构之外更多的生理信息，功能型OCT技术得到了 开发并广泛应用。
[0003] 诸多的疾病（如青光眼、老年黄斑病变等眼底疾病，脑中风等脑疾病）与血流灌注 的病态变化存在密切的相关性。如若能够实时监测病人血管的变化，并提供血管的=维影 像，对该类疾病的早期诊断和控制具有重要的意义。Angio-OCT技术是一项具有巨大前景的 可有助于血管疾病诊断的工具，其能够较好的区分静态的组织和动态的血流信号，实现了 OCT在获取血管内血流信息方面的功能性拓展。相较于传统的造影技术，其优势体现在无需 注入显影剂和利用X射线，并结合前面所述OCT技术的特点，能够实现微血管S维深度分辨 的高对比度的成像。通常，血流对比度的模型建立是运用数学的手段分析光散射信号的时 间统计特性，并用阔值分割动态的血流信号和静态的组织背景信号来实现的。然而由于动 态和静态信号统计曲线之间存在重叠，导致动态和静态信号分割错误，血流对比度受到抑 审IJ。因此，采取有效的方法提高造影图的血流对比，有助于更加清晰地解释说明图像的特 征。
[0004] 基于Angio-OCT信号的时间统计特性研究得出：多个独立的血管造影子图的平均 可W提高造影的对比度。参照一些消除散斑的方法，可W通过波长多样性、角度多样性和偏 振多样性等方法来获取独立的造影子图。Jia等人提出了一种类似于波长多样性的光谱分 割Angio-OCT方法。该方法将OCT干设信号的全波长光谱分割成不同的子光谱，每个子光谱 可产生独立的造影子图，通过复合形成新的血管造影图。然而相对于原来的全光谱，每个子 光谱带宽变窄，导致图像的轴向分辨率下降。
【发明内容】
[000引本发明针对现有技术的不足，提出了基于全空间调制谱分割角度复合的微血管造 影方法与系统。
[0006] 基于全空间调制谱分割角度复合的微血管造影方法，具体包括W下步骤：
[0007] 1)通过在OCT深度(Z)域去除共辆镜像，重构复数值的干设光谱，在空间频率域得 到全空间的横向扫描调制谱；
[0008] 2)利用空间频率域的横向扫描调制谱编码不同入射角度的探测光，通过分割调制 谱获得入射角度分辨的相互独立的OCT干设子光谱；
[0009] 3)结合2)中的调制谱分割技术与光学微血管造影技术，得到基于横向扫描调制谱 分割的多个角度分辨的独立的血管造影子图；
[0010] 4)结合1)中的全空间调制谱和2)中的调制谱分割技术和与光学微血管造影技术， 获得基于全空间横向扫描调制谱分割的多个角度分辨的独立的血管造影子图；
[0011] 5)对3)或4)中得到的血管造影子图进行平均，得到多个空间角度复合的血管造影 图。
[0012] 根据2)中的调制谱分割技术，实现基于横向扫描调制谱分割的多个角度分辨的血 流绝对流速测量；
[0013] 根据1)中的全空间调制谱和2)中的调制谱分割技术，实现基于全空间横向扫描调 制谱分割的多个角度分辨的血流绝对流速测量。
[0014] 本发明的目的是通过如下技术方案实现的：
[0015] 基于全空间调制谱分割角度复合的微血管造影方法：
[0016] 1)在OCT成像中，对应于一次横向扫描得到的干设光谱实数信号，可表示为S(k， X)，其中k代表波数，X代表横向的快扫方向。由于直接对光谱实数信号沿横向X方向作傅里 叶变换，所得到的空间频率分布会产生正负频率的混叠，因此需要重构出复数值的OCT干设 光谱。为了重构出复数谱，对实数谱S化，X)沿着k方向进行傅里叶变换，得到深度域(Z)的空 间结构信息。通常在谱域OCT中，一半的Z空间对应约3mm的成像范围，能够满足绝大多数的 应用要求，于是将待测样品完全置于零光程差的一边，W分辨复共辆镜像。通过去除Z域中 的共辆项，沿k方向进行傅里叶逆变换，最终得到复数值的光谱信号对该复数谱沿X 方向进行傅立叶变换，在空间频率域获得全空间的横向扫描调制谱。
[0017] 2)通常在OCT系统的样品臂中，为获得较高的横向成像分辨率，要求照射到物镜上 的光束直径尺寸尽量大。因此，当一定宽度的准直光束照射到样品臂扫描振镜的转轴中屯、 时，光束距离转轴中屯、存在一定的偏移量8。而偏移量S会引入样品臂光程的调制，导致在空 间频率域的横向扫描调制，其调制频率fm与偏移量S成线性函数关系：
[0019] 其中，k表示光源的中屯、波数，CO表示扫描振镜的角速度。因此，在全空间的横向扫 描调制谱中，不同的调制频率对应着不同入射角度的探测光。对该调制谱进行分割，可获得 入射角度分辨的相互独立的OCT干设子光谱。分割得到的子光谱个数及子光谱间的重叠视 最终的图像效果而定:分割的个数越多会导致横向分辨率减小得越多，此外还增加了处理 算法的复杂度;子光谱间相互重叠的部分越大，会导致子光谱相互之间的独立性下降，影响 编码不同入射角度的探测光。
[0020] 3)对获得的入射角度分辨的相互独立的子光谱，利用目前的微血管造影技术，如 幅值差分法、复数差分法、去互相关法和散斑方差法等，分别产生相对应的独立的造影子 图。对运些空间角度分辨的造影子图进行复合平均，得到新的微血管造影图。
[0021] 传统的多普勒频移测流速法只能测量与探测光束方向平行的速度分量。要测量血 流的绝对速度，还需测得探测光束方向与血流方向的夹角（即多普勒角）。基于全空间调制 谱分割，能够分辨不同入射角度的探测光束。对分割得到的角度分辨的调制子光谱，利用多 普勒血流速度计算方法，可W获得不同入射角度探测方向的流速分量，进而根据几何关系 可确定绝对的血流速度值。
[0022] 基于全空间调制谱分割角度复合的微血管造影系统，包括低相干宽带光源、光环 行器、光纤禪合器、参考臂、样品臂、光谱仪和信号处理模块;其中参考臂包括参考臂偏振控 制器件、参考臂准直透镜、参考臂聚焦透镜和平面反射镜;样品臂包括样品臂偏振控制器、 样品臂准直透镜、正交扫描振镜和样品臂聚焦物镜;光谱仪包括光谱仪准直透镜、闪耀光 栅、傅里叶透镜和线阵CMOS相机；
[0023] 低相干宽带光源经过光环行器与光纤禪合器一侧的输入端相连接;光纤禪合器另 一侧的其中一个输出端口经过参考臂偏振控制器件与参考臂准直透镜相连接;参考臂聚焦 透镜的光轴与参考臂准直透镜的光轴重合，平面反射镜置于参考臂聚焦透镜的焦平面处； 光纤禪合器另一侧的另外一个输出端口，经过样品臂偏振控制器与样品臂准直透镜相连 接;正交扫描振镜的第一扫描镜转轴中屯、位于样品臂准直透镜的光轴上，第二扫描镜转轴 中屯、位于样品臂聚焦物镜的光轴上，待测样品置于样品臂聚焦物镜的焦平面处;光环行器 的另一个输出端口与光谱仪准直透镜相连接，闪耀光栅根据分光原理放置在光谱仪准直透 镜的出射光路上，傅里叶透镜置于闪耀光栅的出射光路上，线阵CMOS相机的采光面与傅里 叶透镜的后焦平面重合;线阵CMOS相机后面接信号处理模块；
[0024] 低相干宽带光源发出的光经过光环行器后入射到光纤禪合器，出射光分成两部 分，一部分进入参考臂，经准直聚焦后照射到平面反射镜；另一部分进入样品臂，经过准直 聚焦后照射到待测样品上；样品臂的正交扫描振镜实现样品臂光束对待测样品的=维扫 描;参考臂反射镜反射回的光与待测样品背向散射回的光在光纤禪合器处发生干设，出射 的干设光经过光谱仪后由信号处理模块采集处理。
[0025] 与【背景技术】相比，本发明具有的有益效果是：
[0026] 1)利用本发明获得的多个空间角度复合的微血管造影图，具有增强的对比度和血 管连通性
[0027] 2)相比于【背景技术】中谈及的类似于波长多样性的光谱分割Angio-OCT方法，其对 光谱分割会导致轴向分辨率的下降。本发明对轴向分辨率不产生影响。
[0028] 3)本发明谈及的全空间的横向调制谱分割技术，可W提供最大化的空间频率域的 横向调制谱，消除由于调制谱分割导致的横向分辨率下降。
[0029] 4)利用本发明谈及的基于全空间调制谱分割的多个角度分辨的探测技术，可W在 单光束单次测量的情况下实现血流绝对流速的测量。
【附图说明】
[0030] 图1为样品臂扫描光束示意图；
[0031 ]图2为偏移量S引入的调制频率fm图；
[0032] 图3为横向扫描调制谱正负频率混叠图；
[0033] 图4是本发明的方法实现流程图；
[0034] 图5是本发明的成像系统示意图。
[0035] 图1中:①、探测光束;②、扫描振镜;③、振镜转轴中屯、;④、聚焦物镜。
[0036] 图2中：偏移量S引入的调制频率fm。
[0037] 图3中：正负频率混叠的横向扫描调制谱。
[0038] 图4中:A表示OCT干设光谱实数信号;B表示深度(Z)域的空间结